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1. FOOD SERVICE MANAGEMENT (U105) 

1.1 Vrhe, Problems of Food Service Management 

)The food service management for an institution j 
such/as a hospital -must concern itself with the problems 
of serving a large number of mea'ls every day. Such a' 
service demands, at the very least: 'l) that sufficient 
quanti ties'of food items be on hand to prepare the menu 
items, 2) that certain nutrient requirements be met, and 
3) that reasonable costs \\e maintained. There are, of 
course, other aspects $f food service, but for the 
purposes of this unit we wi^l concentrate on these three. 
1.2^ Definitions 

We define food items as purchased^ raw food, and'a 
menu item as a' single serving of a dJh made from food 
fit-ems. /f'oT^xample, cake is a menu item whose ingredie- 
nts consist of food items sucfTas eggs, flour, butter 
and sugar.f Nutrients are the properties of food such as , 
calories, protein, fat, carbohydrates, calcium, vitamin 
^ A. . / 

' 1.3 Calculating Menu Costs 

[ Let us suppose tha\ we wished t^ calculate, the cost Q 

of a m^nu item, for. example, a pound cake. The recipe 
calls ,for. 1 unit, given as a weight, each of eggs/ flour,' 
. sugar, and butter. They cost of these units are, respec- 
tively, 70 , 10, 25, and 50 cents. It i<r easy to see that 
the cost of the ingredients in a pound cake i^ : 

' * ' 7Q(1)_+ 10(1) +JS(1).+ Stfd) = 1SS cents = $1.SS. 

. Another recipe, _ one foO^: rambled eggs , 'requires 

j unit of eggs and J-'unit of butter. This cost, is leal- ' 
culat-ed as: , 

•70(j) + S0.(}) = m\ cents. 
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The cost of each of- these »menu items is obtained by 
^taking the sum of several products. The products, in 
each case, are the costs of a unit of food item multiplied 
by the number of ijnits of that item^which is needed in 
the recipe. Thus is very simple arjfthmevt ic . However, 
since we are dealing with a large ^number of recipe^ with 
a wide range of ingredients, there is some efficiency in 
organizing our calculation in a structured- way . We will 
proceect to demonstrate how Sjiis may be done.' - 

1.4 A Matrix Representation of Menu I.tems 

The ingredients of each menu item can be arranged 
a^columns. of a matrix in the following fashion: 



Menu I tern 



eggs 
flour 
sugaY 
butted ' 
beef 

s troganof f 




beef strog. 

0 " 
0 

' (0 
0 



X 



*Beef ^trogarioff is a convenience food, and is 
purchased already prepared." For that reason it 
1 i*sted as both a food i ten) an(fa menu i tern. ' 



-The eotuqna o£ this "ingredient matrix 'contain the 
quantity of each food item neededfcfor the menu i t tem that 
is represented by each .column . - The rows i pf the matrix 
represent the food interns as they appear in various recipes. 
A zero entry indicates that the food item is' not used in 
^a particular recipe.' Tne list of frod items and menu- f 
items is limited here for simplicity". A practical list ) 
woulrd contain hundreds of items. 



K 5 The Price Vector 

A- matrix consisting of a single row or column ij£ 
called a vector. We can arrange the prices of each of 
the food items in the above matrix as such an ordered 
rov of numbers to define a price' vector: 

beef 

eggs *flour sugar butter strog. 
Price =,[70 - 10 25 4 * 50 100 -•/.] 

1.6 Calculating the Cost Vector 

Iff wemultiply *the pr^Qe vector by the ingredient 



)r^ce ve< 
recto rwl 



costs for each recipe 




That 


is, 


for 


the example values 


given: 








/ 




[70 10 25 50 100] 


1 


.5 


1 


(T 




j 


1 


.0 


• 25 


0 






1 


0 


0 


. 0 


( 




1 


.25 


. 25 


0 




_0 


0 


0 


1_ 

4 


= [1SS 47. S 85 100) 
\ 



Notice the order in which we write the left side of this 
equation: row vector first, ingredient matrix second. 

Compare the computation made previously for the first 
two of the menu items in this matrix with the computation 
which is made in matrix multipl ication. We can s"ee that 
the* computations are the same, and that "the resulting 
vector does indeed contain, in order, the cost of each 
menu item. . * 

Some of the advantages of structuring the problem 
this -way -should be evident. For one, matrix multiplica- 
tion accomplished the needed computation. In this form, 
the data c^an be easily entered and* the operations per- 
formed on a computer. Second, a simple change in pVice 
can be entered once and wilP always be "applied to alj 
recipes. Finally, new recipes and ingredients can be > 
added to an^expanded matrix with little trouble. 
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, Note: Sometimes itMs convenient to i?se matrix 
multiplication just to sum a series of columns (or rows) 
of a matrix. For example, suppose we wished to add each 
column of this matrix: 

2.5 10.0 9.4 7.4' 
7.3 7.6 3.6- 2.1 
7.0 3.2 1.9 8.3 

Verify that the addition could be accomplished by 
premultiplying the matrix by a row vector, [1 1 l] . 
This is a convenient way to find the column sums when 
We are using a computer to perform tHe calculations. 
QuIS^tion: If we wished to find the sum of, 
the rows of a matrix by means of 
, .-matrix multiplications', how could 
we do it? 

1 . 7 ♦ Ex ercises , < , > 

\ : — . * * 

1. Suppose the number of servings of each recipe to be prepared on 

siven day for pouVid cake, scrambled e^gs, omelet, and beef 
.strpganoff is, respectively, 5,JO*,*2, and 50- -Using matrix 
multiplication, -calculate the amount of each food^needed on 
that day. The resulting vector is called tbe (input) food 1 
package for the day. 1 . 

2. Find the cost of the food package using the price vector given 
in Section f. 5. " t 

3. Suppose that number of recipes for several days is as follows: 



Day 1 
Day* 2 

Day* 3 



2, ' 
1, 



10, 
3, 
1. 



2, 

3*. 
10, 



50 
6 

0. 



These recipes are' the same as those in Problem 1. Formulate 
these nee<fs as a matrix. Find the total amount of each food 
needed for each of the three days (or find the input food 
package for each day). 

Find the cost/of the food package for each of the three days 
using the^frice vector given in Section 1.5. . ^ 
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1.8^ Sam ple Menu ^T>gms and Costs Per Serving 



Sample Menu 'Items 

Dressing and Gravy # ^ 
Brown Gravy 
Onion Gravy 

Spanish Sauce « 

Tartar Sauce ■ • ^ . 

Bar-B-Q Saufce 

Chicken- Gravy 
- Cole Slaw Vinegar Dressing 

Pimento Cheese * 

Russian Dressing 

French Dressing 

Thousand island Dressing 

Lemon Sauce ► * 

Whipped Topping ^ 

Pickle Salad Dressing ; * 
' ^utrees 

Beef Stew with Vegetables 
Chili With Beans 
Meat Loaf and Gravy 
. Oven-Fried. Steak * 
Roast Bee£ with Gravy 
Country -£tyle Steak 
Smothered Steak \ 
Baked Swiss Steak 
v Baked .Macaroni anfl Cheese 
Baked Haddock * 
Sole Fillet with Tartar Sauce 
Fish Sticks with Tartar Sauce 
Salmon Pattie 
Bxeade4 Pork Chop 
Breaded Pork Cutlet- 
'Pork ChoR with Mexickji Sauce • 
^Deep-Fried Pork Cutlet - 
Rbast Fresh Ham, with Gravy 



\ 



Cost in Cents 
to Management 
- per Serving 

fc.66 

" 0.93 N 
3.75 
1.98 

1. v25 
0.53 
0.36 
4.38 

2. U 
1.29 
1.45 

• 0,.76* 
■ 1.08 



22,95 

18.46 

22.90 

28.16' 

40.93 

48.94 

50.51 

54.52 

10.00 

30.77 

32.00 

12.2*5 

52.49 

56.54 

42.1$ 

48.30 

42.21 

34.29 
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Sample Menu 1 1 ems 
' — 1 — — % — 

Bar-B-Q Chicken 
Fri^d Chicken 
Roast Turkey with Gravy] 
Chicken- Fried Veal Cutlet 
Oven-Fried Veal Cutlet 
Baked Ham Loaf 
Hot Corned Beef ■ 
Smothered Liver with Onions 
Polish Sausage 
Sauteed Chiqken Livers 
Vegetables 

Buttered Whole Kernel Corn 
Seasoned Hominy 
Baked Beans 

Seasoned Blacfcfeyed Peas 
Buttered Green Lima Beans \ 
Buttered' Egg Noodles 
Potatoes Au Gratin 
Baked Potato 
Buttered Diced Potatoes 
Creamed Diced Potatoes 
Hash Browned Potatoe's 
French Fried Potatoes 
Whipped Potatoes ' 
Oven-Browned Potato 
Paprika Diced Potatoes ^ 
Buttered Steamed Rice 
Rice Piiaf 

Glazed Swejet Potatoes 
Buttered Steamed Cabbage 4 
Harvard., Beets 
Cauliflower Au Gratin 
Buttered Broccoli ' 
Buttered Brussels Sprouts 
Buttered Cauliflower 



Cost in Cents 
to" Management 
per Serving 

' 46..21 

44^50 

.38.21 



44.45 
. 42.21 
32.93* 
*34. 50 
34.36 
30.12 
16.80 



§.31 
4.18 
7.90 
7.85 
5.34 
2.23 
4.97 
8.26 
4.16. 
2.47 
8.49 
16.12 
2.51 
8.68 
2.84 
1. 41 
3.06 
8.67 
5.16 
2.61 
6.94 
7.55 
14.29 
10.30 
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Sample Menu Items 

k * * * • 

Buttered Diced Carrots 

Parsley-Biittered Carrots 

Buttered Spinach ' 

Buttered Chopped Turnip Greens 

Buttered French-Cut Green Beans - 

Buttered Canned t Green Beans 

Buttered Mixed Vegetables 

Buttered* Greej* Peas 

Buttered Wax Beans 

Stewed -Tomatoes 

Buttered Onions , 

Seasoned Yellow Squash 

j 

Salads ^ 

k Pineapple Waldorf Salad 
Strawberry Jello with Bananas 
Cabbage Slaw/Green Pepp.ers 
•Mardi Gras Cole Slaw 
Carrot-Celery Sticks 
Relishes (Crts, Dll Pkls, Rp Olives) 
Stuffed Celery , * 

Devile^ -Egg Salad 
Lettuce Wedge/Russian Dressing 
Tossed Salad/French Dressing 
Lettuce Wedge/Salad Dressing 
Tossed Salad/1000 Isle Dressing • 
Marinated Vegetable Salad 
Ambrosia Salad 
Jellied Grapefruit Salad 
Peach-Cottage Cheese Salad 
Jellied Pear Salad 
Pineapple-Cheddar Cheese Salad 
Sliced Tomato Salad 

j 

Perfection Salad , 



Cost in Cents 
to Management 
per Serving 

* 4.58 

4.48. 

6.30 

4.7# 

8.68 
12.14 

8.50 
10.26 
10.28 

6.47 

6.82 
12.86 



8.86 

4.43 

4.90 ' 

5.34 

3.24 

7.52 

4.7,0 . 

7.87 

6.99 

6rS5 
10'.93 
- 5.95 

8.51 
1^.29 

6,98 
10.88 

8.22 
10.37 
10.95 

4.06 
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c i w ■ " ( t 0st in^nts 

Sample Menu Items ^ , . . ^ 0 Management 

Desserts 



per 'ServinE 



Apple Betty? ^0.83- 

Canned Apricots 6.22 

Banana Layer Cake • * 16.75 

Gingerbread with Lemon Sauce 10.67* 

White Layer Cake with Icing * 10.57 

White Bread . * 2.36' 

Milk, '8.00 

Cherry Pie ^ ' 13,03 

Angel -Cake/Whipped Choc Topping * . 10.50 

Gelatin Cubes \ t 2*30 

Lemon Sponge Custar3 , • 4,12 

Canned Fruit Cocktail ' 9,52 

Peach Pie ' - y.96 

Canned Pears 6.72 

Vanilla-Ice Cream' - 8.07 

Rice Custard |B- * 4,22 

Pumpkin Pie . ^ * 15.00 
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1-9 Model Exam 

1. Use Computer Program 1 in Appendix A to calculate AB for 



2. 



3. 



• 


1 3 
» j 


in L 


/ « / 






19 


16.4 


13.7 






.25-4 


22.4 


20.2 




14 


A- 


32 


29.4 


.26.7 


B« 


32 




37:6 


35, 


32-. 4 




0 




42.4 


40 


37'. 1 






46.6 


44 


41.3, 





11.1 J7.5 24.1 



2.5 
6* 



8.5 15 
12'* 18 



21.5 

2r4 



Vegetable* . Salad 



29-8 


39" 


27 


'36 


30 


3i 







Dessert 



Us% the following format.. 

Meat Potato 

Menu 1 
Menu 2 
Menu 3 
Menu 4 

to write a matrix of costs, A, 'foV the items shown in the follow 
four menus. Costs per serving may be taken from Section 1.8. - 
Menu T 



Baked Swiss Steak 

Hash Brown Potatoes 

Buttered Spinach 

Tossed Salad/French Dressing 

Canned Fruit Cocktail 



Merfu 3 

Bar-B-Q Chicken 
Whipped Potatoes 
Harvard Beets 



Menu 2 
Baked Ham Loaf 
Glazed Sweet Potatoes 
Stewed Tomatoes 
Pineapple Waldorf Salad 
Gelatin Cubes 

Menu 4 



Hot Corned Beef 
# Oven-Browned Potato 

Buttered Steamed Cabbage 
Lettuce Wedge/Russian Dressing Relishes * „ 4 
Lemon Sponge Custard^ JCherry Pie ' 

Solve this problem using the 4 menus shown above: 
A cafeteria serves 50* Menu 1, 75 Menu 2, 37 Menu 3\ and 46 
Menu fi 4 orders. * Use matrix multiplication to find the total 
cost for each type of menu." * v 
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2. APPLICATIONS OF MATRIX METHODS: 
FOOD SERVICE AND DIETARY REQUIREMENTS (U109) 

2.1 Challenge Problem 

The manager of food services for a large hospital 
has the use of a computer to assist; him in the management 
Of his service. He*understarids that a computer'is useful, 
in f inancial accounting, but 'he is hopeful that it can 
be used to help him plan ahead and, perhaps, improve his 
services. One of his problems is;/that he must keep the 
cost of food served at a reasonable price and, at the 
same tim.e , meet the nutritional requirements for balanced 
or therapeutic diets. Dietitians, of course, plan meal^, 
but the coordination of this activity With the purchasing 
and preparation of food is needed also. He feels that * 
some of this work could be reduced by using a computer 
to give information which would help in making decisions. 

The computation of the nutrient content of a menu 
i,tem is sample arithmetic once the data has been assembled 
The- sum of the products of the number of units of the 
food items needed for the recipe of the menu item and the 
nutrients, will, of course, give the total for each 
nutrient in the menu item. But, there are hundreds of menju 
items and food items; moreover, new recipes are constantly 
being added and old ones . discarded while the price of 
food items change. There must be some orderly way of 
managing aU this. ' 

It has been suggested to the manager that he should 
.s^t up files in matrix format to handle -all these data. 
He is confused about this concept, and he calls you in 
as an assistant and mathematician to the project. Your 
task is to help him formulate his problem mathematically, 



and to setup the files as 



Question: What is meant 
matrix format? 




trices. 



by setting up 




les in 



10 
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Answer: Data in the files are arranged in row§ 
^ - and columns so that matrix operations 
, may be applied to them. 

2.2 A Recipe Matr ix , ** . % • 

• * k 

• * • « % 

. A recipe matrix is one of the first files constructed. 
Ip order to make this clea,r to the manager, you construct 
a small matrix of menu items .shewing their ingredients. 
With the help of someone such .as the food' preparation 
supervisor, you create the sample reaipe matrix, as you 
did in Section 1.4 of Unit # 105. > r 



Food Ingredients 



Menu Items 

a 



a 


11 


a !2 




a 


21 










a 23 


a 










31 


a 32 


?33 











1»4 



2.3 Assignments (Optional) , % 

1. Make a limited ue'ci^pe matrix^ Consult a competent 
resource person for actual recipes, if needed. 

Questions: What does t eacfr column of the matrix - } 
represent? What does each row, of the 
matrix, represent? What does each 
element of the matrix represent? 

* Since you will be using this matrix in arithmetic 
operations, it is necessary that the numbers be in stan- 
dard units. In quantity food preparation, food is 
measured by weight rather than spoonfuls, cups, etc. 
the number of units in the recipe for each food should be 
sufficient to produce one serving for a normal diet. 

2. Create a. nutrient matrix for each food in your recipe 
file. Suggestion: Let the entries in column i represent 
the nutrients in food ingredient i. (See .Suggested Support 
Materials, inside front cover,) Many books on nutrition 
contain extended tables of nutrient values, for example, 
Nutrition, b'y Chaney and Ross, Hough£pn Mifflin, 1971, 

' . ' ■' x 11 



contains such a table. The standard nutrients to use are: 

1. Food energy (calories) 

2. Protein (grams) 

3. Fat (grams) 

4. Carbohydrates (grams) 
* 5. -Ash (grams) 

6. Calcium (milligrams) 

7. * Phosphorous (milligrams) 

8. Iron (milligrams) 

9. Sodium (milligrams) 
10. Potassium (milligrams) 

11- Vitamin A (international units) 

12. Thiamin (milligrams), 

13. . Riboflavin (milligrams) 

14. Niacin (milligrams-) 

15. Ascorbic acid (milligrams) 

Questions: What do the columns -o'f the nutrient • 
j. matrix represent? 

What do the^rows of the nutrient 6 matrix 
represent? y 

Wha£ does each" element of the nutrient 
^ matrix represent? * 

3„ Using, the recipe matrix and the nutrient matrix, 
find the, nutrient content of each recipe. 

^ Note that if, in constructing the nutrient matrix, 
you failed to make it conformable for matrix multipli-V 
cation, either pre- or post-, frith the recipe matrix, 
an adjustment will have to be made. to 

4. Prepare a price vector (a matrix consisting of one 
.row or one column) for the foocft in your recipe matrix." 
Make your prices conform to "reality" as much as possible, 
keeping in mind % that tie hospital nlay^hlTYe a reduced price 
for some items which are purchasedAri quantity.' 

< ^ Find the cost of each* recipe v using the recipe* matrix 
and the price vector. Which recipes are the most expensive? 

12 
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Referring % to Assignment 3, compare the nutrient content 
of the most expensive recipes with' that -of the least 
expensive. Jk 

5. ^Prepare a serving matrix, that is * a matrix con- 
sisting, of the number of servings of each menu item per*, 
day fox -same period of*time. One week would be a/satis- 
factory period. Keep in "mind that no iT every recHpe need 
be. prepared every day. Using this* information, and the 
previous matrices which you constructed, how would you 
find a matrix showing the amount of eacli food ne'eded per 
day,? How would you find the cost per-day of- the food 
required? Perform these calculations. 



2.4^ Diets 'Meet ing Certain Requirements 

At this point, we t;urh our attention to a different 
type of problem. Suppose that a doctor has ordered that 
a patient f s diet meet a minimum daily requirement ^of 1 
unit of thiamin, 2 'units of niacin, and 3 units of iron. 
If we select three menu items (food items could also be 
used) whrch< contain these vitamins in the following 
quantities , * 




6 




M 2 


M 3 


Thiamin 


" 1 


0 


1 " 


Niacin 


0 


2 


- 3 


Iron 


4 


0 


1- 



.what portion 6T"^5*s£ recipe (pr food) should be served 
to the patient to meet the minimum requirements? 



Solution : 

Let x, 



the pdrtion of the first item needed, 
the portion of the second item needed, 
the portion^ of the third item needed. 



Then for thiamin x, + x > 1, 
for niacin 2x 2 + 3x >_ 2, 



and iron 



4a? j > 3, 
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These conditions in themselves are not enough to 
guarantee that a triple o.f values (x 1 ,'x 2 , x 3 ) will 
give the proportions that the patient needs. This is 
because the system of inequalities may have solutions in 
which one or more of the variables takes on a negative 
value. For example, the.%^iple of values x x = 2, = 3, 
x 3 = sa.tis f ies each of the inequalities. However, we 
can eliminate the possibility of such solutions by 
requiring the variables to satisfy the additional three 
inequalities x 1 > 0, x 2 > 0, x 3 > 0. These inequalities," 
merely formulate the obvious statement that menu portions 
are never negative. 

Now, we can .certainly satisfy all six of our 
inequalities with 'huge values of.Xj, x 2 , and x 3 . But to 
"meet tjie 'patient f s dietary requirements more economically, 
we look for non-negative solutions of the equalities 

This suggests the following system of equations: 



lx n + 3x 

2 3 

4x 



= 1 

2 

*s e 3. 



1 -0 1 






0 2 3 




x i 


4 0 1 








X 



The matrix representation of this system is 

2 
3 

= B * 

Verify that this is true by multiplying A'and X,' and 
showing that the product equals B . The results should 
be the three Original equations. 

* 



To complete the solution, let A" 1 bfe the 
matrix A (if there is an inverse "for A).* 
# 



inverse 



It- Is to be remembered that not alt matrices have Inverses. The 
Inverse exists if and only If |a| (determinant of A) has non zero 
value. Also a matrix B Is called an inverse of*A If BA » I - AB. 
Itideed, if the inverse* exists therw It Is unique. For more details 
see any s-tandard book jon linear algebra. 
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The 
and 

If * 



then 



A "*AX = A - X B 
IX = A - J B 
X * A - a B. 



1 



A-i = 



4 

1 



2.5 Research Project 

Find the recommended daily nutirtional requirements 
for a normal person. You can use yourself as the person 
in questi07r-in oxder to specify characteristics such as 
sex and weight. ' Construct a, recipe matrix of the menu 
items that you eat for a given week. Does your % diet meet 
the recommended daily nutritional requirements most days? 
Does the ave'rage^of the daily^take for a week meet 
these requirements? 
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1.6 Model Exam * ' 
General Directions: , ' 

In tjiis test, do not perform any computations or assign 
numerlcflljralues to the elements of the matrices^you-construct. 
To jtidlcate the structure of a matrix, simply 'label the contents' 
of- the cojumns and rfws as shown below. 

s Name of Items * 
* Represented by 

the Columns 



13 



23 N 



Name of I terns 
Represented by 
the Rows 



1, 



3. 



'4. 



In food preparation, a subassembly is a £>art of a menu item 
which may be prepared separately, for example, salad dressing, 
frosting for a cake, or pie crust' doughy Se* up a matrix to 
show the food content of a set of subassemblies. 

Suppose you wished to compute the cost of each subassembly. 
What sor^t of matrix would you need, and what kind of informa- 
t4oA~sfi?uld it contain? Show the matrix operation which 
would compute this cost.. Represent the matrices as indicated 
in the general directions. * 

Show how you would compute the nutrient content of the sub- 
assemblies. r 

Below are two nonsense matrices called A and B. 



Gooies 



« u 



Fudge I terns 



Unit cost of fudge items: 



B 



Interpret the product BA. 



>..] 



20 . 



16 



ERIC , . ; -V 



In this section we have been concerned with the following 
activities: 

a. Deciding on the kinds of information needed to answer 
, certain ^qdestions^ , 
br Collecting data 

c. Displaying data in a matrix * * 

d. Deciding on, the matrix operations to be applied to ' 
obtain information from data 

e. Making computations 

f. Interpreting results — ** ^ 

* V 
Which of these activities do you consider the most difficult 

and which the v easiest? Give reasons for your answers. ♦ 

(This question is evaluated on'the reasons you give for 

your answer.) * ■ 




The' Project would like to thank Kenneth R. Rebman of ^ 
California State University, Hayward, and George Springer of ^ 
Indiana University, Bloomington, for their reviews, and all 
others who assisted in the production of this unit. 

„-Riis unit was field-tested and/or stude/it reviewed by 
Ellen Cunninghan, SP, St. Mary r of-the-Woods College, "St. Mary-of 
the-Woods, Indiana; Robert, M. Thrall, Rice University, Houston, 
Texas; Philip D.,Straffin, Jr., Beloit College, Belott> Wisconsin; 
Michael A. GrajeR^ Hiram College, Hiram,, Ohio; 'Donald G. Beane, ' 
The College of Wooster, Wooster, Ohio; Dina Ng,,,CaI I fornia Polytech 
,State University, San Luis Obsipo, California, and; Henry J. Osner, 
Modesto Junior College, Modesto, California, aYid has been" revised 
on the basis of data received from these sites. 
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3. . ANSWERS TO' ^XERCISES CU105) 

Answer to question on page 4. N • 

,Postmultiply the matrix by a column vector consisting of I's. 
tn the above case r use: 1 

1 
1 

















/ 




1. 


1 .5 


1 


0 




5 




12 


eggs 




1 0 


.25 


0 


X 


10 




5.5 


flour 




1 0 


0 


t) 




2 


a 


5 • 


, sugar 


t 


1 .25 


.25 


0 




50 , 




8 


butter 




0- 0 


0 


1 






- 50 . 


beef stroganoff 



/ 



Note that it was necessary to post^mul t ipfy the ingredient 
matrix by the number of recipes vector. This vector was formulated 
as a column vector. Verify that the proper sums of products were 
formed in order to give the total quantities of eggs, flour, et$. ** 



2. [70 10 25 50 100] 



3- 



1 .5 
1 6 
1 0 



12 
5.5 
5 



50 



1 0 
.25 0 
0 * 0 



1 .25 .25 0 
0 0 0 1 



5 2 1 

10 3 1 

2 3 10 

50 3 -6 



6420. 



"12 6.5 11.5" 

5.5 2.7S 3*5 
5 2 1 

8 * 3.5 3.75 
50 6 0 



4. [70 10 :ZS 50 100] x 



12 6.5 

5.5 2.75 

5 2 

8 3.5 

50 6 



11.5" 
3.5 
1 

3.75 
0. 



[6420 1307.5 1052.5] 
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ANSWERS TO MODEII EXAM (U1Q5) 



1 . t Input Matrices 



Y 



. 13.000 i 
19.Q00 
25,400 
32.000 
37.600 
1(2. 400 
46 .$00 

34.600 
32.000 
0.000 

Ha^ P 

78{.599 
. 1182.199 
1595.639 
2048.000 
• 2420.959 
2747.040 
3020.36a 



10.400 7'.70*0 

,16.400 13.700 

22.400 20.200 

29.400 26.700 

35.000 32.400 

40.000 37.200 

44JO0O « 41.300 

> 5.1u0 «JVlOO 

2.5-00 8.500 

6.Q00 ' 12,. 000 



roduct 

138.500 
220*059 
306.739 
396.859 
473.6)9 
539.439 
595.459 



2. 



\ 



HI 
H2 
H3 
M4 



^ Meat 
fi4.52 

32:93 
46.24 
34.50 



325.099 
514.^99 
714.739 
925.499 
1103.659 
1257.039 
1386.859 

Potatoes 
8.49 
8.67 

in 



17.500 
15.000 
18.000 



522.099 
825.099 
1144.099 
1481.599 
1766,199 
2011.599 
2213.899 



'24.100 29.800 
21.500 - 27.000* 
24.000 30.000 



39. 
36. 
36. 



poo 

000 
000 



721.699 
1139.300 
45Z8.539 
2044.099 
2.436.259 
2774.639 
3060.259 



899.199 
1420.000 
1967.719 
2548.399 
3037.479 
3459.519 
3816.679 



11 58. 600 
1824.60a 
2524.199 
3267.600 
3892.799 
4432.799 
4888.200 



Vegetables 
. 6.30 

6.47 

2.61 

5.16 



Salad 
6.85 

. 8.86 
6.99 
<?.52 



Desserj^ 
9.52 
.2.30 
4.12 

' 13.03 



3." To find the cpst of one s,erving\o?* an individual menu, we need 
to sum the costs of the component menu items:" i.e., we need to 
find'jthe sum of a row of the matrix oFcostSj A. This can be' 
done i>y post-multiplying A by the vector 



> * The result will be a vector whose components are the' costs of 
the menus: 

[cost (menu 1) cdst (menu 2) c<Jst (menu 3) Cost (menu 4)] 

Ip^intf the total cost to the cafeteria, multiply this vector 
by the number-of ^servings vector ) 
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50 
75 
37 
46 



J 



Notice pre-multipjy ing matrix A by C givers a cost vector 
whose elements wo^u ^represent? the total cost for all the 



menus of: 



[meat potatoes vegetables -salad dessert] 

but not the cost of the individual menus. We see that it is 
important to consider carefully vjhat the entries in a product 
matrix represent and the order of mill tl plication. 




ANSWERS TO MODEL EXAM (U109) 



1-2. Let/,, f t , f )t 



and S, 



Let x?, , c? 2 , c? 3 , 



reprfssenr^uantities of food items, 
be- sub-assemblies "containing' these items* 
be unit costs of these food items. Then 











s 2 


s 3 


* 


to, a 2 o,. . 


.] x 




4 




» 


(Cost vector) 




4 


4 


4 








4 


4 


4 
















I 








i tern 




, Some f . may be 


ze ro. 







■.Cost of each 
/ sub-assembly r 



3. 



» 5. 



) 



Be careful of tftis, problem. It is deliberately vagiie. One 
interpretation would be to'compute the total nutrient content 
of several sub-assemblies. Then the problem may be set upj 
as follows: y 



[a, a 2 a 3 . . .] 

(quant i ty of each 
sub-assembly 
prepared) 















","> 


*3 


fx 


x n 


X 13 


4 






*23 


4 


X 31 




X 33 



Total of each 
nutrient in 
set of sub- 
assembl les 



^(Nutrient vectors showing 
s " unit quantities of each 

.for food items.) 

The product B * A ■ a cost vector for Gooies. 

Individual answers are acceptable. However, for the type of 
problem we are solving here activity £ and activity £ are 
probably the critical, ones. Activity f»\s important, but 
should not be difficult 4f the others are performed properly. 
Computations should not be considered difficult if a computer 
is, used. Deciding on the matrix operations to be applied 
should follow from the matrix structure. 
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APPENDIX A 



C*****APPLICATIONS OF MATRIX METHODS - PROGRAM 1 

C*****THIS PROGRAM PERFORMS MATRIX MULTIPLICATION OM ANY NUMBER OF PAIRS 
<OF MATRICES 

DIMENSION A(20,20),B(20,20), C (20,20) , IHEAD(40) 
C*****NR - LOGICAL NUMBER FOR THE CARD READER. 

NR » 2 <r 

C*****NP « THE LOGICAL NUMBER FOR THE PRINTER. • ^ * 

. -HP ** 5 

C*****READ "fDJNT I F I CAT 1 6N (*RD. LAST CARD SHOULD CONTAIN /* IN COL- W 
100 READ (NR, 1 , END*70) IHEAD 
1 FORMAT 2) 

WRITE(NP,3) IHEAD * . 

3 F0RMAT(IH1 ,40A2//) ' % 

C READ THE -NUMBER OF ROWS AND COLUMNS IN THE FIRST MATRIX. 

READ(NR,5) NROV/, MCOl » ' 

5 * FORMAT(215). - 

C*****READ THE FIRST MATRIX. . P 

DO 10 I - 1 ,NROW % 

READ(NR,7) fA( I , J) \ J»1 ,NCOL) 
7 FORM AT (\0F 5.0) ' k , 

10 CONTINUE / 

C*****READ THE V^JMBER OF ROWS AND COLUMNS IN THE SECOND MATRIX. 
READ(NR, V 5) ffifcw,MCOL 

DO 20 I - 1 ,MROW , > 

C*****READ THE SECOND MATRIX. ' ' 

20 READ(NR,7) (B( I , J) , J - 1 ,MCOL) «" M 

C**A**CALL THE SUBROUTINE TO MULTIPLY TflE MATRICES. k # 

.CALL MATMY(ArtNROW,NCOL ,B,MROW,MCOL ,C) 
C*****WRITE THE INPUT MATRICES. % ^^S*** * 

WRITE(NP,25) • 1 

25 FORMAT (5X, * INPUT^TRICES'//) 

DO 35 I - 1 ,MROW 
"WRIVE (NP,30) (A( I , J) , J - 1 ,NCOL) , 
30 F0RMAT(5X,10F10.3) 
35 ' CONTINUE 

WRITE(NP,40)' - 
40 FORMAT (///) 

DO 50 I - 1,MR0W 

50 WXlTE(NP,30) (B(l ,J),J » 1.MC0L) * *** 

C***f^WR|T^ THE MATRIX PRODUCT. 

WRITE(NP,55) 
55 FORMAT(///55(, J MATRlX PRODUCT 1 //) 

D0"60 l»« 1.NR0W 
60 .V^ITECNP^P,) (C(l ,JhJ-t,MCOL) 
C*****RETURN FOR A NEW PROBLEM. 

GO TO 100 , "" ^ 

70 CALL. EXIT 
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SUBROUTINE MATMY(A,NROW K NCOL,B,MROW,MCOL,C) 
DIMENSION A(20,20) ,B(20,20) ,C(20,20) 
NP.« 5 

IF(NC0L-MR0W)20, 10,20 
20 & WRITE(NP,15) • > 

15 FORMAT(SX,' NUMBER OF COLUMNS IN FIRST MATRIX MUST EQUAL NUMBER OF 
> 1R0WS IN SECOND MATRIX 1 ) 
CALL EXIT 
10 DO 40 I » 1,NR0W 
* *D0<40 J - 1,MC0L • 

C(l/,J) »0.0 « 
DO AO K » 1,NC0L 
40 X(*,J) - C( I , J) .+ A(l,K)*B(K,J) 
RETURN 
; END . 
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1. MARKOV CHAINS (UNIT 1071 

1.1 Introduction * 

Jj^ order to understand what is meant by a* Markov 
chain, consider the following situation. " ' 

In Ascertain class, a teacher Iras observed that 
students' performance on tests is affected by how well 
or poorly the^have done on the last test taken. In 
particular, 801 of the students who did well on the last 
test will rate we*ll on the next one, 15% will be average, 
and only 51 will be poor. For those who were rated as 
average, on a test, 60% will continue to be average on 
the next test, while 10% will do welj and 30% poorly.. 
For those students who were rated as poor,^nly 1% will 
do well, 15% average, and the remaining 84% wili continue 
to<rate low in the next test. We can" think of this as a 
process which will continue through several tests. For 
the sake of discussion, we will ignore. any factor which 
might upset these predictions. . 

1.2 Tree Diagrams 

These probabilities 'can be represented by a tree 
diagram. Let the ratings be labeled as: 

a^ - good 
a 2 




a % = 

3 i . 

where a^ a^ a^ represent the current *test score *for ' 

any student. If we^ start with a student who has received 

a good grade, we can show the possibilities 'for the 

next test with a diagram, like the one below. 
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Figure 1 . 



In Figure' 1, the lines drawn from a 2 , or branches 
from a 2> are labeled with the probabilities that any 
one of them will lead to the next event; that is, getting 
a gratfe of good, average, or poor". Since these are the 
only possibilities, one of them must happen if the student 
takes another test. For this reason, the sum of the 
probabilities stemming from any one point must equal 1. 
Otherwise, some event could happen which is ngt accounted 
for. 

Since we know the probabilities for the students 
who receive average or poor grades, we can extend the 
tree in Figure 1 to show this information., 

Figure 2 shows the"probabilities through a series 
of tests. The branches stemming from the left-most 
a 2 point to the three outcome points for test 2. The 
branches from each of these three points indicate the 
probabilities for test 3. 



Figure 2. J 

Representation of probabilities in the form 
above is called a tree diagram. 

N 

1.3 Calculating Probabilities" From a Tree Diagram 

Suppose we wished to khow the probability of getting 
a good rating on the third test*£ox a student who had 
♦received a good rating on the first test. If we examine 
Figure 2, we see that^ appears three times in the 
right-most column, which indicates the. outcomes for the 
third test* These three paths are along the branches: 

. i 1 2' *» 

>. T * 

According to the rules of compound probability, , 
the probability of one 'event following another is the 

product of their -probabilities . Theref ore,,, the 
■■■ ■ 1 ■ i — ■ 



For the product rule to hold, the events Jn question must be 
Independent. A . * 

9 r ' * 34 



probability for each of t*hese series of events is 

<* 2 d 1 -+ a 2 * ,8(.8) = .6^ ' ^ 

*a 2 a 3 a 2 = .05(.01) = . 0005 

Since one of these events must occur if the test is 
taken, and siase the events are mutually exclusive (cannot 
occur together), the probabi 1 ity ,of receiving a rating 
of good on the third test is'the sum of the probabilit ies^ 
of completing the* paths shown above. 



Thus, the possibility of reaching a^ in two moves 
is .64 +'.015 + .0005 *^j6555. If we were to explore 
this process beginning at any state, we could compute 
the probability for any subsequent state in a similar 
manner . - k 

1.4 The Matrix Representation of a Markov Chain . 

The format .of the computations made in Section 3 
suggests that the information in the tree diagram could 
be structured as a matrix, 



M = 



a 3 



0 


a l 


a 2 


« 




(good) 


(•average) 


(poor) 


(good)^ 


.8 


.15 


,05 " 


(average) # 


.1 


.6 


.3 


(poor) ■ 


.01 


.15 


*.84 



to which we attach the meaning: the probability of 
going, in d one step from \ 



a l t0 a l = Pu = - 8 
a l t0 a 3 = PJ3 = - 05 - 
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Exercise 1 : 

1 k 1 

Complete by giving the ^meaning <pf and probability for 

- • — : " p 2r . 

^ ■ p 22 - , 



= P 33 " — 



^ ' * 



The process described in the preceding sections* i;s * 
-Bn example of a Markov chain. Such a chain consists of 
a series of states, an'd the .probabilities oT passing, to 
a new one in some defined process. As the example shows", 
each state is* always dependent on the one that precedes' 
it. 



1»5 Experiment 1 ^ _ 



• Use the test exanrple and ^he computer to compute the probability 

of being In state or a^ or|a 3 beginning, in any one of the* three 

states after the second test; the thi.rd test; the fourth test* 

JtiiEUL 1 Do yo u see that multiplying the matrix J\, 
by Itself according to £he rules of Wrix 
multiplication will give the desired proba- 
bilities for the'second test? d ^ 

Interpret M-M-M » M 3 -. 

What is the probability that the' third test.will be rated 
average If a student was rated as poor in the first test? 

Extend tjj* tree diagram in Figure 2* to show 'the continuation 
of the process through four steps. Calculate \he probabil i t ies of * 
being Instate a ^ or a^ or aftcff the fourth jtest, using the tree 
diagram and the method shown on pages 3 - Compare this with the 
results^pbtalned by calculating H k ^ Does the matrix method proaVce' . 
the same results? Do you see that use of matrices in this problem, 
• 5 
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makes thejralculations simpler and more iikely to be accurate when 
using a computer than when working manually from a tree diagram? 

1.6 Experiment 2 - 

Assume* <;hat women's occupations could be classified as follows: 
Housewife-, full time =w 

*2 



Housewife,/part time work outside ■ 



Full time-work outside home t = 
Full time professional career = * • 

A sample is taken of wpmen who have at least ^one daughter* ^ 
The 'following trends were noted: Of those daughters whose mothers 
haji been fulj-time housewives, 50* were classified as^, 25% as 
W 2 , r 20fc as W^, and 5% as W^. For those whose mothers were housewp/es 
and w\>rke<( part time outside the home, 60S of the daughters also * ^ 
did so^but^ 1 5% became full-time housewives and 15% worked full time, 

'with 10* having a full professional career. The daughters <rf*the 
full-time workers were distributed as fallows: 20% full-U-me 

.housewives, 25% part-tim^ workers* k0% full-time workers, and 
15* professional women. Finally, the daughters of professional 
women- were distributed in this fashion: 30% housewives, 20% part- 
time workers, ^20% fulj-time.workers, and* 30% professional wdmer>. 

Construct a matrix to represent a Markov chain for these data. 

• Assuming that this trefld continues, find the probabilities 
that a woman will have the" same career as her grandmother. 

Calculate M 2 , M 3 , M\ M 5 , M 6 . ... up to any power you wish 
for this matrix. You are now able to mak* long term predictions 
about this process. What seems to be haffoening? 
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1.7 Model Exam fUnit 107) 

3- In the matrix representation of a Markov chain, 'what do the 
elements of the matrix represent? 

The row sum of any Markov^ matrix must be 1, Why? 

The. following diagram represents a maze. Each compartment 
can be considered a "state" of the sy>*em. if a rat is 
placed in compartment a^ 9 what is the porbabiiity that he ~* 
will escape from the maze after a given number of trail's? 
(A trial consists, of a .move from one compartment to another.) 
Where there is only one way out of a compartment, the , 
probability of choosing that exit is, of course, 1. Clearly, 
if a compartment cannot be reached directly from another 
compartment, the probability of passing between these two 
is zero. Movement when there are multiple exits is considered 
equally probable. 

J 



a (Exit from 
tlje system) 



zero 




Rat Maze 

Draw a tree diagram to represent this system, and then set 
up a matrix of possibilities. Compute the probability a , 
rat leaves the maze after three trials. 



( 
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APPLICATIONS OF M ATRIX METHODS: 

— — 



FlXEtf- POINT AND ABSORBING MARKO V CHAINS (UNIT 111) 



2 . 1 Challenge Problem ^ 

Competition is a way of life for the producers of 
many things, from TV shows to detergents One problem 
stems from the fickle nature of consumers. They tend 
to Switch from one product to another. 

For ex'ample, consider three TV networks which are 
competing for viewers in a given time slot. Three 
shows, SUNNY DAYS, LOTSA GUNS, and MOON dlfcuRES are* 
all broadcast on Tuesdays/ at 7:00 p.m. Su^eys taken, 
indicate that for those who watch SUNN.Y DAYS one week 
theVe is a probability that 601 will continue to watch 
it'the next week, while 30% will probably switch to MOON 
CREATURES, and 10 ^ to LOTSA GUNS. For persons who watch 
MOON CREATURES there is a 501 chance that they will 
continue to do so the next week, with 40% changing to 
LOTSA GUNS, and 10% going to SUNNY DAYS. Finally, those 
who watch LOTSA GUNS have a probability of 70% of staying 
with the show the following week, and a 30% probability 
of switching to MOON CREATURES. 

Let us formulate this information as a transition 
. matrix : , 







I 






.' SD 


MC 


LG 


SD 


.6 


. 3 


.1 




,1 


. 5 


.4 


LG 


.0 


.3 


" .7 



SUNNY DAYS started out wifch 70% of the audience, 
MOON CREATURES had 10% .and, LOTSA GUNS* had* 20% . In 
spite of the good start, the cast of SUNNY DAYS were 
worrying about their jobs at the end of the fifth week, 
and were, definitely out of a job by the tenth ,week. 
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Could this have been predicted? The answer is yes, if 
it is assumed that the trend shown in -the survey continued. 
In fatt-jVit^is possible to predict tKat Eventually SUNNY * 
'DAYS will have 9.4* of the viewers, MOON CREATURES will 
have about 37.5%, and LOTSA GUNS will horti 53.1% of the ' 
audience. JVhen that point is reached, there will be no 
further 4 changes. The probabilities become fixed. 

In Unit 107 we showed tha,t we could predict the 

iSt'ate of a Markov chain as the process went through 

several stages. We did this my multiplying the matrix 

by itself, or raising it to a power. We did not explore 

the possibility of the matrix reaching £ steady> state , 

- that is, that raising, the matrix to higher and higher 

powers no % longer changed the probabilities. We consider 

this situation now. . * 

* 

2.2 Regular Transition Matrices 

a A transition matrix' representing a Markov chain is 
said to be regular if some fower of the matrix has only 
positive .components. 

The transition matrix from the challenge problem is 

an example of a r eg u 1 ar-onat r ix though the .original 
, " * , \ 

qQ^matrix has a zero element, if^we taloe xhe second power, 

we find that all of the elements are positive . Verify 

this by constructing the matrix and u^ing Program 7 in 

Appendix A to find some power o£ it. 

^ 2.3 Fixed-Probability Vectors • $ ^* 

A row vector that consists of non-negative elements 
whose sum is 1 is called a -probability vector. From this 
definition each swingle row of *a transition matrix is** a 
probability vector. If a transition^ matrix is regular, 
then after a number of steps, sometimes a large number, 
the probability vectors (rows) tend, to become the some 
and remain fixed. To illustrate thi^s , we use the" 
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challenge problem. When the transition matrix for TV 
shows is raised to the 20th power, it becomes 

p. 0937 0.3749 0.5312] 
0.0937 0.3749 0.5312 
_0.0937 0.3749 0. 5312_ 

When this happens the process is in a "steady" 
•£tate. The probabilities will not change in future 
steps. The row vector which gives these probabilities, 

[0.0937 0.3749 0.5312] 

is called the f ixed-proba-bility vector. *i 

2.4 Calculating a Fixed^robability Vector 

We can, of course, search for a fixed-probability 
vector by raising a regular transition matrix to a power,* 
continuing until the fixed state i's reached. With a 
computer this is not particularly difficult, although it 
may converge slowly and the result be only approximate. 
There is , however, a direct way of obtaining the fixed- 
probability vector. 

Ifsp is -a fixed-probability vector for 'a matrix A, 
then it can 'be shown that pA = p . If we 'use this theorem, 
we can set up a system of equations which can be solved 
for the vector p 

2.5 Experiment 1 

Verify that multiplying the transitioo matrix for the challenge 
problem in Section 2.1 by the vector ohtained in 2.3 gives the same 
vector as the product. Thus* 





0.6 


0.3 


o.f 


[0.0937 u„.37**9 0.5312] 


0.1 


V'5 






0.0 


0.3 


0.7 



- 0.937 0.37*9 0.5312! 



Use Program 1 . v 
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2.fr A F ixed Probability Vector from a System of Linear 
Equations 

* * • 

Since, Tor the fixed-probability vector p and the * 
regular matrix A, pA ■ p, then, if 



and 



A = 


_.6 


r 

.4 




P = l> 2 

r 




1*1 


[.: 


f 
.4 


- [* 2 *,] 



If we carry out the indicated, multiplication , we 
obtain 

** 



or 



**1 + Ax 2 " *2 



and because the sum pf\,any row probability vector must 
equal 1 ; 



*1 + X 2 ' l - 



2~. 7 Experiment 2 * 

• Use Program 6 in Appendix A to solve the three equations in 
two unknowns which were develop Section 2.6 We r«state them as 



~ x l + ' ° 

*1 ' " 0 

« 3 •+ « 2 - 1 



Interpret the result. 

* 

Show that the same result could* be found by raising matrix A 
to a sufficiently high power. Use Program ? . u 
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2.8 Experiment 3 



Four companies are competing against each other with products 
in toothpaste. A Survey shows that the shift 'from one brand to 
the A other qan be presented by this transition matrix. 



Brand 



C 
. 1 



.7 .1 
■ T .7 
.2 ♦l 



What is the long , term prediction for each company'^ share 
of the market? * fc 

%i 

# What change would occur, if any, if company D changed its 
product and a new survey showed the transition matrix to be: 



B 

.2 
.6 
. 1 



C 
. 1 
.2 
.7 



2.9 Absorbing Markov Chains 

Some Markov chains contain states from which, once 
entered, departure is no longer possible. This state^ 
from which there is no return is called an absorbing 
staie. We might have considered the rat maze problem 
in "Unit 107 as having such a' state if we, assumed that 
once the rat left the maze it could not go back in. 

We can recognize 'an absorbing state from a transition 
matrix. Any state, a^ 9 for which the element a., is 
equal to 1 and all other elements of that row are zero * 
is an absorbing state. As an example, recall the rat 
maze problem on page 7. The transition matrix is" 
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It should be immediately evident that a, 



is an 



absorbing state. The probability of going from a 5 
that a, 



2 > uy or a 4 is zero in each case. We note 
ss equals 1 and all the other elements in that, 
row are z^o. * 

For a Markov chain to be an absorbing chain it 
must be possible to get from^any non- absorbing state 
to an absorbing state. One way to recognize this from 
the matrix representation of a Markov chain* is to 
examine the columns which contain the 1 for the absorbing 
spates. ^For e v ach such column the remaining elements 
mu^t not be all zeros if there is to be a transition " 
to this absorbing state. For example, the following 
Markov chain contains one absorbing ' state but is not 
an absorbing chain. 



.3 
0 
0 
'0 



0 .2 .3 

10 0 

0 c .8 .2 

0 0 0 
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The state a $ is an absorbing state, but there* is no way 
f t0 reach u from any other state:* Since it is the only 
absorhing state in this chain the- chain in not an absorbing 
chain. \We could verify this by drawing the tree diagram 
for this chain. ♦ ♦ 

2 ' 10 Exercise for Absorbing Markov Chains ' \ * 

State whether the following transition matrices are for 
absorbing or foi 1 nonabsorbing Markov chains.^ Why? . 
a. 



o r 

.5 .5 



i 



d. 



F T 0 



, 2.11 A'Second Challenge Problem 

. , The Ace Collection Agency decides to add a service 
for its department store customers, and, perhaps, 'improve 
its own business. The °pres ident of Ace has observed 
that some department stores turn over their bad accounts 
for collection at varying times, while .other companies 
rarely use the agency. The latter companies simply 
write-off unpaid bills after repeated attempts at 
collecting on their own. The president of Ace proposes 
that, for a reasonable fee, his agency\ili analyze . 
the paying habits of customers who have charge accounts 
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♦with department stores. This analysis will produce, it * 
is claimed, information that will enable *a store to 
decide on a policy for turning over bad accounts to a 
collection agency. At the .same time, the analysis .will 
give the 'store a\ray of calculating how 'long, on the 
average,' it takes for accounts to be either paid up or 
classified as bad. 

The manager* of Homer "Department Store, after seeing 
/ this analysis service advertised, decides to try it, 
but he insists that the' method applied to determining > 
any policies for the store be made clear to him. before 
they are effective. He asks that a representative from 
Ace give an explanation of how it will be determined that 
a debt will probably end up as bad, or how longNdebts 
are likely to stay in various stages of being overdue. 



The ^^representative agrees to give an explanation. 
He begins with^^ypothetical case. Suppose, he says, 
that after studying\vour accounts it was found from past 
history that Jrour customers' paying habits could have N 
probabilities attached to them. These probabilities of - 
changing status from mo/ith to month are shown in Jable I. 



V 



Prob 



>^fYt-tf 



TABLE 



ies of Future Debts of a Typical Customer 



Future*Months in Arrears 



Present 
, Months 
tn' 

Arrears 



\ 



0 
1 
2 
3 

■ h 
5 

Paid-up 
Bad 



,.60 


.15 


0.0 


0.0 


0.0 ■ 


• 0.0 


■25 


0.0 ' 


.20 


.35 


.25 


0.0 


0.0 


0.0 


.20 


0.0 


.10 


.20 


.10 


.27 


0.0 


0.0 


.13 


0.0 


• 05 


.10 


..20 


.18 


.37 


0.0 


.10 


0.0 


.02 


.03 


.07 


.30 


.28 


.15 


-.15 


0.0 


.01 


.Oh 


.04 


0.0 


.25 


.45' 


.06 


.15 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0* 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 
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.In this table,, the status of accounts is given in 
months overdife. If there are only current charges, this 
appears in the 0 column. The entries in the table are 
the probabilities of changing status from one month to , 
the next. For example, a customer who is two months 
in arrears (the row labeled 2) has ^ probability of 
.It) of having paid sufficient amounts o^n his account 
to be classified as having o^ly current charges next 
month. The same customer has a probability of .30 of. 
being stUl 2 months behind in the next month. 

Some members of the Homer Department Store had 
sufficient mathematical training to recognize that this 
fcaMe (Table I) could bejcons idered an absorbing Markov 
|chak. However, they wSre not advanced enough to know 
how the questions that- were. a*ked by the president could 
be answered from this information. More explanation 
was needed. 

2.12 S^arnterd Form for an^Absorbing Markov Chain 

Before answering some. of the questions raised, it 
is necessary to rearrange the absorbing Markov chain to 
standard form. TJiis, requires^ jnterchang ing some rows' 
and columns of the matrix so that absorbing states are 1 
placed first. We illustrate tjhis. 'using fc a simpler 
matrix than the one deriued ftom Table* I. 

Given, the absorbing Markov c^ia^n in mafrix form: 



3* 





a l 






a 4 




a l 




0 0 


0 


(T 


i 
















.5 


• - 0 


.5 


0 




a l 


0 


4 .a 


0 


.2 

•> 




% 




0 


0 # 


1_ 





we look^for the absorbing states . From the % disfcussioV 
in Section 2.9 we should recognize states <x 2 and a* as • 
absorbing states. To obtain i standard form for the * 



matrix, interchange the columns and 'rows so that there 
is an identity matrix in, the upper left hand corner. 
As you can see, this can be accomplished by placing 
* the Absorbing* states first. This does 1 not change any 
relationship, i.e., the probabilities of going from 
one statp to another are preserved. 



*2 . 



a 4+ 

"1 

'0 
.5 
0 



*4 
0 

1 

0 

.2. 



0 
.5 
0 



2.13 Partitioning the Standard Form 

Once the matrix is in standard form, we can proceed 

to partition it in such a way that four matrices are 

formed Jfln the original. Later we will see that we 

can use^hese new^matrices to help answer our questions 

about the charge accounts? We partition the matrix in 

4 this example so that there is an identity matrix in 

upper left hand corner. Thus * 

• * 

1 0 | 0 . 0 



0 1 I 0 0^ 
,5 * 0 "t « 0 - .5 



I 

.2 J .8 



-.0 



'or 



and we label each new matrix as follows, 



C'3-rC 



ffl 
0 



"5. 
0 



.3 -[ 



.51 
o 



m If we review the original matrix, we can see that the 
entries in £ are the probabilities of being absorbed, and 
the entries in f are the probabilities of being in non- 
absorbing states. This is trute because .5 is the 



probability of going from a 2 to , and .2 is ^the proba- 
bility of going from a 3 to a 4 > (q 2 and are the two 
absorbing states. Similar statements can be made for 
the entries in # .» s 

There is a theorem (which we state but will not 
prove herej) thaft ,is useful for our. purposes. On the 
average, the number of times a process will be in each 
nonabsorbing state can be* found by calculating N = {1 - T) 
where I is an identity matrix, and T is the matrix formed 
by the partition just made. 

» 

Program 8 in Appendix A can be used to calculate 
N - (I - T) for our sample problem., However, since 
-this is a very simple matrix we will do the calculations 
by hand in order..to illustrate the intermediate steps. 



I -. T 



0 
.8 



-.8 



r 



(I - Tj 



Therefore N 



i~c..s 


- 1 


8 1 




a 2 




10 


5 


X 




8 


io; 


F 


T 



io s; 

8' ' 10 



2.14 .Making Decisions Based on Probability 

All the information in a Markov chaiat consists* of 
probabilities, but in the case of the xharge accounts 
these probabilities werd^ based on the past'history of 
a Jarge .number of people 1 s paying habits. They are 
likely to be fairly predictive of the future. In the 
absence of any other knowledge, past. history forms\he 

-49 



18 



best basis for making decisions about future events. We 
Kill interpret the entries in' matrix N and see how they 
# form a basis for additional information upon which some 
decisions might be* made. 

^JJie interpretation of the entries of matrix N is 
thi^. Starting in one of th'e npnabsorbing states, say, 

the mean number of times in state before absdrp- 
tios is -g- and in state a $ is |-. The total time before 
absorption for state a 2 is i£ + | or i£. A similar 
interpretation could* me made for. the other entries v 

Question: If this matrix had come from the 
charge account problem, how would 
you interpret the" entries? 

Answer: If a customer was* in state a pt we 
would predict that on the average 
in Z3$ months he would either be 
paid-up or beco"me a bad account. 

The row sums of N.give the average time for each - 
state before it is absorbed. If we wish to find the 
row sums for some larger matrix, using the computer, 
we can use Program 1 in, Appendix A and multiply N by v 
a column matrix consisting of JL's. For example, 



10 
IT 



5 



8 10 
6" V 



15 
IT 

9 



This is a Convenient way to calculate row sums for a 
large, matrix* The column mat-rix should have enough 
l's to be conformable for multiplication. 

2.15 The Probability of Reaching a Given Absorbing State 

We still have the question concerning the probability 
of a give/i absorbing state'a? *hje final one. If "you will 
accept another theorem, we can answer this quest ion . 
According to the theorem, the product of N, ^ust computed, 

. % 19 
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and the matrix S from the partition on 'page 17 is a 
matrix which gives the probabilities of ending up i^ 
.a given absorbing stat/e. From* our example s 



/ 



Then 



10 " 


-5 




5 




TT 




' S = 




0 


8 


10 






2 


I 


6 




0 















NS = 



We interpret the entries in matrix A as follows. 
Starting in state a 2 there is a probability of | of 
absorption in state a 2 , and a probability of ^of 
absorption in state d 4 . A similar interpretation is 
made for the other entries. 

In the original example of Section 2.12 tfll of 
the matrices I, 0, S/aniT turn out to be square. 
This'will always be true for I and T, but is not 
generally the case for .0 and S. Consider the matrix 
below along with its -standard form 



The standard form 



0 .8 



"3 
0 

.4 
0 
0 
0 



U 
0 
0 

. 2 
, 2 
0 
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a 1 a 
1 ' 



a 2 a 3 



0 1 1 0 0 

.5 .1 | a. .4 . 

0 .1 ( .7 0 



0 



0 



.8 



"4 
0 
0 
0 
. 2 
. 2 



and we label each new, matrix as J^llows 



0"" 








0 


0 


0" 










0 














1_ 








_ 0 


0 ^ 


q_ 






.r 








" 0 


.4 


<T 






.i 




T 




.7 


0 


.2 






0_ 








-r 8 


0 


.2_ 




* 



0 



Question: If the matrix A had" come from the 
charge account problem, how would 
you interpret the entries? 

If a customer was in state a 5 , we 
f '-would predict that there was a 
probability of 5/6 that he would 
end up in state aj (which might be 
the state "paid up 1 '), and there was 
a probability of 1/6 that he would 
be absorbed in state a 4 ("bad account") 



Answer: 



2,16 Experiment 4 

Form an absorbing Markov chain from Table 1 on page 15. Put 
the matrix in standard form and then partition it as shown in 
Section 2,13* Use Program 8 in Appendix A to calculate the matrix 
N, and Program 1 to find the row sums of N and the product N x S . 

Write a report on the informat io4 you can give the Homer 
Department Store as a result of these computations. 
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2,17 Model Exam (Unit 111) 

I. *Show that this matrix is a regular matrix. 



0 1 

1 4. 
7 ¥ 



Find a fixed-point for the fol lowing matrix. 



3 1 

7 7T 



If P in question 3 is raised to. the 100 th power „ what is 
the approximate value of the entry in the first rod, first 
column? V 

<* 

Joe, as a student, is not very regular in completing 
assignments. ^However, if Joe is late with an assignment 
on one^lue date, he i s 70% sure to have the next one in 
on time. If he finishes an assignment on time, there is 
only a 20% chance that he wi 1 1 finish the next one on 
tftne. In the long run, what percent of time does Joe 
miss due dates for his assignments? 

Does the following matrix represent an absorbing Markov 
chain? Give the reason for your answer. 



0 
0 
0 
.7 



0* 
.5 

0 
.3 



Put the following absorbing chain.in standard form. 



"4 i- 



*1 

; i 
A 

.2 
0 



"2 
0 

.2 

.3 

0 



0 
.5 
.3 

0 



"4 

0' 

.2 
.2 
1 



Which are the absorbing states in this chain? 

i y * 
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If we start In. state a^ f how many steps will there be on 
the average before absorption? 

What Is the probability that if we start in state a^ f 

absorption* will occur in states,? 

4 * 



\ 

< 



* 
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3. ANSWERS TO EXERCISES (UNIT 107) 



a 2 tP a l -->23 " '- 1 



a 2 t0 a 2 a = " ^_ 

a 2 to a 3 z = P*3 * ' 3 



■jLi°, a j • ; " p 3 i m 01 



a 3 t '° a 3 " P 33 " M 
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4T ANSWERS TO MODEL EXAM (UNIT 107)' 

1. \ Each element, a^ t represents the probabi 1 i^ty 4 that a process 
• which starts In the i th state will go to the j th state In 
one step, j can be equal to i. 



Each row in a Markov matrix* represents *the probabilities 
for all possible next states. The sum of these probabili-L 




The Markov chain is: / 





a l 


a i 4 * 


a 4 


V 




















a l 


0 


0 1 


0 










"0 


0' ' 1/2 


1/2 








a 3 " 


1/3 


Ul 0 


1/3 








a 4 


0 


J^3* .1/3 


0 

o 

% 


1/3 






ft . 


0 


0 0 


0 


0* 







The probability that* the rat leaves.. the maze in three 
trials iitol/9. Irtdeed, from the tree diagram above, the 
bnly^ssibillty-for the rat to leave the maze is to travel 

through the branch rf- f^ — a /™'a fi *-. By the* compound 

^robabPlity rule, the probability of thfs 4 event is r t^ie product 
* ' * * ' ' ' 25 




V i * 

of probabilities from flj to a y a 2 to ^, to a 5 - Thus 
the probability that the rat leaves 0 the maze is 1-1/3*1/3 
1/9. 
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5. ANSWERS TO EXERCISES (UNIT 111) 
Reason 

absorbing There is an absorbing state, an<J 

there is a way to reach it. 

There is.no absorbing state. 



nonabsorbing 
absorbing 

absorbing 



There is an absorbing state and 
way to reach i t . 

There Is an- absorb i ng state and, 
way to reach it. 
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6. ANSWERS TO MODEL EXAM (UNIT 111) 



P> = 



1 3 



1 3 



3 9 
1Z TS 



Since a power of P is pos i t i ve, * the matrix is regular. 

7 



2 r 

3/~3 



2 
3 



/ 



Let represent "assignments on time" probability, and 
represent "assignments late" probability. Then 

*1 S~i \ 



c 



.2 

.7 
L \ 



.8 
.3 



is the Markov chain for this problem. Solving for the fixed 

point for this matrix, wc find a = ^ , and a n . The 

long run probability that Joe's assignments will be late is 

53.333% (~). 

1 -> * * * 

The, matrix represents an absorbing Markov chain. It has two 

absorbing states, a 1 and a . 

it o * 



1 
0 
.1 
.2 



0 
1 

.2 
.2 



"3 

b* 
o 

• 5 
.3 



States cij and are absorbing. 
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N - (I - T)" 1 - P 8 --5T 1 a p. 70 1.221 ; 

. , L--3 -7 J L -73 I.95J 

The number of steps before absorption, beginning in a 3 is 
approximatgTy jfi p 1^05 = 2.68. 

N x S P -70 1.22\ r.l .21 ' VMk .5841 

, L-73 1.95^ L-2 .2 J [.463 .586 J ' 
The probability is 58.4% 



f artral"/Bsv< er for Expef4ment 4 
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0 * 
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e 


^.25 


. 0.0- 
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, .15 
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0.0 


0.0 


0.0 
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,20' 


0.0 


.20 


.35 


.25 
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0.0 


0.0 


. 2 


.13 


0.0 


.10 ! 


.20 


.30 


.27 


0.0 


0.0 


3 


.10 


0.0 


'.05 


.10' 


.20 


.18 


.37 


p.o 
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.15 


0.0 * 


.02 


.03 


.07 


.30 . 


.28 


.15 
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.06 


.15 


.01 


.OA 


.OA 


0.0 


.25 • 




1 - 


"1 
_0 


0~ 

1__ 


0 f 


~0 
_0 ' 
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0 0 
0 0 


• 

0 
0 


0_ 



.25 


0.0 




.*0 


.15 
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0.0 


•0.0 


0.0 


.20 


0.0 




.20 


.35 


.25 


0.0 


0.0 


0.0 




0.0 


T « 


.10 


.20 


.30 
.20 


.27 


0.0 


0.0 




0.0 




.05 


.10 


.18 


.37 
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.15 


0.0 




.02 


.03 


.07 


.30" 


.28 


\M5 


.Q6 


• 15 _ 




_- 01 


.OA 


..OA 


0.0 


.25 


.A5 
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. . ' APPENDIX A 

. PROGRAM 1 . t )■ 

L*»»*»APPLICATIUNS OF MATRIX METHODS - PROGRAM 1 

C****#THIS PROGRAM PERFORMS MATRIX MULTIPLICATION 0*N ANY NUMBER OF PAIRS 
C4; - OF MATRICES < 
♦ DIMENSION A (20f 20) f B(20 1 201 f *C ( 20 » 20 ) t IHEADt 40 ) •/' 

C NR, '= LOGICAL NUMBER FOR THE CARt READER. '* ' 

"HR = °2 > • - 

L»«*'»«NP = THE LOGICAL NUMBER FOR THE PRINTER. M v 
* 4 NP__ =^5. ' 



C\»»»»READ .IDENTIFICATION CARD. - LAST CARD SHOULD CONTAIN /-* IN COL.' i-2 
'100. READ { NR , 1 ,E NU=7'0) I HE AD . 

I FORMAT 140A2 ) ■ * *. * 

/ WftlTE(NP,3) I HEAD * ' , . 

3 FORMA TUH1,40A2A/) * * . ✓ 

C»»*»»KEAD THE NUMER OF ROWS AND COLUMNS IN THE FIRST MATRIX.' 
READ { NR » 51 NROW, NCOL 

FORMA T { 2 I 5> I . 
*»*£AD THE FIRST MATRfX. 
00 10 I = 1 ,NROW 

t -' ' READ CNR » 7 ) ( A ( I , J ) J= 1 , NCOL ) - -.- 

• 7 ,FORMA T { 1 0F5. 0) ' 
, 40 . t . CONTINUE 

.C*»* *».READ THE .NUMBER. OF ROW$,AND COLUMNS IN THE SECOND MATRIX. 
READ ( NR i 5J MROW,MCOL - • 
, • DO 20 .1 =*1 ,MROW «*' * * • . wf- 

C t #i»#*READ THE SECOND MATRIX. # • • 

20 READ'( NR ,7) (BUfJ-)tJ * 1,MC0L )< 

C»»##»CALL THE SUBROUTINE TO MULTIPLY THE MATRICES. .. 
, : CALL' MA TMY ( A » NROW » NCOL » B » MROW » MCCL » C ) 

C#»»*» WRITE THE INPUT MATRICES. . • * * ♦ 

« . -WRITE <NP.,25). * * i 

25 FORMA iTSX, 'INPUT MATRICES'//) ^ » 

DO 35* I =1 ,NROW ' - - * 

WRlTE(NR,30r (Ad , J) , J - 1,'NCOL) ' . ' - 

30 . FORMAT{5X,ft)F10.3) * * - - « 

35 CONTINUE* • '*.<• ' 

• WRITE (NP, 40) 
AO ' 1 FORMATt///) 

00 50 I = 1%MR0W . . ~, ' ' • 

50 'WRITE(NP,30). (B(I ,J) ,J-= 1,MCCL) 
C»*»»*WR.I TE' THE MA TRI-X "PRODUCT. '* 

WRITE (NP, 55) ( ( 
55 F0RMAT(///5X, 'MATRIX PR0DU6J'//) 

DO 60 I = 1 ,NROW * 
60 WRITE(NP,30) (C ( I , J) ,-J= 1 ,MCOL ) 
C*»»»*RETURN FOR A NEW PROBLEM. 

GO TO 100 
70 CALL EXIT 
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PROGRAM 6 r -~ • 

C«#««APPLICATIONS OF 'MATRIX METHODS - PROGRAM 6 

OIMfcNSIBN A (20» 4) » B (20 ,4 ) , C ( 20 , 4.) , IHEAD( 40 ) , OBL ( 20, 4 ) ' • < 

NR = 2 

NP = 5 ' * 

C*#«»READ, PAGE AND WRITE HEADING 

READ(NR,1) IHEAO ' 
i w FORMAT (40A2 ) 

WRITE(NP,3) IHEAD . . ' 

3 ^ FORMAT ( 1H1 ,40A2//) 

C*»«.»CREATE OBLIQUE TRANSFORMATION MATRIX FOR PLOTTING. • 

DO 4 J = lU 

4 OBL (I » J) =0.0 
0BL(1,1J= 1.0 * 

0BL(2,2) = 1.0 ' > 

0BL(4,4) = 1.0 ' 4 " . ' , 

. 0BL(3,1 ) = -0.4333 
5 0BL(3,2) = -0.2500 
C*»*»READ THE NUMBER OF POINTS » MAX = 20. - ^ 

RE AO (NR ,5 ) N . * « v 

5 FORMAT (1 2) '. ; 
'C***»REAO THE COORDINATES OF THE POINTS CN THE FIGURE. 

00 10 I = 1,N ' • ■ " * 

READ (NR ,6 ) (All, J), J = 1,4) 4 • l * 

6 » FORMAT (4F5.0) ' 

10 . CONTINUE „ 

C»»«* THIS TRANSFORMATION fS CARRIED OUT TO GIVE AN OBLIQUE PROJECTS 
C^.ALSO, THE Z AXIS IS FORESHORTENEO. UHLiyUE PROJECTION. 

CALL MATMY(A,N,4,0BL,4,4,C) 
WR1TE(NP,15) 
15 FORMAT (5V, 'ORIGINAL FIGURE* ,// ) 

CALL KPLOT(C,N,1.0,0,1) ' 
C«**««READ, PAGE ANO WRI TE 'HEADI NG* ( 
100 READ,(NR,1,END=30) IHEAD 

WRITE(NP,3) IHEAO * 
C»(»*«REAO THE TRANSFORMING MATRIX 
1 DO 20 I = 1,4 * 

READ (NR ,6 ) (B (I , J) ,J = 1,4) 
CONTINUE 

CALL MATMY(A,N,4,B,4,4,C) 
CALL MATMY(C,N,4,0BL,4,4,B) , k ° 

CALL KPL0T(B,N,1. ,0,1) 

GO TO 100 - ' • 

30 CALL EXIT 



I- 
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•C **#•* PROGRAM 7 APPLICATIONS of. MATRIX METHODS 
• C **»*t POWERS OF MATRICES ' 

• DIMENSION *A (20,20) ,B<20-,20 ),C(20,2Q), IHEACUO ) 

NR = 2 . . ' 

NP = 5 * '. 

READ (NR, 5) I HEAD • .* 

5 . FORMA T( ( 40A2* ' * •* ' V 

WKITE(NP,'6.)IHEAD 

6 - FORMAT(IH1,40A2-,//) ' 

C*.»»» READ THE NUMBER OFC ROWS AND CotUMNS IN THE MATRIX. 

C»«#»» ,M = THE POWER TO 'WHICH THE MATRIX IS TO BE RAISED.' 
_ READ (NR, 101 NR0W,NC0L, N ' 

10 ' FORMAT (31 f> J . 7 ' ^7 

DO 20- I '= .1 ,NROW v ; *' ' ' 

KEAD(NR,lf>) (A(I.,J) ,a,= t»NC'OL) . . 

li> ^0RMAT(F6F!>.0) . 

fO CONTINUE ' s " '*«5r . , 

' Pl3 # I = 1 Jnrow , " "' . • 

DO 30 J I ,NCOL v * 
30 . 6 (I , J) = A(I r J) 

DO* 50 K = 2»N ' \ ."J* 

CALL MATMY ( A 

DO '40* I 

DO 40 

40 A(I ,J) = C(I ,J) * * 

i>o continue .. ; , d 'fr * 

% • DO 70 I = 1 ,NROW ■ ' j 

WRITfr(NP,60) ( A tl , J) ,J=1,NC0L) * 
60 FORMA T (5 X , 1 1F8. 4)' 
10 t CONTINUE 

CALL E XI T 

6NU 



K = 2,N e . . ■ ». * 

1ATMY(A,NR0W,NG0L,B,NRCW'.NCdtv|C) A 

I I ,*IROW , A * »%] 

j = i.NcoL • * ' . V* 



r 
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PROGRAM 8 



C*»*»*APPLICATIONS OF MATRIX METHODS - PROGRAM 8 

DIMENSION A<20,20) ,T(20,20) ,C(20,20), IHEADI40) 
NP = 5 

NR = 2 » 
C*»*»*READ HEADING AND IDENTIFICATION &F INVESTIGATOR, 
1 READ (NR,I0»END=50) IHEAD 

10 FORMAT(40A2) 

WRITE (NP, 20) IHEAD 
20 FORMATdHl ,40A2//> 

C*»»»*»READ THE DIMENSION 0F C^UST BE SQUARE 

READ ( NR ,25 ) N ( 
25 FORMAT (I 2) 
C*»*»»»READ THE ARRAY T 

DO 35 I = 1,N 

READ(NR,30) ( T( I » J) , J=I ,N) 
30 FORMAT(10F5.0) 
35 CONTINUE 
C»»»»»(-ORM IDENTITY MATRIX 

CALL IDN(A ,N) 
C»*#»CAbCULATE . I - T. : 

CALL MA f SB (A,T,N,N,C) 
C»**»CALCULATE THE INVERSE^flF I - T. 

CALL INVER(CN) - i 

DU 45 I =--l,N 

WKITE\NP,40) (C (I , J) ,J=1 ,N) 
40 FORMAT\5X,10F10.4) 
45 CONTINUE 

GO TO/ 
50 CALJ^EXIT 

El 
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SUBROUTINE MATMY 



20 
15 



10 



'SUBROUTINE MATM* ( A , NROW , NCOL , B , MRGW , MCOL ,K) 
DIMENSION A(20,'20) ,B(20,20) ,G!L20,20) 
NP = 5 * ^ 

IF(NC0L-MR0W)20*10,20 
WRITE(NP,15) 

FORMAT (5X,«NUHeR OF COLUMNS IN FIRST MATRIX MUST EQUAL NUMBER OF 
IROWS IN SECOND MATRIX') 

„ CALL EXIT * 



00 40 I 
DO 40 J 
C (I ,J) » 



*• 1 t NROW 
■ I tMCOL 
0. 0 



00 40- K = 1 , NCOL 

C(I,J) = C(I,J) + A ( I f K ) *B { K f J') 

RETURN' 
END 



SUBROUTINE MATSB 



J 



t 



♦ONE WORD INTEGERS — 
SUBROUTINE MATS3(A, 3, WOW, NCOL , C ) 
DIMcNSlUN A(20,20), 6(20,20), C(20,20J 
DC 10 1=1, NROW 
DO 10 J=r,NCOL 

10 C(I,J) „= A( J,J)-8( I,J)" 
RETURN r 
CND » : 



r 



SUBROUTINE IDN 



♦ONE WORD INTEGERS 
♦LIST ALL 

SUBROUTINE IDN ( A » M ) 

DIMCNSION A(20,20) ' 

DO 20 1=1, M 
. DO 10 J=1,N 
10 A(I,«J) = OiO 
20 A J I , I ) ='1.0 

RETURN 

END 




SUBROUTINE INVER 




20 



♦LIST ALL 

*CNE W.CRD 'INTEGERS 

SUBROUHNE INVER(X,N) 
DIMENSION X(20,20), A120,40) 
DO 10 I = i,N 
DO 10 J = i,N 
A(I,J>J = X( I, J) , 
M = J + N ** 
IF(I-J) 20,l5t20 
A( I,M) =1.0 x " 
GO TO 10 , 
A( I ,M) = Q*0 

ML". C ON TINUE 3__ 

^0 55 K a l,N* * 
PIVOT = A( K,K ) 
IF(PIVOT) 35,30,35 
WRITf(5,101) . < 

FORMAT {///,5A»ZCR0 PIVOT*) 
CALL r EX I T 
A { K , K ) = 1 . 
IR = K + l. 
M = K+N 

DO AO J = n,H 
A ( K, J ) =. A( K, J )/P IVOT 
00 55 I • U N — 
IF(I-K) 45,55,45 
PfVOT = A'( I « K ) 
A( I,K ) = 0.0 
00*50 J = IR,M , 
MI, J) = A{ I, J) - PIV0T*A(K, 
CUNT I NO C 
IR ? 2*N 
K a ,M + L 



30 
ICl 

3 5 



40 



45 



50 
55 



it* 



60 



Du 60 I 
JO 60 J 
f = J-N 
X ( I ,M ) = 

RETORN , 

end 



= l,N 
= K, IR 

Al.lJ, J ) 
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SUBROUTINE KPLOT 
* 



ICHR/ • A 1 , • B» , 



•C I ,».D», , F«, 4 F», 



* H • . • I 



SUBROUTINE KPLOT (C, IROW,S, IR, IP ) 
C*****THIS SUBROUTINE IS THE SAME AS IPLOT 
INTEGER PLANEI41, 71 ), ICHRI20) 
REA1 C(20,2) 
LATA IPLK,IX/« «,♦.»/ 
DATA 
1 •N% , 0« 

S IS, A SCALING FACTOR TO BE USED 
RANGE IS FROM -20 TU +20 
SET' IR = 0 TO BL ANK * OUT GRAPH FRAME 

PUT ' NEW GRAPH IN WITH PRfcVIO'iS 
SUPPRESS PRINTING OF THE GRAt'H 
PRINT THE GRA"*H 



EXCEPT T^A f IT SETS UP 3 AXES. 



c***#* 

c ***** 

c***** 



SET- 
SET 



C***** SET 
C***«* mp 

.«P = 



IR 
IP 
IP 



1 
0 
1 



TO 
TO 
TO 



If COOKDl nML r > ARE 
U <E. - 



OUT 



OF RANGE 



\ 



NUMBER FUR PRINTER. 



J 
6 
B 

10 

2fc 



25 

30 
15 



35 
36 

37 
39 

40 . 



IF(S - 1. ) 5,6,5 v 
CO 7 1=1, IROW * 

DO f J = 1,2 
C(I,J) = C( I , J ) * 
I F( IK) 15,8, 15 
DO 10 I = 1,41 
00 10 J = 1, IV 
PLANCH, J) = I3LK 
UU 20 J = 36,71,2 
PLANE(2I,J) = IX 
J = 34 

00 25 I = 1, 10 
PLANEl 1+21, J ) = IX 
J = J-2 

UO 30 I = 1,20. ' 
PLA.mE I 1 , 36 J = IX 
DC AO K = 1, 1 ROW 
J = 36 + (5./3. * C(K,1) + .5) 
IF(J) 40r40,35 . 1 

1 F ( J-7 1 ) 36 , 36 , 40' ^ 

I = 21 - <C(K,2) - .5) 
IF< I } 40,40,37 
IF( 1-41)39, 59,40 . * 
PLANC(IVJ) = ICHR(K) 
CONT INUE 
RESTORE MAIRIX C IF 



r 




IT HASBlEN SCALED. 



42 



44 
4i 
45 

50 
60 
65 



IF(S-1. ) 42,43,42 

SI = l./S 

DO 44 I = 1, IROW 

DU 44 -J = 1,2 

C( I, J) = C( I, JJ *SI 

IF(IP) 45,6S,45 * , 

00 50 K = '1,41 ' 

wRlTE(NP,-60) ( PLANCCK, j') f J =r,7l) 

FORMAT i 5X ? » 7 1 A 1 ) 

RETURN • ' '. 

END * - 
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STUDENT FORM 1 
Request fen; Help . 



\ 9 "i 



'Return to: \ * 
EDC/UMAP 
55 Chapel St. 
Newton, MA 02160 



Student:, If you have trouble with a/specific part of thia unit, please fill 
out this form and take it' to your, instructor for assistance. The information 
you give will help the author to revise the unit. 



Your Name 



Page 

0 Upper 
^QjWiddle 


OR* 


Section . 
Paragraph . 


• 

OR . 

> 


0 Lower 









Unit No. 



Jlodel Exam 
/ Problem No , 

'Text • o 
Problem No, 



Instructor : Please indicate your resolution of the difficulty ih this box. 
Corrected errors in ^materials. List corrections here: / " A 



o 



Gave student better explanation, example, or procedure than in unit, 
Give brief outline of your addition here: 



o 



Assisted student in acquiring general learning and* problem-solving ^ 
skills "(not 'using^ examples .from this unit.) \ 



88. 

• * * 
' Instructor's Si£natu$e_ 



Please use^ reverse if necessary. 



Return to: 

STUDENT FORM 2 EDC/UMAP 

Unit Questionnaire 1 J 5 ? hape * S *' „ 

flewton, MA 02160 



Name Unit No, Date 



Institution s Course No, 

Check the choice for each question that comes closest to your personal opinion. 

1. Hpv useful was the amount of detail in the unit ? 

Not enough detail to understand the unit 

Unit would have been clearer with more detail 



Appropriate Amount of detail 

Unit was* occasionally t N oo detailed, but this was not distracting 

T oo much detail; I»was often distracted 

♦ * 
2. How helpful were the prdblem answers? 

JSample solutions were* too brief; I could not do the intermediate steps 



^Sufficient information was , given to solve the problems , 
^Sample solutions were too detailed; I didn't need them 



3. Except for fulfilling the prerequisites, how, much did you use other sources (for 
example > instructor, friends, or other books) in order to understand the unit? 

A Lot Somewhat A Little Not at all 

4. How long was this unit in comparison to the amount of time you generally spend on 
a lesson (lecture and homework assignment) in a» typical math or science course? 

Much ^ Somewhat About Somewhat Much 
Longer Longer the Same Shorter Shorter 



5 • Were any of the - following parts of the unit confusing or distracting ? (Check 
as many aa apply.) 

Prerequisites 



S tatement ofskills and concepts (objectives) 

Paragraph headings 

Examples > n 



^Special- Assistance Supplement (if, present) 
.Other, please explain 




6. Were any of the following parts of the unit particularly helpful? (Check as many 
" as apply 

. P rerequisites <T > 

Statement of skills an^ concepts (objectives) 

3 Examples ■ * 



_Problems 

^Paragraph h,eadings v . 4 
JTable of -Contents' , 
^Special Assistance Supplement (if present) 
Other, please explain 



Please describe anything in the unit tffkt you did not particularly like. 



Please describe anything that y$u found particularly helpful. (Please use the back of 
this sheet if you need more space. )• * 
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1. ELECTOItAL CIRCUITS OJ108) 

l.lp Introduction * - * 

You may be awar.e that matrix methods play an important 
part in solving systems of linear equations.' We will 
examine a few aspects of this problem which aVe treated 
more completely elsewhere/ We turn our attention to the 
way in which a system of'equations migfit arise^ from a , • 
simple problem in physic* dealing with' an electrical 
circuits * * " 

Most people have a general notion of what is meant 
by an electrical current flowing in a wire. The, flow 
of Electrons in a wire is somewhat like the flow of 
water in a pipe. To produce a flow of current, some 
source of power is needed, such as a battery/ There is 
also a part of a circuit which consumes power. This is 
a resistance. Current is measured in amperes , the source 
of power in volts and resistance in ohms. 

1,2 Laws for Electrical Circuits * 

In studying circuits we will use' three laws. We 
first state^ the laws, yien .show how they are, applied. 

1. The sum of alll^rrents flowing to a 
point equals ,the sum of the currents 
flowing away from the point. 

2. The algebraic 'sum of the voltage drops 
around any loop of a circuit is zero, 

3. The voltage <tr op between two points of 
a circuit* algebraically equals the 

- product of the current and the resist- 
ance between the points*, 1 

*To illustrate these laws we use the cifcuit shown 
in Figure 1. The symbol \/\J\ represents a resistance. 
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10 amps 



3 ohms 



5 ohms 



> VNAA/V 



Figure 1. 



Ten amperes of current are flowing ffto point A; 
therefore 10 amperes must be flowing away from point A. 
If and % 2 represent the currents in the upper and lower 
branches of the circuit, respectively, then 

i x + i 2 » 10. (1.0) 

Across the resistor in the top branch there is a 
voltage drop, which we will lafiel * r The thir£ law tells 
us how to calculate this voltage drop from the values of 
the resistance and the current; ' 

Similarly, in the bottom branch 
, e 2 = Si 2 . 

Finally, we note the opposite directions of the current 
at point B and apply the second law, to obtain the > 
equation 



3 *1 " 5 *2 = °* 



(2.0) 



1.3 Solving Linear Systems Using an Inverse Matrix 

To find and % 2 in the circuit shown in Figure 
1, we now have the system of equations 

i 1 + i 2 " 10 . (1.0) 

(1) 

0 „ (2.0) 



A 
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which can be solved easily by substitution. The solution, 
can also be expressed in terms of matrices. The mati^x tf>ft> 
of coefficients of the system in fl) is * * 



/ 



1 1 



3 -5 



If we write the unknown values as the column vecW)r>* 



*^^and the constants on the right hand side of the equations 

% rioi 

*"\ « • B = 

S 

then the system of equations can be' written as 

, (2) AX = B. ' 

One way to solve for X is to multiply both 'sides of 
Equation (2) by 'the inverse matrix A -1 ^ 4 to obtain 

A" 1 AX = A^B 

$ IX = A^B 

X = A^B. / 

The matrix I in this calculation is the identity matrix/' 
Since the inverse .of A in our particular case is 



.-1 



we have 



A" 1 !* 



5 


1" 


8" 


T 


3 


-l 


jS 


8" 


~S 


,l" 


8" 


B" 


3 









10 



25 
T 



15 
T 
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From this we conclude that 



25 



15 



t l * T amps and i 2 = X am P s « 

1.4 Consistent and Inconsistent Systems of Linear 
Equations 

A system of linear equations is safid to be consistent 
if there is at least one solution for the system. An 
inconsistent system has no solution. To illustrate this 
geometrically for a system of two linear equations in two 
variables, we may represent each equation by a straight * 
line, as in Figure 2. - * • , 

1.5 Existence Theorems 

There are theorems about systems of linear equations, 
called existence theorems, that allow us to determine 
whether a system has no solution, a unique solution, or 
an infinite number of solutions. As you continue your 
study of linear algebra, you will learn about thesl * 
theorems ^and their proofs. Program 5 in Appendix A applies 
these theorems to systems of linear equations . You can 
use this program to .investigate the nature of the systems 
igith which you work even though you have not studied the 
theory. 
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Consistent, one point in common, 
uniqUe solution. 




1 ♦ 6 An Examp 1 e 

Let us apply the three laws for electrical circuits 
to the following circuit* 




Using Law 1 we have the following 
For points : A 
B 

e 

D 



% 1 + 

i 2 + h 



*3 



Using Law 2 we have: 

* For- loops : & 

ABC 

* ^ * BCD ' 9i 



6i 1 + 2i 3 --3i 4 



ABCD 6i 1 + 9t 2 - 7i $ 



3t, 



We can arrange these equations in the form AX = B. 

Since we have £i>Te variables and seven equations the 
system seems. to be over-determined. A system of equations 
*is said to be oVer-determined when the number of equations 
I 

P y 6 



is larger than frhe number of variables involved in the 

equations. It is possible to mdtfce substitutions which 

eliminate two of the equations. However, it is^not 

necessary to do this. In many practical problems, 

« 

systems 6f equations which are derived £rom physical 
situations may consist of 50 or even IOjO or rflore * 
equations. One of the benefits of matrix theory is 
that we can use it to find out whether an over-determined 
systelh is consistent, a*rtf*wheth4r or not the .system has 
a unique solution. t. 

1.7 Experiment I 

Consider the system 

2x + y - z = S 
x i 4y '* 3 

^ Sx - 2y - 2z =*13 ■ 

. 6x - Sy - 3z = , ' . ^ 

Is the system consistent? Is there a unique solution? j 

1.8 Model Exam for Unit 108 

1. A system, of equations which has more equations than unknovtfi 
variables is called * 

> : * 

2. If a system of equations has two equations and two unknowns, 

the system is if the graphs of the 

equations intersect, the system is 

if the graphs are identical , and the system is 

if t^ graphs do not intersect. 

^. Consistent systems may or may not have 
solutions, 

4. Find the values of iy and for the electrical circuit 
on the following page. 
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/s 



1 ohm ' 



5 ohm 



^ 2 ohms 



< > ' 




V * / 





\ 



Figure for Model Exam Problem 4. 



.Use Program 5 to investigate the nature of the solution, if 
it exists, for the following system; * 

2x - y + 3z + w - 1 
x + Ay - 3 + 2ii?=:4' 
5x + 2iy + 53 + 4u = 0 

Note that this is an "under-determined" system in the sense 
that it has more unknowns than equations, but it is still 
possible -to investigate it With our program. 



'80. 



ERIC 



2. APPLICATIONS OP MATRIX METHODS: 
ANALYSIS OF LINEAR CIRCUITS (Ull2) 



2.1 Intr oduction 

; — V- * 

In Unit lOfr on electrical circuits, we considered 
ho.w a system* ofVlinear equations could be used to 
represent some of ^the ^relationships *j.n an electrical 
circuit. Such a system of equations is, useful in linear 
circuit analysis'.* \ We then explored the use o f> matrix 
methods to^solve systems- of xaguat ions , which,* in some- 
cases,, were over*- determined. We also referred^to the 
fact that systems ^of equations can*J>.e <onsistent or in- 
consistent, and can \ave^a unique solution, an infinite * 



number of solutions, or no solutio 
simple method to solV 4 e such systems^ 

2 T 2 Efementar^ Row Operations 



We. now explore a 



I- 



We are interested Vn*the following three types of 



elementary row operations which may be performed on a 
matrix: A ^ * 



fo^me 



1. the interchange of any' two distinct rows; \ 

2. the multiplication of any row by a nonzero, 
scalar ; .„ \ • • - J 

3. the addition of a scalar multiple of one rdw 
o± a matrix- to, some oth^n row of the same '<> 
matrix. * ... 

If you think of ordinary linear equations , Uhese are the. 
usual ways in which you manipulate them. ' 



*A linear circuit for a direct current is, one containing ' 
elements ^hat obey Ohm's Law, such as metallic conductors.' Ohm's 
Law was given as fc^e third law of electrical circuits on page 1> 
There ajre many devices in electronics that do not obey Ohm's Law. 
They are cfalled "nonlinear." ^ 



x 
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2.3 Exercises Using SLejr.entary Row Operations 
1. Which elementary row pjterakion transforms 
1 8* 



1 8 



to 



0 • -17 



Transform 




-1 

7 



to 



1 



if possible. r v 

2.4 Row Equivalence ^ / ^ 

V matrix is said to be row^equ-iva^ent to another 
matrix if the f irs£ matrix, can be transJormed into the * 
second by a sequence of eleirfentary row operations. .Fo^ 
example, in 'the exercises in Section 2v3 the matrix 




♦ is .row e^uivalerit to 



is row equivalent to 



n 



-17 
0 
1 



ami the matrix 



^ Why we are interested in j^ow equivalence will be 

evident, in the discussion which follows. -We will see 
that a transformation of a matrix of coefficients of a 
system of linear equations which leads to a particular 

^ row equivalent matrix is a means of 'obtaining solutions 
for the system of equations , ;if they exist*, or in' 
determining that the system is inconsistent. . 

2.5 -Row Echelon Matrices ^ , 

i • <& 

If you look ^p the term, echelon in a» dictionary , v you 
■ will £ind that it refers to arforirfation, often us,ed for 



V 
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c 'V; 



\ < 



(A % 



X 



airplanes'or ships, in which there is a lead plane or • 
"vessel with the others arranged in step-like fashion 
slightly to the right or left and to the rear. We use 
the term echeion .here to refer to matrices of a form 
that, in a way, suggests the meaning of the term as just 
^given. More precisely, we defiae a matrix to be in 
row Echelon form if it has the follow^itig properties: 

■ 1*. in any row of the*matrix the first nonzero 
* element >at the left must be a 1 unless the 

row consists of all zeros; N „ 

rows of all zeros should follow nonzero 'rows ; 
•the column containing the leading <i.e., the 
leftmost) 1 has zeros elsewhere* in the column; 
the leading 1 of any nonzero tow must appear 
*to the left of the leading 1. of the nonzero row 
that fojlows it/*" 



3. 



The" following examples should help to cla*rify this 



idea: 



1. 




1 


0 


*2 












r- 




















8- 


1 


0 




t 


& * 








0 ' 


0 


0 






^ techelon 


form) 


* 2. 




1 


• 0 ' 


2 














1 ■ 


0 


0 


• ✓ 












0 


0 ■ 


0 




* 


tnot in echelon 


form) 


3, 




1 


0 


0 


- 3 




* • * 








0 


0 - 


1 


2 












0 


0 


0 


' 0 




, C e chelon 


.form) 


4, 




1 

0 


% 0 ; 
1 


0 

. '6 










*? 






2 


■ l 






(not in exhelon 


form) . 

i 














t 

* 
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2.6 Using Flow Echelon Form to Solve Systems of Equations 

For a given system of linear equations, the 
coefficients of the variables and the constant terms 
A can be represented as two matrices. F.or example: 

2x - ly = ?6 / 
Sx + y * -3 





2 


-3 




6 


V » 






and C - 






5 


1_ 




^3 



If we write the matrices V and C as one matrix, by 
writing the constants as a new column on the right: 

[vie] , ' 

i 

and trans ft>rm this augmented matrix- to echelon form, we . 
can find "Che solution of th*e original system of linear 
equations represented by V and C, if the* solution exists. 
Further, ,this method will expose inconsistent systems, 
and systems with many solutions. We will demonstrate 
^ this with an- example. For the above -system 



.-3 



multiply tqw 1 by j,. 



3 
'7 



multiply row 1 by -5. and add to row 2 



4 : 



1 


3- 


. 3 


0 


17 
T 


-18 



: j 
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v 



imultiply row 2 by ^ 



n 3 1 

1 T 


3 






0 1 V 


. 36 


_ s. 


"T7 



multiply 'row 2 by j and add to row 1 to achieve row 
echelon form * 



1 , 


0 


, 3 
"17 


0 


1 


36 
"T7 



The solution for this system can be read from ,the 
row echelon form as 



and 



".3 
•T7' 



36 
"T7- 



) 



i From this problem we can see, that ^transforming a 

•matrix of coefficients to row echelon form, an'd at the 
*. same t.ime applyi*ng the$row operations to the 'augmented 
, matrix of coefficients and constants can produce- the^ 
solution* to, \the system, if it'eys.ts. , Sp 

2>7 Examples. and Exercise's 

Use the method of transforming thfe augmented matrix to. row * 
echelon form to sol ye each of these systems,, if possible. State r 
whether or not the system is inconsistent 5 . The first t example shows 
how this can be detected- * » ' 

a. x - y =3 • . / \ . 

• , 2x + Zy t + z - 1 . * , 
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Q . 9 . 



The 0 answer is 


1 


0 " 


-l 


2 




0 


1 


i 


-1 




0 


b 


0 


-2 



This system is inconsistent. Look. at row 3. There is no 
* value of x, y t or z such that 0 = -2. 



% b. x * 2y - z * 


* 4 








• /-x - 


1/ + 23 = 


-s 




* # 




4x ♦ Ify - 3 = 


14 








c. 3x ♦ 2y 


+ 3 = 5 










2x ♦ 5z/ 


+ 43 = 8 










x-+ 4z/ 


+.63 =4 










d. tr + z/ + z = 3 
x - y - s = -4 




• 






3x + y + z = 2 




ft 






Answer: 


1 0 


0 


1 

"2 




V 




0 1 
0 o 


1 
0 


" 2 
0 




♦ 



From the first row x = -i, and from the second row y V z * j. 

The^third revindicates £hat 0 = 0. This system is consistent, but 
there is no unique solution. We can, however, find a particular 
solution if we assign some value to i/ or 1 For example, let 2 - ^ 

.7 c 
then y ~ 2 For = y s f- •* • 

In this case, the system is consistent but- has an* infinite 
. «. . * • 

number of solutions. We refer to such systems as consistent systems 

\ With parameters, ^ntf the ^example above k is a parameter. A system 

may have more* tha'n one parameter. *. ^ * 




14 



e, x - y ».!*•. Does thisystysteia have a parameter? Why? 



3x - y = 3 

x - 2y = ~2 
x - z/ = -2 
2x + zy = S 

x ♦ 42 = S 

-2x + *y ♦ z = S 



•V- 




• 2x ♦ */ ♦ 17s = 13 

h. ' x ♦ ly - s = 4 
2x - zy ♦ 22 = 3 
2zy ♦ 2 = 7 

2.8 Electric Circuits Revisited 

In Unit 108 we used matrices to analyse an electric 
circuit.' The* progranr is quite adequate for even a large 
system. However, at the time, we**started by formulating 
O equations based on the laws of circuits and then 
constructed a matrix from this system. Tf a system is" 
more complicated, this may be a nontrival task. 

It is possible to formulate the etectric circuit 
problem in terms of matrices from the beginning without 
writing the equations. This will be illustrated^ in°the 
following example. 

» * ■ *# 

, 2.9 An Example, of an Electrical Circuit 

m * "* 

. Before we start the example, it is necessary to 

discuss, brief ly; network branches and sign convention^. 
Each branch of a network -£n an electrical circuit can be 
represented as shown'in Figure 1. 



— v r - Z>|... 

-Q- 1 — U/VVV^ ■ i 



Figure 1. 

The voltage drop across the branch is given by 
» r - where * r is an electromotive force in series 
with v T . Figure 1 shows the sign conventions used for - 
a branch. * 

Figure 2 shows an example of an electrical circuit 
usinf the conventions illustrated in Figure^. 




s 



. Figure 2 . 
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Each hranch is connected to the rest of the network 
at precisely two points f or nodes . We number these nodes 
in Figure 2, arbitrarily from 0 to,4\ We are interested 
in only the voltage drop across a branch 1 , that is, < 
differences such as V~2 - which is .dqual to v 2 " e i> 
so '.we can set the voltage at one arbitrarily selected 
node equal to 0. We choose V Q s 0. 

From the- laws previously stated in Unit 108, the 
net current at each nodte must be zero, arid we can write 
the equatioTTSP^or this circuit as 



I 



0 



I, 



-t 1 + t 5 + t 6 



l 3 + > l 4 



"1 

• h ; h ' U ' h - H * 

% H&wever, instead of writing the,se equations, especially 
if the system is" large afid complicated, - we can construct 
the matrix S whiqh. preserves ^ie si£ns of the system. - 
This matrix can be constructed directly from the diagram * 
without the 'necessity of writing the equations. 





i l 








*5 




. r l 


1 


0 


i 


1 


I 


0 


■ s " h 


0 


1* 


-l 


0 




^0 




0 * 


-1 


0 


-1 


0 


-1 


1 






V 


\ 







X 4 



> 17 



to 
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Verify that the equations "could be obtained from 
the product 

" * -* " 

J - St ^ » 



where 



From the ^iagram in Figure 2, we can construct the 
followi ng matri ces . * 

















or 


Tv l~ 






f 


°K 




e l 


V 
















V 2 


J ■= * 


h 




v 2 




e 2 


< : 






h 




V 3 




e 3 










v = 


V 4 




e '< 












W S 







































Since the law of electric circuits 

* v ■ Ri • 

holds for .the voltage drop acr.oss each resistor in the 
circuit shown in Figure 2, we have" c ' 



\ 



or, in matrix notation, 



v 

v l 




T> 

R l 


U 


u 


u 


U 


u 




£ 1 


V 2 




0 * 


R 2 


0 


0 


0* 


0 






V 3 




0 v 


0 


R 3 


0 


0 


0 






V 4' 




0 


0 


0 


R 4 


0 


0 






V S 




0 


0 


0 




R 5 


0 




*s 


V 6 




0 


0 


0^ 


0' 


0 






£ 6 



Satisfy yourself that this equation holds. 

All our' information is now organized in the matrices 
S, i , V, J t \v , e y and R. 

From the laws of electrical circuits the 'foi lowing 
relationships are true: 



v = Ri 

T 

J = Si 



(2.1) 

(2.2) where 
(2.3) 



indicates the 
transpose * 



where all the variables in Equations 2.1, 2.2, and 2.3 
represent the matrices constructed above. 

2 . 10 Experiment I . . 

Using Figure 3 construct £he matrices S and R for this 
system. - Let J = 0 and solve' for 4V. 



V 

r 
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1*2 = 1* 




Figure 3. 4 
The values for the e f s not shown are all zero. 



2,11 Model Exain foy. Unit 112 
1- Transform this matrix * to row echelon form* 

Y o ; 

-1 /-3 2 

2, Solve this system of linear equations, if possible.. 
x - y + Jjg = 2 
2x + 2y + 23 '= 2 ( 
3r. + z/ + S3*= -4 



,3. 



Is the system in Probjtfm 2 consistent? Does it have 
parameters? ' ~ , * 



20 



J- 



0 
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Xreate the necessary matrices for the analysis of the 

electrical circuit in Figure 4 such that the formula 

-1 T -1 
SR £ V = J - SR e can be used, Do not solve the system 



i , 




»R, 



^ ^(Figure 4. 
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3, ANSWERS TO MODEL EXAM (UNIT 108) 



5, 



Over-determined . 

Consistent, consistent, inconsistent, 

Unique. 

N b 4, M = 3,' 

AUGMENTED MATRIX 



1.00 
1.00 
0.00 
1.00 



1.00 
-5.00 
5/00 
0.00 



I. Off 
0.00 
'-2.00 
-2.00 



30.00 
0.00 
0.00 
0.00 



UNIQUE 'SOLUTION VECTOR X IS 

X(l) = 17.64 
X(2) = 3.52 
X(3) = 8.82 

N = 3, M = 4. 

AUGMENTED 'MATRIX 



2.00 


-I.OjO 


3.00 


1.00 


1.00 


l.do 


4.00 


-l.do 


2.00^ 


4.00 


5.00 


2.00 


5.00 


4.00 


p. 00 



\ 



EQUATIONS ARE INCONSISTENT. 

4. ANSWERS TO SOME EXERCISES FROM UNIT' 112 
i *■ — 

Exercises from Section *2. 3: 

The multiplication of row 1 by a -2 and addition of jthe 
results to row 2 will transform 



1 8 




i 

1 


* 8 




to 






2 -1_ 




0 


-J7 
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C ■ .7 > * 



To transform 



2 ~1 

1^7 



1 Q 



to 



0 1 



multiply row * by 2 •» 



multiply row 1 by j and 
add the results to row 2 



1 -i 
1 7 



1 



°-. 7 i 



multiply row 2 by — r— 
' 7— 

O n 1 / r 

multiply row 2 by j and ^ 
'add the results to row 110 

0 1 



5. ANSWERS TO MODEL EXAM (UNIT 112), 



The row echelon form is: 



1 
-1 
0 



2-3. 



Using Program 9 we g#t the ^following result 
ROW ECHELON FORM 

1.0000 o.oodb 
0.0000 : 1.0000 
0.0000 $0,0000 



0,0000 *0.0000 



2^000 
-1.0000 
0.0000 
0.0000 



1,5000 
0.5000 

o.ooc/o 

0.0000 



1 



From the row echelon form we cap see that the system 
consistent, but that it has parameters. 



9 
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Using Program 5 we can find particular solutions 
might want to discuss this result -in class, 
PARTICULAR SOLUTION VECTOR IS 
X(l) = 0.00 
X(2) = 0.*4 

X(3) = 0.75 ' : ; 

LINEAR INDEPENDENT VECTORS ARE 
U(l) to U(l) 

-1.00- r 
0.50 
' 0.$9 





























} 2 








0 




*3 


* 


V 2 




°/ 


i = 




V = 


V 3 


If* 


0 




e:J 




_ V 4 




0 










; i 


0 


t 














-r 


0 


' 0 


• 0 


0 ' 




0^ 


-1 


1 


,0 


-1 


S = 
















o — 


-r 


1 


0 






0 


0 


M 


1 






* 








and^ J 











0 

R 2 
0 

0 

0 

0 



1 

0 
0 

-1 



0 
0 

R 3 
0 

0 

0 



0 
0 
0 

R 4 
0 

0 



0 
0 
0 
0 
0 



0 
0 
0 
0 

R 5 
0 



'is 

r 
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APPENDIX A 

PROGRAM 5 4 ' 

INTEGER 'HEAD (40) \ , 

DIMENSION AA(10,10),BB(10),X(10),U(10,10) 
C****READ ONE LINE HEADING WITH SYUDENT NAME 
C**** A /* TERMINATES THE RUN ' * I * 
100 READ(2,20,END*50) HEAD- 

20 F0RMAT(40A2) " I 

WRITE(5,21) HEAD • 

21 F0RMAT(1HU5X,40A2//) • . f 
C****READ THE DIMENSIONS OF THE SYSTEM 

READ(2,1C) N,M 
,10 * F0RMAT(2;5) 1 

WRITE(5,33) N,M 5 
'33 F0RMAT(5X,'N = '.y.SX.'M = ' ,13//) 

WRITER, 34) 

34 F0RMAT(5X' AUGMENTED MATRIX'/) 
;D0 1 I = 1,N ' 

C****READ THE COEFFICIENTS AND CONSTANTS 

C****THESE ARE PUNCHED IN FIVE COLUMNS EACH WITH A, DECIMAL POINT. 
C****CHANGE THIS PROGRAM IF THIS FORMAT IS NOT SATISFACTORY 
C****THIS IS A SHORT CALLING PROGRAM AND CAN BE ADJUSTED EASILY 
READ(2,2) (AA(I,J),J=1,M),BB(I) 

2 F0RMAT(11F5.0) / 
WRITE(5,35) (AA(I,J),J= 1,M),BB(I) 

35 F0RMAT(5X,11F10.2) 
1 CONTINUE 

CALL SOLEQ(AA,N,M,BB,X,K,U) 
IF(K)100,41,42 

41 WRITE(5,5) 

5 , F0RMAT(//5X,' UNIQUE SOLUTION VECTOR X IS'/) 
GO TO 36 

42 'WRITE(5,43) „ „ , 

43 F0RMAT(//5X,' PARTICULAR SOLUTION VECTOR IS'/) 

36 DO 32 I = 1 ,M * 
WRITE(5,31) I, X(I) * 

31 FORMATISX/UMS,') = ',F8.2) t 

32 CONTINUE ^ * 
IF(K)100,10b,40 

,40 WRITE(5,6) ' t 
6" F0RMAT(//5X 1 LINEARLY INDEPENDENT VECTORS ARE'/) 

WRITE(5,7)K 
7 F0RMAT(5X'U(1) TO U( , I2, , )7) . 

DG 4 I = 1 ,M 

4 WRITER, 3) (UCI,J),J=1,K) * 

3 F0RMAT(5X,10F8.2) * * t 
GO TO 100 

50 CALL EXIT 
END . 
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PROGRAM 9 * 

C***** PROGRAM 9 - APPLICATIONS OF MATRIX METHODS 
DIMENSION AA(20,20),IHEAD(40),A(20,20) 
NR = 2 
NP = 5 

5 READ(NR,10,END=60) IHEAD , / 

10' FORMAT (40A2) * 
^WRITEtNP^O) IHEAD # 

20 F0RMAT(im,40A2//) * 
WRITE(NP,25) • 

25 FORMAT ( §X, ' INPUT DATA',//) 
READ(NR,30) NROW,NCOL 
- 30 FORMAT (215) 

• WRITE(NP,21) NROW, NCOL i 

21 F0RMAT(5X,!ROWS = \ 1 3, 5X," 4 COLUMNS = ',13,//) 
NCOL = NCOL + 1 

DO 40 I ^ 1 ,NROW 

READ(NR,35) (AA(I,J),J. = 1,NC0L) - - 
' 35 FORMAT (10F5.0)\ 

WRITE(NP,45) (AA(I,J)J = l.NCOL) ' 

40 CONTINUE K 
CALL ECHEL(AA,A,NROW,NCOL) 

WRITE (NP, 41) 

41 FORMAT (//5X,' ROW ECHELON FORM'//) 
DO 50 I = 1',NR0W 

WRITE(NP,45) (A(I,J),J»1,AC0L) , 
45 F0RMAT(5X,10F10.4) 
50 CONTINUE 

60 TO 5 
60 CALL EXIT * 

END 
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SUBROUTINE ECHEL 

x*LIST ALL * 
- *ONE WORD INTEGERS . * 
i C***** RETURNS ROW ECHELON FORM OF A MATRIX. 

SUBROUTINE ECHEL(A,AK,NROW,NCOL) 

DIMENSI0N-A(20,20) t AK(2O,20) 

DO 10*1 = 1,NR0W 

DO 10 J = 1 ,NCOL - \ 

10 J(IJ) = A(I,J) 

' " NC = NCOL - 1 X 
C***** RETURN IF COEFFICIENT MATRIX ROWS REMAINING' ARE ALL ZEROS 
1 5 



DO 20 I = K,NROW 
DO 20 J = 1 ,NC * 
IF(AK(I,J))30,20,30 
20 CONTINUE 
RETURN 

C***** FIND THE FIRST NONZERO COLUMN ENTRY 
30 DO 40 J = K,NC 

DO 40 I = K,NROW 
IFtAK(I,J))50,40,50 
40 CONTINUE 
50 IC = J 
IR = I 
V DO 55 J = 1 ,NCOL 
X C = AK(IR,J) 

^AK(IR,J) = AK(K,J) 
55 1\k(K,J) =,C 
X<£J<(K,IC') 
DO 60* N 4 = IC.NCOL ' 
AK(K,J) = AK(MJ/X 
60 CONTINUE* ' 

Da 70 I f 1 ,NROW 



W 



.0OOl*ABS(W))66,67,67 



IF(I-K) 65,70,66 
D = AK(I,IC) 
DO 67 J = IC,NCOL 
W= Afc(K,J) * D 
AK(I ,J) = AK(I,J)* 
LF(ABS(AK(I;J)) - 
AK(I,J) = 0.0 
CONTINUE 
.-, CONTINUE • 
C***** RETURN IF LAST ROW HAS BE€N PROCESSED. 
IF(K - NR0W)80,75,75 
RETURN 

K = K + 1 . <f 

GO TO 15 
END 
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SUBROUTINE SO LEO 

*LIST ALL ' v 
' ** NUMERICAL ANALYSIS SUBROUTINE SOLEQ 
SUBROUTINE SOLEQ(AA^NI ,M,BB,X,K,U) 

DIMENSION- AA(10,10),BB(10),A(10,11),X(10),ID(10),U<10, 10) 
N S NI . 

MM=Mfl 

.DO 200 I=M 
A(I,MM)=BB(I) 
DO 200 J=T,M 
200 A(I,J)=AA(I,J> 

IF(N-M)15,1,1 
' 15 IT=N+1 ' . 

N-M i ' 

DO 16.I=IT,M ' , 

DO 16 J=1,IW 1 . 
16 A(I,J)=0 

1 ' CONTINUE ' 

DO 21 1= 1,M 

21 ID(I)=I r ' * 

2 CONTINUE 
KK=K+1 4 

* ' . IS=* f ^ 
•» , IT=K 
* B=ABS(A(K,K) ) 

DO 3- I=K,N 
"DO 3 0=K.M 

IFtABS(A?U))-B)3,3,31 
31 l&l 
IT=J 

B=ABS(A(I,J)) 

3 CONTINUE 
IF(IS-K)4,4,4T \ 

41 > DO 42 J=K,MM 
.-/ X C-A(IS,J) ' ^ 

A(IS,J)=A(K,J) 

42 A(K,J)=C 

4 CONTINUE 
IF(IT-K)5,5,51 ' 

51 . IC=ID(K) 

ID(K)=ID(IT) 
' ID(IT)=IC 
DO 52>1„N ' - J 
C=A(I^T) 



A(I,IT)-A(I,K) V 
52 A(I,K)-C 
5 CONTINUE 

IF(A(K,K))71 .61J1 

61 KK=K 
K = J<-1 

DO 62 J=KK,M 

62 A(J,J)=-7 
GO TO 6 

71 IF(K-N)87 ,72,120 
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SUBROUTINE SQLEq^(ContJ 




72 A(N,MM)=A(N,MM)/A(N,N) 
" GO TO 7 

81 DO 8 J=KK,MM- 
A(K,J)=A(K,J)/A(K,K) 
DO 8 I=KK',N 
W=A(UK)*A(K,J) 
A(I,J)=A(I,J)-W 

IF(ABStA(I,J)>.0001*ABS(W))82,8,8< 

82 A(I,J)=0. 
8 CONTINUE 

v IF(K-M)22,6,120 
22 • K=KK • & 
60 TO 2 A 

6 CONTINJJE J 
DO 73 I=KK,N ' 
IF(A(1,MM))120,73,120 

73 CONTINUE- 

7 CONTINUE 
KJsK-V ' 

D09IS=l,Kl . • 

I=K-4S < 

DO 9IT=II,K~~? \ . 
DO 9J=KK,MM 

A(I,J)»A(I;J)-A(I,IT)*A(IT,J) 
CONTINUE • 

D010I=1,M . i 
DOlOJ^l ,M ' , 

if(id(jH)io, 111,10 

X(iy=A(J,MM) 
CONTINUE 

IF(K-M)10l,10,101 ' 
DO 1D2 4 IS-KK,H ' " 
ISUB=IS-K - i 

UM,ISUB)=A(J,IS) 
CONTINUE • • \" 

K-M-K 
RETURN , * 
120 K =.-1 ^ • 

• « WRITE(5,1000) ' 
. RETURN 

MOOO- FORMAT (27H EQUATIONS' ARE INCONSISTENT) 
END 
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1. CHALLENGE PROBLEMS 

* 

1.1 Introduction 

You have used calculus to solve problems that wduld 
have been either impossible or much more difficult 
'witho&t calculus . # In«the problems posed here you will 
find that calculus will -lead to the solution, but 
that' calculus jnust be applied to trigonometric functions. 

Read ^through each problem carefully.. Deci<£e 
which concepts and procedures £rom calculus are needed 
to solve eacift ^problem. Then, af£er your study of 
differentiation and integration_-appl ied to trigono- 
metric functions in'Unit 159 - Unit 161, you should ^ 
be able to find the Volutions t'o the problems. 

1.2 Out Fishing • 

Jack Jukes is out fishing on a spring . afternoon . 
First, there is no wind and his cork is perfectly 
still in the water. Later in the afternoon a wind 
comes up causing the cork to bob upr and down. 

_ From his physics course of the previous semester 
Jack knows that the vertical position of the cork 
plotted as a function of time will be a sine curve. ^ 
The graph- of the position of Jack's cork with respect 
to time 'is shown in Figure 1. . ^ 

y/ (Note: If there were no wind., the position of 
/the cork would remain stationary ,at y' * 0' as t increased. 
"Also, t = 0 is exactly 2:00:00 p.m.) With all this 
information magically* at his disposal Jack asks him- 
self, "Wha't is the position of my cork and how fast 
is its position changing at h second and at lh seconds 
after 2:00:00 p.m.? Also, at what point during the 
first 2 seconds is my cork falling fastest?" 
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h t(sec) 
2.0 



Figure 1. Graph of the position of Jack's cork. 

1.3 Gutter Gutters 

The Putter Gutter Company is planning to make gutters 
from 14- inch strips \>f galvanized steel. They are" to be 
designed as shown in Figure 2. 



i" 






X 


■ i V" 






V 



(8 - X)" 
Figure 2. Section of gutter. 



T 
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pAs illustrated, one inch on the outside edge will be 
used for the lip and one inch o*n the inside edge will be 
used for securing the gutter to a building. It is afcso 
desired that the total length of the side against the 
building be five inches and that the bottom be perpendicular 
to the side against the building as shown. The final 
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* A • x * * 

.consideration is to design the gutted such that it will 

shold the most water possible when filled' to a" depth of 

*four inches* The q'l^estions. that jieed to be answered are, 

"Where should th£ behjl between the bottom and outer side 

be?" and "What will the angle that the outer side ig to 

be bent up (angle e in* the illustration) be?" 

1.4~ AveragV Power ' * , 

Electrical power, measured in watts, is the product 
of the Impressed voltage an$ the resulting* current in 
amperes.* We have • * 4 

pfwa^ts) .= v(velts) x i(amperes). 

, When resistance is measured in ohms, we als'o have 

When the voltage for alternating^ current (AC) is 

graphed, wi.th rested! to time, the result is a sine curve. 
* » * . 

Suppose AC voltage is given by the equation, 

* * v A = 1 7*0 s in Y2 t , \ * 

where t is in microseconds^ * « 

» 

Also, suppose the resistance in a circuit is 17 t ohms 



Jhen 



Now the equation for power is 



HOsinit' 1 
1 = 4, — = 10 sin^t. 



p = (170 sin^t) (10 sin^-t) , 
where p is in watts; 
\170-*V 




t (micro sec 



-170 4 



Figured?* Graph of v» 170 sin ^ t. 
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The average power P, where p is . a periodic function 
of time t with period T, is defined^as follows: 



, 1 f T •' * 

( (average power) * „ P = £• l p dt. 



What is*tlve average power where g = 17X)0 sin f2 /12 t? 
What is the geometric interpretation of P? 4 

1.5 Pulling a 'Box " - - 

Jason Baxter and Sam Jones are having* an argument 
concerning pulling a heavy box across" a long' room. 
They have a rope tied to the box and Sam says, "We 
shoyld pull parallel to t f he floor. ".' Jason says, 
"It is better to pull at an angle.' 1 




\ \\\w\\ 

Figure 6a. .Sam's Proposal. 




\\\\\\\ 

Figure 6b.' Jacob's Proposal. 



j Figure 6. Proposed pulling angles. 

You are called in to settle the argument. You begin 
by recalling your recent physics course. First, you 
recall "coefficient of friction.'" , Friction 6r resistance 
varies f or^dif f erent surfaces. If it requires a force of 
magnitude F kg directed paralled to a horizontal surface ~* 
to pull an object of weight ^W kg steadily across the sur- 
face, then the coefficient of friction K is the ratio -of 
F to W. That is, v 
F 



K =: 



For example, if a horizontal force o£ 6 kg will move 

a box weighing 45 kg steadily across tfye floor, then the 

^coefficient of friction between the box and the floor is 

ERIC v 112 . 
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Y - 6 fr? 



,133, 





W M5 Kg 


F « 6 Kg 
















Figure 7- The coefficient of friction is K ■ -j^- = 0.133- 



. In considering theSpfoblem you assume the weight 
is concentrated in, a single point and^hen the force 
is applied to an angle 6 as suggested by Jason, 
Figure 8 illustrates the situation. „ 




Fsin 9 



S ^ | W (weight) 

Figure 8. Magnitudes of forces acting on box. 



In this case, the upward component of the applied 
"force nullifies par^ of the downward force*of the box, 
giving (W - F sin 8) as our replacement for W in deter- 
mining the coefficient of friction. Since the magnitude 
of the applied forck parallel to the floor is given 
% by F cos 6, -this is our replacement for F to determine 

6 
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tjhfe coefficient of friction. Thus, the coefficient of 
fepiction is 



Y F COS 6 

K " W - F sin 8 



Solving this equation for F we obtain 

F ± 0! Va 

K sin 6 + cos 0 - 



Now'' your^problem is to find the value of 6 
that minimizes F, and resolve the argument. 
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MODEL EXAM 



* Read each of the f oMowing^problems carefully 
and decide whether calculus is needed to solve tfiem. 

1. The angle of elevation of the top of a television ' 
tower from a point 1200 meters away is 0,3 radians. 
What is the height of the tower? 

2. Suppose that the resultant souod from guitar 
strings vibration has a voltageV(v) given by 

v = sin 2t - 2sin(t + J), What is the 'maximum ' 
voltage where 0 < t < |? 

r 

3. Suppose that a 14" (diameter) pizza is cut 
through the center in such- a way that a parti- 
cular piece forms an angle measuring 120°. 
What is the area of this 'piece of pizza? 

4. Suppose that' owners of a store want to put a 
triangular sign on top of their building which 
is^Bfer 6 meters long. They want the sign to be 
an isosceles triangle and have 6 meters of molding 
to put around the 2 sides that; "aTfe above the 
building. What should 9, the measure -of the * - 
base angles, be to get a triangle of maximum 
area? » 




Describe, a problem that involves the calculus 
of trigonometric functions.^ 

f 
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1. . ^TANGENT, METHOD APPLIED TO y - sin x- AND y = cos x 

r * 

1.1 Tangents to y = sin x • - ^ 

You may be familiar *with the so-called Tangent 
Method fo,r measuring the slope of-a line graphed in a 
coordinate plane*. In this unit, we are going to use 
the method to measure the slopes of lines that are tangent 
to the curve y = sinx at various values of x. This will 
give us numerical information about^the 'instantaneous rate 
of change of ~the function at these values of'x. - 

\ 

We recall from trigonometry that y = sin x, where. , 
y is the sine*of the angle .whose radian measure is x, 
is a periodic function witjv a period, of 2tt . It then - 
seems reasonable to consider x values such that 
0 £ x 2v . -Let us try to pick x values at .approxi- 
mately .5 unit intervals, recalling that the x values 
of 0, tt/2, it* 3tt/ 4 2, and 2 it are of special significance 
in graphing trigonometric functions. With these 
considerations w,e chose the/xva lues* 'that appear in 
Table* 1. ». 

' Notice in Figure 1 Which is the graph of y =. sin x 
that^ each small subdivision represents 0.1 unit and that 
each- large subdivision represents one .unit. It ,is 
instructive to use a common 'reference point to ocompare 
the slope of the tangent 4ine to see hGw the slope . _i 
changed as x increases'. Slide your triangle along 
your stationary ruler (procedure is explained'in 
Appendix 1) to translate from thr tangent to the 
curve to a parallel. line through the. point labeled 
* A .to compute* the value of the tangent for the *x values 



See Appendix 1 . 

y 




of 0, 0.5/ 1, tt/2, 2, 2,5, aird ir. Use the same procedure 
-to translate from the tangen^ line to a parallel line 
through the point labeled B tfo compute the value of the 
tangent for x vaiaes of ji7, 4.2,'zr/Z. 5.2, 5.7, and 
2ir. Record 1 your values to nearest 0\*05 unit. 

s 

The wdrk for x ? 0 is done for you. Line i "is 
tangent to y = sin^x-at x = 0. Line TTTTs^parallel to 
Z and contains point A. Line m intersects the^vertical 
lineup which is 1 unit to the right of point A,, atrpoint * 
P. rSince P is one unit above the horizontal line through 
A, the slope of the tangent line is 1/1 *or 1." 

Note: Although thfe tangent line Z to the' curve 
at x = 0 is drawn in foV illustrative purposes, it i§ 
advisable not to draw in other tangents. The many ' 
lines may cause confusion. * 

1.2 .Graphing Derivative of y = sin x 

After competing Table I we will plot the points 
with coordinates (x,y) where each x is an x value from 
t£e table and the y value corresponding ?o each x is 
the slope of the tangent to y = sin- x. . Plot these ' 
points on the coordinate system provided in Figure 2. 
The first point plotted will have coordinates (0,1). 

l^avnig plotted the" thirteen points using Table I, 
sketcty a smooth curve through (or very close to) all 
theseypoints. The* curve we- now hav^Tis^The graph of * 
the rate of change in y = 'sin x. ^> 



1.3 Making a Guess /fc^ : 

Let us hypothesize that tljfs curved also the . 
gragh of a trigonometric ^function of Y$u are'now 
asked, ywhat is the .trigonometric function of foriji 
y = f(x) whose graph this- curve* most closely approxi- - 
mates?" You may wish to refer to any trigonometry book * 
handy to refresh your memory about the graphs of 
trigonometric functions.' 

'■■ ■ ' ' . * 3 



\ 



My §uess is , ., 

Before proceeding further refer to page 13 . 

Hopefully, your table and graph, closely approximated 
those on page 13 and you c guessed th^ trigonometric 
function that was given. If you missed some values by * 
0;2 or more it is advisable to review the procedure in 
Appendix 1, and try to do the exercise again. 

/ c 

1.4* Taqg,ents to y = cos x 0 

* * J 
Let us use the Tangent Method again to* see if we r 

<can guess the function of x which" represents the rate 

of change of.iy = cos x.' Usei the procedure "yOu used 

for y = sin x to complete Table II oil ppge 6\ The work 

is again illus trat£d\ for x = 0. \ , ' 

Graphing the .Derivative of y =• cos x \ 

• With Table II complete we will plot the thirteen 
points with coordinafes (x, y) where each x is ah x 
vaLue^ from the .table and the y value corresponding to \ 
each" x value is the f slope of the line ta*ngent to \ 
y = cos x. Use^ the coordinate axes (Figure 3)* on £age £ 
to plot these points and sketch the curve. 



1.6 Guessing Again 

t ■ 
The question is again, "What is the trigortometrfc 

^function of the form y = £(x) whose graph this curved 



most clrosely approximates 

My guess is _J ^_ 



Refer to page 14 before proceeding. 
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" 2. NUMERICALLY CALCULATING DERIVATIVES FOR 
% * y 5 sin x AND y = cos x 

2.1 ^ Introduction 

In'Section 1 we used the Tangent Method of Appendix 
1 to approximate the instantaneous rate of change of 
y - sin x and y ■ cos x for various values of x. , Then 
for each function we plotted points P(x,y) where the 
y value was the instantaneous rate of change of the 
orij^nal function flbr the given x value. Next, we 
sketched a smoot'h curve determined by the points for 
each of the original functions. Recognizing that this 
curve represented a function in each case, we guessed 
an equation far this function. As you know, this new 
derived function is "called the derivative of the ori- 
ginal function. Thus, we are led to guess that the 
derivative of y = sin x is y ± cos x (Notation: 
dy/dx = cos«x when y = sin x) and the derivative of 
y, = cos x is y f -sin x (Rotation: dy/dx = ,sin x 
when y = cos x) . 

Now that we have formulas for the derivatives 
that may be correct, let us check further using 
numerical calculations. For each function, let 
us numerically calculate the average rate of change 
over various intervals with a fixed x value (call 
it x^) as one end point and numerically approximate 
the value of the; derivative at x^^ by letting the 
lengths of the intervals approach zero. 

2.2* The Procedure Explained 

Again, consider y ■ sin x, where y is the sine 
of tfhe angle whose, radian measure is x, and approximate 
the, value of the derivative .at x = x, . Our intervals 



See Appendix 2. 
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along the x-axis will have x^ as one end point and x^ + Ax 

as the oth'er end* point. We'want to calculate the ra,tio 

of the change in y to the change iji x as we move from 

? 1 to P 2 where ? 1 has coordinates (x , sin % and ? 2 

has coordinates (x 2 ,y 2 ) where x 2 = + Ax and 

y 2 = sinCxj* + Ax)..' Now, 1 the change rn-*x is Ax and the 

change in y is Ay = sin^ + Ax) - sin^x. 

'We will f^rst choose positive values for Ax (Refer 

'to -figure 5) and pick them so tjiat each successive 

choice is closer to zero than the preceding one. & We 

will then choose negative values for Ax (refer to 

figure 7) again packing them so that each successive 

choice is closer Jo zero than the preceding one. By 

observing Figurts 6 and 8 we see that* in either case 
* 

^y/Ax should approach the value of the derivative at 
the point, wheje x = x^ . In this way we will get a 
decimal approximation of the value of the derivative 
of y = sjn x at x = x^ We will then find the value 
of cos x^ and if the approximation is close to the 
value of cos x^^ we will have further reason to believe 
that out- formula is correct. * # 

» * 
, 2.3 Applying the Procedure \ 

We will now use the procedure just discussed to« 
approximate the, value of the derivative *of y ■ sin x 
at x = 0.5. -Figures 5 through S illustrate the material 
just discussed. 

\^ We will give 0.8776 as our approximation since we 
;get this value as we* approach for both the left and 
right. We find on our scientific .calculator that 
correct t*o four decimal plages* cos .5 = 0.8776. Thus 
we have further reason to believed that dy/dx = cos x 
when y = sin x. 

Next, we consider y ■ cos x v *and estimate the 
derivatives at x = tt/3. •» In recording values in the 
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Figure 5- Ax is positive. 




Figure 6. Positive values of 
Ax approach 'zero. 




. x «.5+Ax 



Figure 7. Ax is negative. 




x 2 ».5+Ax k » 

Figure 8. Neqative values of 
Ax approach zero. 



Our table to approximate the value of the 
'derivative -follows . 
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TABLE III 

Approximating the derivative of y * s i n x at x * 0.5 



* 


Ax * 


Ay = sin(0 


.5 + Ax) - 


„ s in 0 . sS 


Ay/Ax, 


* 


0.1 


*0 


. 08521' 




0.8521 




ni 

^.01 


* 






.87 51 


- * 


'•001 " 




.0008773 




~ -v&773 




. 0001 




.0000877$ 




.8776 




-.1 * 




.09001, 




.9001 




-.01 




.008800 




.8800 




-.001 




.0008778 




.8778 




-.0001 




.00008776 * 




.5776 



y column we list»four decimal places plus the number of 
decimal places in Ax. In the Ay/Ax column we will record 
four decimal places. 

< Exerc Ises 

* 1. You should complete Table IV and the. Sentence following 
the table. 
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TABLE IV 

Approximating the derivative of y ■ cos x at x * tt/3 



f 



Ax* 


Ay = cos(ff/3 + Ax) - cos (tt/3) 


Ay/Ax 


• 0.1 • < 


-O.8896 


-O.8896 


.01 . 


-cu 008685 


-0.8685 


.001 


t 




.0001 






-0.1 


0.08396 


-0.8396 


-0.01 


.008635 


-0.8635 


-0.001 > 






'-0*000,1 










4 - 



With our "approximation of 



we find that -sin(Tr/3) 
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2. Use this numerical method to estimate the derivative of 
y ■ cos x and x * 2. ' Compare the result with the value 
of cos 2. . . 



Ax 


Ay = sin(x + Ax) - sin x 


Ay/Ax 


0.1 






.01 






.001 - 






.0001 






-0.1 






- .01 






4f -ooi 
- .0001 


• 

9 





3. Again, use numerical method to estimate the yJeVivat 1 ve 

of y n cos x at x = tt/6. Compare the result with tl# value 
of - sin tt/6. , v 



. Ax 


Ay = cos (tt/6 + Ax)' - cos 


Ay/Ax 


0.1 






. .01 






.Q01 






.0001 






-0.1 

- .01 

- .001 

- .0001 







2.4 Using Degree Measure t 

You may check your results to, Exercises 1, 2, and 3 
with those on pag£ 12. You now may be reasonably 

.11 

V 

4 
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convinced «that we have chosen ,the correct formulas. 
In all these calculations and in the work in Unit 
we were evaluating- the sine or cosine of an angle • 
given its radian 'measured- 
Let us consider y = ,sin x where we are evaluating* 

the sine of the angle whose degree measure- is x. 
/ 

Now, we will use the same procedure to approximate 4 
the derivative of y = sin.x at x = 3s! This time we 
•''will take the sine of the angle whose degree measure 
is x. The results appear in Table *V. * 



V 



TABLE V 

Trying to approximate^the *deri vative of 
y = sin x at x = 3S usirtg. degree m^sure 



Ax 


Ay = sin(35° + Ax) < 


• sin 35° 


Ay/ Ax 


5 


0.0692 




0.0138 


3 


.0420 




.0140 


1 . 


.0142 




.0142 


.1 


.00143 




* .0143 


.01 


.000143 




.0143 


-5 






0.0147. 


-3 


-[ .0436 J 




.0145" 


-1 






.0144 


-0.1 


Vyooi43 




.0143 


-.01 


fl .000143 




.0143 



Thus,, our approximation to the%erivative of y = siiKX*' 
and x » 35° is 0.0143. The value of cos 35° = 0.8912". 

This is not x at all. close to what we may have expected 
from our \/ork in Section 1* In trying to salvage some- 
thirfjg we recall that all the previous work used radian 
measure^ Maybe we should have stayed with radian measure, 
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In fact, in Unit 166 we prove that dy/dx = cos x when 
•y f sin x and x is the radian measure of tfie angle. 
The problem raided by measuring the angte in degrees . 
has yet to be resolved. - "j^ 



4 



j 

\ 



\ 
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TABLE IA 




Figure 2a. Graphj^f the derivative of y « sinx . 

* . Your values in Table P should be close 
to those listed^here in Table IA and your* 
graph in Figure -2„sbould 'be .similar to 
Figure 2a above. ■ v 

The curve graphed in Figure 2a looks 
■like the graph of y = cos x . Was that your 
guess?* j 



X 


•Slope of 
tang'ent to 
y = .sin x 


' 0 p 




.5 


,9 


' 1 


.55 


tt/2 


0 * 


2 


-.4 - 


2.5 


-.8 , 


o 

7T 


-1 


3.7 


' ^-.85 


4.2 


-.5 


.3tt/2 


0' 


5.2 


.45 


5.7 


.8 


2ir 


1 




3.7 *f.2 3tt 5.2 5.7 2tt 



Figure 3a. Graph of derivative of y » cos x° 



You should ckeck your entries in 
Table II with those listed inSTable IIA ' 
at the right. Your entries should be 
close to these.. Yoilr graph in Figure 3 
should\look like Figure 3a above. 

H\e curve graphed in Figure 3a looks 
like theVgraph of'y = -sinx. Did you 
guess thi9tf 

* 



TABLE IIA 



X 


Slope of 
tangent to 
y,'^ cosoc 


0 


0 a 


' .5 


-:s 


1 


-.85 


ir/2 


-1 


2 


-.9 ' 


2.! 


• 

. -.6 


7T 


* 0 


3.7 


- .? 


4.2 


.9 


3tt/2 


1 


. 5.2 


.9 ' . 


5.7 


,55 


2ir 


0 
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MOffEL -EXAM 



*• Based upon the results of your graphical work and 
numerical calculations, complete the following 



statements : 



a- ^ (sin x) = 
b. "g^- (cos x) /= 



2. Complete the following statements. ^In determining 

the derivative of y = sin x graphically, the y valifc> 

; . wa§ the sine of the angle whose * % ^Jli 

measure was x. - 

3 I * Is iz * im Portant to use a particular unit of measure 
for angles to get the results that you listed in 
answering problem 1? 

4. From the graph below, determine geometrically the 
, ^ value (to nearest tenth) of the derivative of 
y » cos-px at x 1 = .8 and at x ? = 2 (radians). . 

9 - 

,y A 




iai 



a£ x = t).8, the- value of"^ cos x as approxi- 
mately ~ * . * 
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b. At x = 2, the value of ^ cos x is approxi 
mately * 

Complete the headings, then, use a scientific 
calculator to> complete the following table from 
which you will aprgoximate the value of the 
derivative of y = sin x at x = 0.4 radians. 



Ax 


f — ■ 

ay 








.1 

""".01 

!qoi 

.0001 






-.1 

* -.01^ 
-.00Io 
-.0001 


9 





The value of the derivative y = sin x at x * .4 
is approximately ,> . 
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OF y 88 sinx AND y « cos x AND APPLYING THE RESULTS 
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1. PROVING. THE FORMULA FOR TUE DERIVATIVE OF y = sin x 

l'.l Applying the Definition of Derivative 

^ ^Let us now try to prove our conjecture tha-t-^ 
j£ sin x = cos x where radian measure of angles is 
used. We begin by applying the definition of derivative 
where the Ax used in previous work is replaced by h. 
We must show that for arbitrary x, 

i A m s i n (x + h) - sin x 

lim * H = cos x. 

h+o h * 

Now, using sin(a + b^) = sin a cos b +^P>s a sin b, 

lim s i n (* * h) : sin x 
h ? o h 

becomes 

1* sin x cos h + cos x sin h - sin x 
h-fo h * 

^Upon collecting the sin x terms and writing as a sum, 
we have 

lim pin x(cos h - 1) + cos x(sin h) "1 

Applying laws for limits and keeping in mind that 
x is fixed, we rewrite the previous expression as' 

sin x lim (cos £ : ♦ *os x lim ilgJL . 

In order for our conjecture to hold ujf, the first limit 
must be zero and the second limit must be one. 

1.2 Some Numerical Calculations 

Before undertaking attempts at $proof., let us use 
a hand -calculator to compute values of these expressions 

N 

«w - 1 
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for {* values, close to zero Complete the following 
ta^'e^ remembering wa are using radian measure, for 
angles, * * / . %. } 

TABLE I . 

Considering ^J— 1 and cos £ - 1 for small values of h 



J, h 


sin h 


siiv n 


cos ft 


cos n - 1 




h 


h 


,2 










.1 






• 




.05 ' 




c 






*01 












« 




ft 




* 

-.1 






m 




-•05 










-.01 




t 







The values just recorded should lead us to believe 



that 



and 
v 



Urn cos t : 1 * o' 



as we had hoped. 



1.3 Proof that lim ^^r-^^ 1 
h-o , h . 



4* 



We will now attempt to prove that 
lin SiaJl - i . , *. 
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It is clear that for this acute angle we have 
Area AOPQ < Area Sector OPR < Are'a LOKSr." 
X 





0 Q(cosh,0) R(1,0) / 

Figure 1. Considering iROS with radian measure h. 



Since OP = 1, we' find that Q has coordinates 
(cos h^O) and P has coordinates (cos h, sin h) directly 
from the definitions of sin h and cos h. Thus 
Area AOPQ * (%) (base) (height) = (%)cos h sin h. 

Next, from " 

.]•' Area of Sector m Rdn measure of angle of sector 
Area of Circle ~ ~ " 2ir — 



we have 



Area -0 P R h . 
I ,2 r 2 



Now, in order to find Area AOSR we need to find RS. 
Since AOPQ =.A0SR,« we have . 



So 



RS m 
sin h ~ cos^ h * 



Area AOSR = &) (.base) (heigtit) 

- CD (fip4> . • 
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Substituting in our inequality involving these two 
triangles and the sector, we have 

(Wcos h sin h < \ < • 

• • 2 
Multiplying by gin which is positive since h is 

positive, we obtain 



cos h < 



h . 1 



sin 



TT < cos h " 



Next, we use the'fact that where a,b,'c,d are all 
positive, | < j if £ > ^, Using this on each half 
of the compound inequality *just obtained we get 



1 _ sin h ^ , 
cos - n" > li > cos h 



lim — i-r- =|=1 
h+0 cos h 1 



lim cos h = 1 
h+o- 



Now, 



and 



30 

li, linh 
hU h 

* 

must be 1 since si P h is sandwiched between — 

h cos h 

and cos h. 



Now, we need to show that 

l im S3*Ji = i. 
h*o h 

To do this let h = -t-> where t > 0. With this* sub- 
stitution we^have 

l im '^2JL = i im sin (-t) m 
Recalling that sin (-t) ■ -sin t, we get 
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li» iMlll . lim isin_t m Um sin_t . lt 

.Thus 



h-o h 

and combining this with the proof for h-*o , we have 
proven that 

i a tti sin h _ - \ 
lim — r — ■ 1, ♦ 

h-o n 

wherl radian measure of angles is used. ** 

1.4 Proof that lim cos {? I 1 = 0 

h*o n - ' 

^ We will now. try to show that . - ' 

l im cos h - 1 a 

using our last result. Knowing that -sin 2 h = cos 2 fr : 1 
we will multiply to "obtain an equivalent fraction with 
cos 2 h - 1 as numerator.' ^ USt 

lim cos *h ' 1 (cos h + 1) 
h-o • - h * (cos h + 1) 

= iim cos 2 h - 1 . .. . -sin 2 h 
J™ h(cos h - 1) " ft ,h(cos h - 1) 



: lim s * n h . -sin h 
Vo °~ h(cos h - 1) 



Urn . lim " sin h 



IT" ' hHt0 h(cos h - 1) . 



= o • Um ' si g h \ = 0 
h*o cos h - 1 

1..5 Conclusion 

Recall that in Section 1.1* we found that 
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sinOc + h) - sin x 
h+o . R . * 

= - sin x . lim (cos j| - U ♦ cos x . i im sinh t 

With the resul,ts of Sections 1.5 and 1.4 the previous 
expression becomes sin. i. (0) + cos x (1) or cos x 
which completes the proof that g^- sin x = cos x. 



2. DERIVATIVE OF y = cos x 
c 

> % • 

2.1 .Introduction 

Now, we know that \|~ sin x = cos x where x is 
/ any real number and we take the sine of the angle 
whose radian measure is x. By the Chain Rule 

d d . , ydu 

ax sin u ■ air sin u '/3x 

where u is a dif f erentiable function of x. We use 
the Chain Rule to obtain derivatives for the other 
trigonometric functions. 

2.2 Proof 

> * * 

We Will now prove that we jiere correct in our 
guessabout the derivative of y = cos r ' We use the 
identities j 

cos x = sin - X) 

and* 

cos (x - y) = cos x cos y + sin x sin y. 
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d d . r ir v 

^ cos x « ^ sin ( 7 - x) 



« (cos j • cos x t sin ^ • sin y) 
'= (0 • cos x ■+ 1 • sin x) (-1) = -sin x. 



3. WHEN, DEGREE MEASURE IS USED 

3.1 t Geometric Consideration 

Let us now attempt to resolve the problem that 
arose in Unit 159 where we computed the derivative of 

'y = sin x (at x = 35°. , We hoped to get cos 35°. Let 
us use the notation sin x° if we-* are taking the sine 
of the angle whose degree is x and the notation sin x 
or sin x (radians) if we are taking the sine of the 
angle whose radian measure is x. We consider the 
following graphs. Observe in Figure 2 that -0.801 
the value of cos ,2.5 does not seem to disagree with 
what the slope of the tangent to y = sin x at x = 2.5 
looks to be. Now, in Figure 3, does the slope of the 
tangent* to sin x at x = 2.5 appear to be 0.9996? 

"It should be if ' ^ sin x\ = cos x°, Now we see 
geometrically that we should not expect the result 
we obtained when we were using values of trigonometric 
functions whose radian measure was x. 
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y « slnx 




Figure 2. Radian measure. 



.999 



y ■ slnx 



027 



- IT 



.999 



y » cosx 



997 



Tangent at x ■ 2.5 
distinguishable 
from curve. 



not 



.082 



2.S 



—i — 

3tt 



Figure 3. Degree measure. 



\ 



3.2 Obtaining a Formula When the Angle is Measured 
in Degrees 

Let us obtain a formula for 

AKsin x°) 
^Tx * 



Where u is a dif f erentiable -function of x, the 'Chain 
Rule gives 



'■ERIC 
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- d[sin u jradians)] = cos u (radians) < fa 
Now, sin x° = sin x) . 



Thus, 



d(sin x°) _ _d_ 
3x " dx 



[sin \^ x)] 



by Chain Rule. 

3.3 Reconsiderations 

Returning to ^tlr geometric consideration of 
y = sin x°, where' x = 2.£ we have 

dCsin 2.5°) _ . v - 

*■ 

^-"(0. 9990) C 3 ^ 16 ) = 0.0174 . 

This certainly looks lojce* a much more" believable value 
for the slope of the tangent line at x « 2. 5 "which is 
sketched in Figure, 3. ** . * * 

We now can express d ( sin ^ 35 °> in terms of c 0 ^35°. 
If ^ 
■ d ( sin 35 °1 = c8s 35° > -5i 



0.819? (%^*% 0.0143 



Our approximation of o/oi43 that we obtained in Unit 
159 on page 12 now l'ool£s good. 4 ' . * ; 

' K • - 

fl?rom^ou on ufcen ve ^differentiate trigonometric 
funo-Uone we. will always' uee rai dan "measure. 
/ 
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4. PRACTICE PROBLEMS INVOLVING sin u AND cos .u 

4,1 Finding Derivatives 

'* Let us combine our new knowledge with ^previous 
techniques for finding derivatives to work some problems! 
By the Chain* Rule, where u is a dif f erentiable function 
of x?, we have * # 

d d du du ' 

37 s u = 3u sin ° u ■ 37 - cos u * 37 



and 

37 cos u - atr cos u ' 37 = - sin u 37 



d ^^V., d du . du . 



These formulas are used in the following examples, 

^ ^ Example 1: ^ sin 3x = 3 cos 3x 

'Example 2: ^ sic (x 2 + 1) = 2x cos (x 2 + 1) 

Exampld^: ^ c^s (2x - 3) = -2 sin (2x - 3) 



Exercises 

For each of thefoUowlng find . ^ 

1. y»sln2x ? X 

2, y = cos 2x • r 



3- y ■ cos (y 2 - x) 

4. y - sin (x/3) 

5, y » cos x 



In 'the following examples, we use the .formulas for 
taking derivatives where a sum,, p'roduct or quotient is 
also involved* 'These formulas are given in Appendix 3* 
- if you need to review them. 

Example 1: y ■ sin 2x + cos x 
dv 



* ax a 2 cos 2x * sin x 

is 
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Exampijg^2: X = 2^ sin cos x . , . 

dy ^ 

* 3x " 2 sin x (" sin x ) + "cos x (cos x) 



3x"* = 2 [(- si n 2 x ♦ cos 2 x)]. 
' ' Example 3: y? = \ in 2x 

r < ♦ 1 + COS X • ' » 

[ djr ^ (1 + cos x) 2 cos 2x - (sin 2x)(-sin x) 

(1 + cos x) 2 

* a 2(1 + cos x) cos 2x + sin 2x 'sin x 
x (1 + cos x) 2 

Exercises For each of the following find g£ . 



6. 


y 




sin x + cos x 










7. 


* 




x 2 cos 2x 


8. 


y 




cos 3 (2x) 




y 




sin x + x 


10. 


y 




cos ,2x - 2 cos x 


11 . 


y 




sin 2 x cos 2 x 



4.2 Finding Antiderivatives 

We recall from our work w^th^antiderivatives that 
/f(u). « f!\0 ♦ c 

where 

Thus » . . * 

* /cos u % du = sin u + c 

since » 

jj sm u = cos u, j, * ■ 

and 

* r • „ > 

J -sinvU du 3 cos u + c 



since 

d 



-T— cos u = - sin u. — 11 

au _ • > . ^ 
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The following are ^examples using these formulas' 
plus theTorinulas 

• a/f (u) du = /a f(u) du 
and • ' * 

~J /f(u) + g(u)du = /f(u) du+ ig(u) du. 

Example »1: / sin xdx = - /-s*in x dx = -cos x + c. 

Example 2: /cos 2x dx = Jjcos 2x (2dx) 

■ h sin 2x + c. 

Example 3.: /(2 - sinJjx) dx = {Z dx +/-siti J*x dx 

= /2 dx + 2/-sin(J5x) (Jjdx) 
= 2x + 2 cos Jjx + c. 



Exercises 




/ 


t 


12. 


'/cos (-2x) dx 








13. 


/sin (£) dx 










/2 cds x dx 






* \ 


15. 


• /2 sin 2x dx 






> 


16. 


/3 cos (j) dx • 








17. 


/is in (£) dx 








18. 


/(3 - sin x) dx ^ 

> 




4 




19. 


/(cos x + sin" 2x^x 








20. 


/sin 2 x dx 








•21. 


/cos 2 2x dx 






> 


Hint 


for 20 and 21 i e sin 2 


x - id 


- cos 2x)i 






* cos' 

i 


x - id 


+ cos 2x) . - 
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5. CHALLENGE PROBLEMS REVISITED 

5.1 , Introduction to '^S olutions- 

' — - . 

Now we are armed with new knowledge and' skills in 

calculus where sin u and cos t u are involved. With this 

145 . ' . ' 
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additional- ammunition let us return to battl e w ith the 
challenge^problems hoping for a successful ^outcome . 

/ • N 

If -you have trouble getting started on a problem 
or hit a snag -refer to the discussion for helpand then 
tVy to continue on your own. After completing your 
solution to a challenge problem compare] your work with 
the solution given. 

5.2 Out Fishing Again 

Our first step is to use the information from the 
graph in Figure 1 of Unit 158 to d&^mine the equation 
for the position of Jack's cork. * 

From our knowledge of trigonometry it is clear 'that 
we have the graph* of an equation of the form y = A sin Bx, 
where A is th£ amplitude and ^ is the period. Thus, 
from 'observing Figure 1 of Unit 158, we see that A = 4 
and y » 1.6 or B * j. Substituting, y = 4 sin ~ t 
is the equation for the position of the cork. The 
^^sition of the cork at t = 0.5 sec. is 

y s 4 sin (1.5)] ='3.969 cm. 
# At t ■ 1.5 Sec. we have 

y = 4 sin [^L (1.5)] = -1.531 cm. 

Now, the speed (or instantaneous rate of change of 
position with respect to -change in time) is Thus, 
we have , * 

■ ~ at 4 I— J cos ^ t = 5tt cos ^- t. 

So, the speed at \ = 0^.5 sec. is 

g£ , = 5tt cos flrL (.5)1 = -6.011 cm/sec. 
\it=.5) J' 

:he 



and the speed a,t t = 1.5 sec is 
« 



13 
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14.512 cm/sec, 



The last question to be answered is "At what point 
during the 2 second* interval is the cork falling at 'the 
fastest rate and what is the rate?" We wish to mitii- 

mize ,at s - 5lT cos TT t on the interval (°> 2 )- The 
minimum value of cos y t on (0,2) is -1 which is 

attained at-J^L t = tt or t = 0.8 sec. This will 

• • • dy VVc 

minimize ^^-^Thus 



(mill) 



5ir (-1) = 



$ 

•15.708 cm/sec. 



5.3 Maximum Putter Gutter 



Referring to Figure 2 of Unit f 58, we see that 
if we maximize the cross-section we will maximize the 
capacity of the gutter. We also observe that 6 can be 



chosen so that 6 < 



T 




From Figure 2 of Ujnit 158 we have esc 6 
x = 4 esc 8 and cot 8 = J or y = 4 cot 8 . 
A(Area of cross section) 

» Are^> o^f Rectangle -< Area of % Triangle 

• 4[(8 - x) + y] - Js(4) y 

= 4[(8 - 4 esc 6)'+ 4 cot 6] - *s£4 cot 6) • A 

(substituting 
^ for x and y) 

- 8(4' - 2 esc 6 + 2 cot 8 - cot 6) N 
■ 8(4 - 2 esc 6 + cot 6) . 

An We don't 

know how to find the derivatives of esc 6 and cot 8 
at, this time. What we can do is express esc 8 and 
Cot 8 in terms of sin 8 and cos 8. Doing this we have 



da 

The area will be maximum when jg- - 0. 
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q , A 2 . cose. 

8 (4 -iTJTe + iTire ) 



dA v s g r f sin 8 (0) - 2 cos e ] + f sin 8 (-sin 8) - cos 8(cos 8)| 

L I -sin 2 8 , i I sin 2 8 • J 

a Q pcos 8 + (~l)(sin 2 8 + cos 2 8) 1 ■ 
Lain 2 8 sin 2 8 „-J 

. L <s in 2 6 -I ^ 

since sin 2 8 + cos 2 8 = 1. For ^ = (► we must have the 
expression within brackets equal to zero, but thismeans 4 
that the numerator must be zero. Thus, 2 cos 8-1-0 
or cos 9 « 3s. Since 8 < j<we have 8 - ~ or 60 p , We 
also observe that for 8 such that 0 < 8 < ~ we have 
cos» 8 > »V-which implies 2 cos 8 - 1 > 0. Mow looking 
at our expression we note that when 2 cos 8 - 1 > 0 
we have 44 >4k In a similar way we can conclude- that 
j£ < 0 when y < 8 < |. " 

We have now verified that 8 = ~ is in fact the 
value of 8 for which the 8 area is maximum. 

Since x = 4 esc 8, we have 



and 



4 esc j = 4.619 



8 - x = 3.381. 



Thus, the bend should be at approximately 3.381 inches 
from the bend for the right angle and the metal should 
be bent up ~ radians or 60° for a gutter with maximum 
capacity. ' 

5.4 Average Power Computed 
^ We have (Average power) 

P = 4 / p dt, 

0 . 

, 148 * < 
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where 1 

p » 1700 sin 2 ^ t . 

Observing- Figure 5 of Unit 158, we see that the graph 
of p repeats itself every 12 micro -seconds/ Thus T = 12. 
' Substituting, gives . < 

P - A ^170'0Wina^ tdt = ^ \ %dt 

This problem is like Exercise ^O.'in section 4?3. We* 
will replace ' * 



br 



• 2 ^ 

sin 2 ^ t 



Ji(l - cos i.t) 



P ■ if Ad " cos I t) (ft 

0 ° 

= 1700 r r 12 .. , 12 
24 



r[/ o d, : _/ o c. 8 f,'d.] 

1700 [% 12 a 12* i 

;:#[{, d vf / 0 l?r : l tdt ] 

■ [|12 - 0J - | (sin 2* - sin 0)] 



1700 - 
24 



1700 n9 . * 1700 

( 12 ) " —J— = 850 watts. 



Now, consider the geometric interpretation, 
^ T7 



12 * ' , . 



represents the $rea, (A) under one hump of the curye 
\ shown by //// in Figure 4. # We k^ovt*— A = 850, so 



i * 
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A = 12(8j£p. Thus the rectangle shown by \\\\ with 
height equal to the average power P has the same area 
as the area* under the hump. 



1700 



w (wajts) > 




t (mi cro sec) 



6 12 . 18 2k 

Figure k. A geometric interpretation of average power 

- j. 



r 



5.5 Pulling a Box Correctly 



dF 



* Finding ^ using the^ quotient rule gives 



dF = -KW (K cos 6 - sin B) 
37 / (K sin 6 + cos 9) 2 



if 



or 



= 0, we must have 

K cos 8 - sin 6=0 



in ft 



17 
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Thus, the first derivative is zero for 9 such that 
tan- k we can also show the first derivative changesi 
from negative to positive at this value of'e which, 
verifies that this is tjie angl o ■ w ke-re the force is 
minimum. N0w, K is always some positive number, (go 
Jason is correct and the angle depends on the coefficient 
of friction. 




Figure 5'. Coefficient of friction. 



r 



If one object has a higher coefficient of friction 
th^n another object on a particular surface then the 
.angle at which the force should be applied is greater. 
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6. MODEL EXAM 

1. What technique is used to show lim sin h = 1? 

h+o n r 

2. : What limit other than lim sin h - 1 was used to 

h+0 h 

d N 
prove (sin x) = cos x? v 

Complete the statements in 3 and 4. 

3. The value of ^ sin x is b ,when 

x = 17°. 

*d 

4. The value of -r- cos x is ,when 

x = 39° , 



Mn problems 5 thru 10, find ^ 



5. 


' y 




sin (x 2 - 


3) 


6. 


y 




2 sin x + 


Cos 3x 


7. 


y 




sin 2 x 




8, 


y 




si^vx cos 


X 


9. 


y 




cos(3x 2 - 


x) 


10. 


v y 




cos 2 x sin 


X 



Find the antiderivatives in Problems 11 

11 | sin x dx 

' ♦ 

12. | cos (3x) dx ' 

13. | (5 - 2sin x) dx 



14 • Giv f n F = - K 'sin e" cos 9 • find § 



> 1 
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7. ANSWERS TO EXERCISES 



/ 1. y * sin 2x 2 

' * g£ e 4x cos 2x 2 

2. y = cos 2x 

37 = " 2 sin 2x 

3. , y = cos (x 2 - x) 

* ^ ;}£ " -(2x - 1) sin (x 2 -x) 

4. y = sin ~ 
dy _ 1 „ x 

_ -» .c-i cos i 

5. y = cos x° = cos (^ x) 

TT 

180 

6. y = sin x + cos x 
g£ = cos x sin x 

7. y = x 2 cos 2x • 
g£ - -2x 2 sin 2x + 2x cos 2x 
dv 

g£ - -2x (x sin 2x - cos 2x) 

8. y = cos 3 (2x) 
-.6 cos 2 (2x) sin (2x) 
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9, y ^ sin x + x 
• j£ e cos X + 1 

10, y « cos 2x - 2 .cos x 
dv 

g£ = -2 sin 2x + 2 fcin x 
j£ * 2(sirf x - sin 2x) # 

- 153 



. y sin 2 x ct)s 2 x 



» -2 sin 2 x cos x sin x + 2 cos 2 x sin s cos x 
= 2 (cos 3 x sin x - sin 3 x cos x) 
= 2 cos,x sin x (cos 2 x - sin 2 x) 

-hi -2 cos (-2x) dx = \ sin (-2x) + c 

; 

-2j-h sin (|) dx = -2 cos (|) + c 

2/'cos xdx = 2sinx + c i 

-J-2 sin2x dx = -cos 2x + v c y 

9/J • 3 cos (J) dx = 9 sin (|) + 

-/-% sih (J) dx = -cos (J) + c N 
J(3 - sin x) dx = /^dx + /-sin dx = 3x + cos x + c 
/(cos x + sin 2x) d.x = /cos x, dx -Js/-2 sin 2x dx 

* = sin x - h cos 2x + c 
/sin 2 x dx = is/ Cl - cos 2x) dx 

= ^[/dx - Jj/2 cos 2x dx] =| - % sin 2x + c 
/cos 1 2x dx = hf(l + cos'4x) dx 

» <' 

■ h[jdx,+ h /4 cos 4x dx] 

x . 1 - 
= T F SJLn + c 
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1. OBTAINING 'FORMULAS v 4 . • 

1.1" Introduction 

In our discussion of the solution to the .Putter 
Gutter Problem (see Section 5.3 in Unit 160^we' 
obtained the cross sectional area A of the^ gutter. 

from the equation A = 8(4 - 2 esc 6 + c^tfr*). 'To * 

dA ' 
obtain ^g- we had to Express esc 8 and £ot.6 in terms 

of sine and cos 6 since we°did not have* f ormiXLas ' 

for the derivatives of these -trigonometric functions. 

To eliminate this extra work in the future, we will * 

now derive the formulas for the derivative's Q,f the *\ 

other four trigonometric functions. We will follow' ' 

the same plan used to find and express ea'ch of 

^the other functions in terms of the sine an^\:c 

functions . 



1.2 The Derivative of tank — ^ 



To obtain the, formula for the derivative of tan x \ - 

we express* tan x /a^ and use the w/t^for quotients 

which guarantees that ^ '„ ' * 

2. . * 

d ( a v (du/dx) - u (du/dx) ■ -» ■ 

j * » ■ - 

• "where u. and v are di£f erentiable functions of x. 

* f • - " " - 

d + v - r s A nx ^ - cos x • cds* - sin x (-s in 4 x) 



COS 2 i X 



= cos 1 x + sin 2 x 1 * * , 2 . v « 



o COS, 2 X 

1/ 



sec ' x. 



1.3 The Derivative of sec x ** * 

We will^ express— s-ec-x-as cq ^ x and u^^the^quotien^t - 

rule again to obtain the .derivative of sec x. 

t 
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* sec x = 4- f 1 ) - cosx *.(y- 1 (-sinx) 
dx 3x\ ^cos X J 2 



s in- x * l ' s in x 

- sec x tan x . 



cos 2 x cos * cos x 



Exerc i ses a 

1. Expressing cot x as c ? s * , follow the procedure used' to 
obtain the derivative of tan x to obtain the derivative 
of cot x. 

1 



2., Find the derivative of esc x by replacing esc x by 



sin x 

and using the same technique as used for the derivative 
of sec x. 



l. # 4 Complete List 

* With the solutions to Exercise 1 and Exercise 2 
in the preceding section we now have obtained formulas 
. for the derivatives of all six of the" trigonometric 
/functions of ^x. ^ You should check your<pr'ocedures and 
results to .Exercises 1 and 2- with the solutions gwen 
at the end of this unit.. The formulas for the deriva- 
tives of the six trigonometric functions of x are 
listed in Appendix 4 for your convenience. " 

2. PRACTICE FINDING DERIVATIVES - 

2.1 Using Chain Rule Again - 

Having the formulas for the derivative *of t>an x, 
cot x, sec x and esc x, we wish to obtain the formulas 
for these trigonometric functions of u where u is a 
dif f erentiable ^function of x. We apply the Chain Rule 
in each case just as we^ did for sin u and cos u. 
For tan u we have • * 

d ^ d ^ du ? du 

tan u = gjj tan u . gj = sec 2 u g- , 



\ 
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The formulas for the derivatives of^the remaining 
trigonometric functions of u follow* in^the same 
manner. The formulas for all six trigonometric 
functions of u are listed in Appendix 4 for your 
reference. • * * 

2.2 Applying the Formulas 

These formulas. are applied in the following 
examp les . . q 

Example 1: * gjf tan (J) « J sec 2 J. v 
Example 2: ^ esc 2x = -2 esc 2x cot 2x. \" 
Example 3: ^ cot x 2 ' = -2x'csc x 2 . 



.U161 



Exerc j ses 
dy 

Find in each of the following. * 

3. y,« tot | . . 

4. y * esc (2x - 1) - 
5« » y ■ tan x 3 . 

6. y = sec -(-x) . v 



2 , 3 More Iryolved^Applications 

Jnwthe following examples sum, product or^ quotient 
may also be involved. These ^formulas are listed un 
Appendix. 3 for reference. Derivatives involving sin u * 
and cos u which we considered previously are included* 
in the examples and exercises that* follow. 

Example 1; y = tan 2x + ^sec x 

3x a 2 sec 2 2x + sec x tan x. 
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. Example 2: 



Example 3: 



y » tan 2x. 
g£ m x (2 sSc 2 2x) + tan 2x 
* 2x sec 2 2x + tan 2x. 

COS X 



y s 1 + tan 2x 
dy _ (1 



+ nan 2x) (-sin 



x) - 2 cox x(sec 2 2x) 
tan 2x) 2 



Exercises 



For each of the following find & . 



7* y - sin x + -tan-x-r — 

8. y - x 2 esc 2x. 

9. y » tan x sin 2x 
10- y « cot 3 (2x) . 

11. t y a sec x + x . 

12. y » cos 2x - 2 esc x, 

13. y « x 2 tan 



1*1. 



sin 2x 
- tan x 



2 .4 Puttering, Around 

. *~ We "return to the Putter GuAer Problem of Unit 158 

one more time. Let us find jg^Kow that ye have the 
formulas for the derivatives^^of esc 6 and cot 8. 



A = 8 (4 - 2^c 8 + cot 6). 
g£ = 8 £-2(-csc 0 cot 6) + ( -esc 2 8)J 

« 8 ^[2 esc e cot e - esc 2 e]„ 



We see that having the formulas* for the derivatives 
of all the trigonometric fjinct^ons' available>made 
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dA 

finding easier, and they may > be of use to us in 
future putterings. 0 ^ 
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3. MODEL EXAM v 

Differentiate the, functions in Problems 1 and 1 

by first* expressing them in terms of sin x and cos x. 
* * 

1. & Vtan x) - 
2. 



37 (esc x) = 5 

Find g£ inVeach of the following 

3. y = cot x 

k - . 

4. y ° sec x. 

6. y = sin 2x + tan x 2 , , 



8. 



problems . 



1. y = x tan x. 



cos x 
1 + tan 2x* 



i f 
It* 
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4. ANSWERS TO EXERCISES 
£ cot x - £(§§5-7). By quotient rule, 

d» ( cos x , n sin x (<-sin x) - cos x (cos; x) 
dx^sin x J : ; ~ 

sin* 4 x p 

= -sin* x - cos 2 x _ -1 (sin 2 x + cos 2 x) 
sin 2 x m sin 2 x 

- -1 (1) 2 
= -> - — *- = CSC X. 

£in 2 x 

4 CSC X = ^iTTT^- By quotient rule, 



a * ( , 1 ..' sin x CP) - 1 (cos xJ _ -cos X 

# , dx v sm x' / ~ ~~2 



cos x 

" -esc x cot X. 



s in x . sin x 

v > « 4 

Note: We essentially did these exercises in our * v- 

dA ' 
9 calculations for in 5.3 of Unit 160*- 

g-J-c.c« $)(*)'-- J c.c«f. 

4. y = esc (2x - 1) . 

dv \ • k 

3x = " csc ^ 2x " ^ cot ( 2x " x ) \ 

= -2 esc* (2x - 1) cot (2x - 1). 

5. y tan xK ' \ * * , 
% ■ (sec 2 x 3 )(3x 2 ) = 3x 2 sec 2 x 3 . 
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y = sec (-x). 



g£ ■ [sec (-x) tan (-x)3(-l) = -sec (-x) tan (rx). 
y a sin x~ + tan x. 



cos x + sec' x, 



.8. 



X CSC ^x . 



.(2x) cot (2x)(2)] +.[cscf (2x)]2x 



9. 



'10. 



11. 



* 2x esc (2x)(-x cot 2x + 1). 
y. %taif x sin. 2x. 
g£ =* tan x cos (2x)(2) + sin (2x) sec 2 x. 

y = .cot 3 (2x). 

3 cot^(2x) [ -esc 2 (2x)](2) 
\ * -6 cot 2 (2x) esc 2 (2x). 
y = sec x + x. * * ' 4 



l2* y = cos 2x - 2 cse x. 

-sin (2x)(2) -,2(-cse x cot x) 

QX % * ^ 

'= \2 sin- (2x) + 2 esc x cot 
13. ,y = x 2 tan-' (|). 



'=-x[| sec 2 (§) + 2 ta^h (f) ] . 



2x 
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sin 2x 

1 - tan x ' , 

(1 - tan x) cos (2x)(2) - sin (2x)(-sec 2 x), , 
(1 - tan x) 2 

> 4 

2(1 : tan x) cos (2x) + sin (2x) sec 2 x 
(1 - tan x) 2 



j 



o 
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APPENDIX 1 

. THE TANGENT METHOD FOR ESTIMATING DERIVATIVES* 

Objective 1 : To be able to use a triangle and ruler to 

'estimate the slope of a line on a graph. 
Graphical Method for Finding the Slope of a Line ' 

. Many times in our. work we want to measure the slopes 
of lines plotted on graphs. We can always calculate the" 
slope o'f a line by reading the, coordinates of two points 
on the t line^ and applying the forjnula 

sloped change in vertical' units » 
v change in horizontal units 

. y 2 - yi _' Ay ' 

But there is an easier way that saves the effort ~of read- 
ing the four numbers from the gr^aph necessary to calculate* 
each slope. For this method, you will need r a straight 
edge or fuler, and a small drawing triangle. 6 

Figure l/l is a graph that shows the profile of the 
Union Pacific Railroad. The problem is to find the slope 
of the railroad betwjeen Lakeside, Utah and Wells, Nevada 
directjy from the graph using as little arithmetic as pos- 
sible. The following steps provide ar&easy method to 
measure this slope. 

Step 1 (see Figure 1.2). Place the triangle with one 
edge along the line whose slope you wish to measure. 

Step 2 (see Figure 1.3). Place the ruler Against ' the 
other side the triangle. Check t;hat the first edge of * 
the triangle is still along the line you wish to measure. , 

■ — - ! ^ i_ • 

*fcdapterf by the UMAP Project Staff from Differentiation, Second 
Edition , 1975, Project CA^Q, Education Developfoent^Cenjter , inc.. 
Newton,, Massachusetts, pp. 27-60. / ' « ♦ 
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Step 3 '(See Figure 1.4). Slide the triangle along, 
the rifler (holding the ruler firmly so it will not slip) 
until the first edge of the triangle passes through an ' 
easily read intersection of the- graph paper. (In this 
example, the tiT^angle was moved until its edge passed 
through the intersection of the 1,000-mile line with the 
5 , 000- foot'line.' Since the v tri angle was slid along the 
ruler, this edge is still parallel to the lfne whose slope 
is to be measured. The slope of the edge of the triangle 
is therefore still the same as the slope of the original 
line . ) 




Step 4 (see Figure 1.5)*. Read the; slope of the edge 
of the triangle at the point where 1 the triangle cuts the 
next major vertical line on the graph paper. Here the r 
next major vertical line is at 1,100 miles or i-s 100 mifes 
further from the first easily read intersection. The tri- 
angle edge intersects this line 1,100 feet above the first 

' ' 168 % - 



read intersection. The upward slope of the triangle, and of 
the track is therefore* 

1,100 feet • - 1 feet 
100 miles ' 11 mTTe 

or 11 feet per mile. Note that i-f you choose the hori- 
zontal distance to be 1, 10, 100 or 1,000 miles, the 
division can be easily done in your head! 




i 



Figure 1.5. 



Exercise 1 (Short Method "for Grades on the UP-SP RR) . Using the 

method above, find the slope of the track between 

Reno, NV and Tru-ckee, CA 

Greei> River and Evanston, WY 

Midsection between Omaha, NE 

atid Julesburg,. CO * 

Wells and Love'lock/NV 




0 



Improved Method for Finding Slopes 

It is possible to improve this method to avoid the 
division and the placing of the decimal point. Try the 
improved method on the same graph of the UP-SP RR you have 
been. using (Figure 1.1). As with many "how to do it" direc- 
tions, it takes much longer to describe than to do, so 
follow along and your patience will be* rewarded. If you 
have trouble following^he^irextliJiTS", . have ^our instructor 
give you a quick demonstration. 

Setting Up a Scale fqr Reading Slopes * 

Step A (Figure 1.6), Mark a standard interaction 
one major division in from the right-hand^edge of the graph 
paper. (Such an intersection has already been marked with 
a O in Figure 1.1 at the beginning of this section, and in 
which you may continue to set up a slope scale and make 
measurements.) 

Step B (Figure* 1.6) . Temporarily mark the first major 
division above the center of the slope scale with the num- 
ber ^of vertical units it separates. Here it is marked 
+1,000 feet since it represents an increase in elevation of 
this amount . 

Step C . Calculate the value of the^slope for this 
first division by taking the ratio of the vertical increase 
(1,000 feet) to the horizontal increase (100 miles) for one 
major division. % * 

slope = *L = h™* feet * ltf ft ' 
v Ax v iuu miles m^ 

or 10 feet per mile/ Write this number of the* scale in 
place of the temporary mark of 1,000 feet. 

Step D . Mark the slope values +*10, +20, +30, etc., 
opposite the main divis ions, going upward f^om zero. Place 
-10, -20, -30, etc., opposite" The main divisions going down- 
Ward from zero. (See Figure 1.7.) Write the units in which 
the slope is measured (ft/mi) at the top of the scale. You 

are now ready to use the scale to measure slope. 

• • • 

A~6 

"170 , 



To itt up a teak for direct itaditf of slop* . . . 




Figure 1.6. ' Figure 1.7. 



Measuring the Slope with the Slope Scale * 

To use"the scale just marked to measure slopes in a 
convenient and direct *way, set the triangle to, the line** 
you wish to measure and's^de it by means of the ruler ' 
until its edge passe*s through the standard intersection Q . 
Read the value of .the slope at the- point where the edge of 
the triangle crosses the slope scale line. .(See* Figure 1.8.) 

NOTE: You cannot always slide the triangle .to a ° * 
position* where its edge passes through both the standard 
* intersection^ and the slope scale in one slide along, the 
ruler. When;this happens, hold, the 'triangle, firmly in the 
^'Q position, shift tKe. ruler so it is along a different 
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edge of the triangle, and t then proceed to slide the 

triangle in a new direction. By a series of such parallel 

slides of the triangle, it is possible to position the 
triangle so that the slope may be read. 

Exercise 2 (Direct Method for Grades on the UP-SP Railroad). For 
practice, recheck the slopes you measured before and then try these: 

Section ' Slope 

1st section after Jules£erg 

♦ - .1st section before Green River 

1st section after Truckee 

! u : 



V Objective 2 : To be abla to calculate numerical values of 

the average rate of change between values of 
a function: (a) from a tablf of vaVues > and 
^ * • (b) from a graph'. 

A function in mathematics is a rule or a recipe that 
relates two quantities such as distance and time. We now 
begin a process that will lead eventually to a method fo* 
calculating the rate of change of a smoothly varying func- 
tion. f As an example we use the curve shown in Figure 1.9 
which shows the distance a small iron sphere (a ball bear- 
ing) drops as a function of the time since its release. 
The curve in Figure 1,9 was drawn by passing a smooth curve 
through a set of data points . * 

Calculating Rate of Change at a Point 

n Eventually we will want to be able to calculate the 
i*ate 'of change of a smooth function at a specific pojnt. 
For example, we may want to know how fast the distance is 
changing (that is, we want to know the speed of the ball 
bearing) 'at, say, t ■ 1.5 sec. To answer this, let us 
back off a bit aftd answer a related, but different, ques- 
tioni namely, what is the speed of the ball bearing during 
the time interval from t ■ 1,5 sec to, say, t ■ 2,5 sec? 
% T^iis question does not mean very°much without further 
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remarks. As -you can see from'Figure 1.9, the slope of 
the curve changes gradually and steadily from t = 1.5 to 
t = 2.5.* What then do we mean when we ask what the speed 
is during this interval? 

Calculating Average Speed Numerically 

To" avoid this problem, we define what is called the 
average speed* over the interval. Figure^ l.io' shows in 
detail the portipn of figure 1.9* from t = 1.5 sec to 

= 2.5 sec. The "average" speed is obtained by finding 
how far the ball bearing fell during this time interval, 
and then dividing the distance fallen by the length of the 
time interval. From the graph in Figure 1.10 we see that 

when t = J.. 5 sec, h = 36.0 ft 
when t = 2.5 sec, h = 100 ft. 

The average speed over the "inte'rval from t ■ 1.5 to t ■ 2.5 
is then: ^ 

= distance fallen = 100 ft '- 36 ft 

av time to fall this distance 2.5 sec - 1.5 sec 

. ='f^ = 64 ft/sec. ' « 

In symbols this calculation may be written: 

v - h 2 ' h l _ Ah 
av t 2 - At 

and is, of course, just the formula for the slope of the 
line from point A to* point B in Figure 1.11. Point A has 
coordinates t^ = 1.5 sec and hj - '36 ft while 'the point B 
has coordinates t 2 ~ 2.5 sec and h 2 a 100 ft. 



*lt Is vfery important to realize that "average 11 here does not 
mean what it usually means. The average speed is not found by add- 
ing together a number of speeds and then dividing by the number, of 
speeds. The average speed in the sense used here is that constant 
speed at which the ball bearing would have to move between t » 1.5° 
sec and t ■ 2.5 sec to cover the distance It actually does move. 

/ 




'A-ll 





13 



Cafrculating'Average Speed Graphically * ' 

We can also find thjls slope (average speedj using the • 
"sliding triangle 11 method described a shQrt while earlier, 
which leads quickly to reasonably accurate* results . In * 
Figure 1.12 'this method as used to find the average spewed 
^of the ball beating between t = 1.5 sec and t = 2.*5 sec.. 
The result compares very favorably with^the -computed value 
of 64 ft/sec. - - * , * , • ' 

• ' -*±* 

Objective 3 : To be-able to_estimate the^ instantaneous rate 
t I of ^change of a function 'by graph-Leal means; 

that tSj measure the slope o*f a line tangent 
v to the curve, 

ftate of Change of >a Smooth Function at a Point 

When we found the average speeji of a fairing ball- over * 
an interval of one second, the average speed was not the 
actual speed at either the beginning or the 4'nd of the 
interval. Rather, it represented that ctfnsfant speed with * 
which Ihe ball' wouW have covered the 64 feet fallen rn the 
same one second of time." Suppose now that instead of want- 
ing tj^e average , we wanted the instantaneous spee'd at tne' 
moment the clock read J. 5 seconds. The average speed over * 
an i^teryal can be measured with a tape measure and a stop 
watch: we measure the distance traveled *an3 divide ,by the 
time it took to travel that distance. But -obviousLg^an 
instantaneous speed cannot be measured or qalculaxed in the 
same way; We would need to t meafsure the distance traveled in 
a time interval of zero length., 1 

Taking the Average ovejr Intervals 
b — , # , 

We wilL find £lle instantaneous speed,, not by measuring 

a time interval of zero length, but by "sneaking up ! \ on it: 

we find the average speed oyer shorter and shorter time j 

intervals beginning at t * .1*5 sec. The average speed, 

starting at r 1.5 seconds, but measured over only p. 8. seconds 

instead ^of one second; is measured in Figure 1.13 by the 

"sliding triangle" method and comes out to 61 feet/second. 

".. " A : . . 178 i A - 14 




Figures 1.14, 1.15, and 1.16 repeat" the measurement' of > 
average speed over intervals of '0.6 S 4 ec^0.4 sec, and 0.2 
sec. The results of these measurements off the average s 
speed are shownHn Table 1. We would like to continue mak- 
ing such 'graphs in order to extend the table to include 
-even shorter. time intervals. . Rut here a practical diffi- 
culty arises: For valu^ of At smaller than about 0.2 sec, 
it is impossible to read such graphs accurately enough to 
obtain reliable results. *At the moment, the best we can do 
is 'to say that the instantaneous velocity at t = 1.5 sec is 
something near 50 ft/sec. In Appendix 2, however,\ We will 
see how this method can be used to arrive at a more accurate 
answer. 

u TABLE 1 

Average Speed of .Falling Ball Figured'over , 
Intervals Starting *at t = 1-5 Seconds 



* Time Interval 
At (sec) 


Average Speed 
v av (ft/sec) 




i:o 


64.0 




/ 0.8 


61 .0 




0.6. 


57^5 ' 




• O.k 


54.5 ' 




0.2 1 


< 51.0 » • 



: - ; - 

Exercise 3 . Use the graph of Figure 1.17 to estimate tKe instantaneous 

speed at t * 5.0 sec. Follow the procedure «descri bed above, calculat- 
ing Ah/At for a succession of^smaj leV and smaller vajues^of At,' ~* 
begiShTng at t » -2* sec.** . . _ [ . »* " 



ftateof Change of a Point from the Slope of a Tangent Line 

If- the method just used in examined in detail from a % 
graphical point of view, it*can lead to "a more accurate 1 
estimate of the instantaneous^ velocity,. Moreover, this 
more accurate ansiTfer ca*h be found without the need to draw* 
as many graphs as in Figures 1.13-1.16, or to read as many 
successive values of the average s*peed such as those in 
Table 1. *li us, look at the series of diagrams shown in 
Figure i.JL8 and notice what happens to the successive lines 
whose slopes give the average speeds. As the time interval 





e) 



;: j«;H ::::: H:HH::U::H: ::L: HH: ::H? 
» J: HH: :v£i H::: :SH: :H:: H;H :HU till: 



liililllilillll 

illlijiOiiliillllo^Siiii 



::H:bn::n:::3£H;r.:H:s::r*^ 

HH: H::::H:: 5t£ ::::: 

HIh hII I Hln ,r?s Iim IIIU HHI 



Figure 1.18. *T*e line tnat' cuts 'the curve becomes the line 
* \ that touches the curve as the two points move 

* together. 
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gets shorter and shorter, the points wfiere the line cuts* 
the curve move closer and closer together. As At ■»■ 0 (Fig- 
ure" 1 . 18a) , 'these two points unite to become one point; 



the line then touches The curve in one point only; grazing 

the curve. Such a line gi.ves the slope of the curve at 

that point and *is called a* tangent line (from a Latin'word 

meaning "to touch"). The slfcpe of the tangent line is 

defined to the 'instantaneous 'Speed at the point (that is, 
♦ 

at the instant of time) where the tangent line touches the 
curve . . * » 

This conclusion gives us a simple way to estimate the 
instantaneous speed (or rate of change) from a smooth 
graph- We 'merely draw a tangent line (which can usually 
be* done quite accurately by eye) and measure its slope. 
To measure the slope of the tangent' line accurately, either* 
draw a long tangent line and read off widely separated 
.points to compute its slope as in Figure 1.19, or* use the 
sliding' triangle method as shown in Figure 1.20. % 

Note that the value determined from Figure 1.20 is 
48 ft/sec. This, then, is our estimate of the instantaneous 
velocity at T.5 sec. * ) 

Exerc i se k . Estimate the instantaneous speed of the falling ball at 

t = 1, 2, and 3 sec by the tangent method. Use the graph of* Figure 1.17 

Exerc is,e 5 . The graph in Figure 1.21 shows the area of an open wound 
Versus time; In doing the following, use the smooth curve in the 
figure, not the dots. The healing rate is. the absolute value of the 
rate of change jjJPthe area of the^ wotmd ^ rp^asured J n cm^/day. 4 

(a) Sketch a graph of the rate of chancje of the area of the wound 
(cm /day?. 

(b) When is the healing rate'the fastest? The slowest? 

(c) Use the 'tangent method to calculate the instantaneous rate of 
^ healing at 8 days and at 1|.5 days. 
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Exercise 6 . Sse the graph shown in Figure 1.22 to find the rate o 
increase in quantity of antibody 2 days, 5> days and IS days after 
both the 1st and 2nd injections of antigen. Use the tangent method. 

Exercise 7 - Draw a pair of coordinate axes, and sketch a graphNof a 
curve that has positive slope everywhere. Do the same* for a curve\ > 
whose slope is negative everywhere, and whose slope is zero everywhere. 




O 5 10 * IS 

T>9yS after l rt injection of antljtn 




Amount of 20 
-antibody 
in to rum 



O '5, 10 , \ 15 

* Days »fttr 2** injection of came antigen ' 

Figure 1.22. The slfcw production of antibody (A)* after a first 
injection of antigen is followed by a prolonged 
decline in quantity of antibody (noj: shown/) . If, 
after this decline, a second injection of the same 
antigen r, booster shot'*) is given, the rise 1 in 
antibody level is rapid (B) . 



190 



ANSWERS TO EXERC fcSES 



Reno, NV and Truckee, CA: +35 ft/mile ■ 

Green River and Evanston: +7 ft/mile 

Midsection between Omaha, NE 

and Julesburg, CO: * +8 ft/mile* 

Wells aricbjjovelo£k, NV -6 ft/mile 

Section after Jul esburg:^ +16 f t/mi Ve (£2) 
Sect ion, before»>Green River F -5 f t/mi^ (±1 ) k .- 
Section after Truckee: *** \*jk7 Vt/mi le <±5)* 

The table sftow,s tha>t the instantaneous speed at t 
about 64 ft/sec* » . " 



= 2 sec is ^ 



4, 



At 
(sec) 


'Ah/At 
(ft/sec) 


t 




1.0 


80.0* 


f 

0.5 . . 


72 :° 


'i>.i ■ . 


fe5.6 \ „ * 


0,05 • 


*• 64.8 ' 


0.01 


' . ' 64.16 • 


0.005 


64.008 






i 


t 

t <s 
? (sec} 


, <# ec > 




e * 


* «H 


32 * 


'3 


96 



5'. (a) See ttte graph at the top .of the next page. * ' ** /' 

(b) The healing rate is. fastest at about 15 d£fys, when the wound 
Is nearly healed.. It is slowest at the outset when the 
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wound Is newly formed. 



* ?(c) At 8' days! about 0.4 cm 2 /day. 
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aboW 0.5 cm /day. 
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1st injection: "2nd day: 0 units/day * * 

5th day: 0,3 units'/day 

15th day: 2.2 units/day 

2nd injection: *2ti'd day: 15 to 20 units/day 

5th day: O^un its/day 

15th day: 0 units/day 

Note ^hat' after the 2nd injection, the amount of antibody may 
reach a higher level than it does after the first Injection.* 
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APPENDIX 2 
RATES OF CHANGE* 



Objective 1 : To be able to estimate numerically the 

average rate of change of a function given 
by a formula. 

In Appendix 1, all our information about rates of 
change came from graphs. We, now explore now to find the 
average rate of change for functions given' by a simple 
^formula.,) The function we use is one you' have probably 
seen before. If an object is dropped, the distance, h,' 
which it falls in a time t is given by 




h = 16t 2 . 

It is important to keep in mind when we 'use this formula 
that t musfbe in seconds and that h" comes out in feet. 
For example,, , to- find out how* far the object has fallen at 
t - 1.5 sec", we calculate « \ • • 

h = 16 x (1.5) 2 = 36 feet. 
Figur^ 2.1 shows a graph of this furtction. 
^ Calculation of an Average Speed of Fall 

Let* us compute the average rate of change of flp from 

t l = >fc 1 - 5 sec t0 t 2 = 2,5 sec * Th , is * s reall y the average 
sp£ptf of fall and can be calculated fr6m the formula 

(average speed) = — • J 

.where Ah is the distance fallen in the time interval At. 
See Figure 2.2. \ 

*■ f 
*}\ TVs usual, we can write 

Ah = h 2 V hi ^ 

1 ^Adapted by the UMAP Project staff from Differentiation. Second 
Edmon, 1975, Project CALC, Education Development Center, Inc., 
Newton, Massachusetts, pp. 63-75. 
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Figure 2.1. Distance-as a function of time for a falling object. 
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Figure 2*2. Enlarged portion of Figure 2.1 
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where h^ is the value of h when t = (which is 1.5 sec 
,in this example) and is the value' of h'when t = t^ 
(= 2.5 sec in this example). We also have 



» At - t 2 - t x 

so thaV 

(average speed) 



h 2 " h l 
thijf 



1 



We have seen equations like thife in the previous unit. 

What is new here is that now we* can calculate h^^ and h 2 

from a formul'a instead of estimating them fronua graph. 

Thus since h^ = 16^^ ,* we have ^ 

\ • 

h : = Ht 2 = 16(;,5) 2 .- 36 feet. 

t 0 *V 

In ju§jt t the same way, 

"* = 16t 2 2 = 16 (2. 5) 2 = 100 feet.V 

"^Accordingly^ 1 1 ' 

' h h 

" . (a^rage^speed) = f^= ^1—^= fl^= 64 ft/sec. 

Objective 2 : To be ab£# to approximate the instantaneous 

\#ate* of\changt of a function gii\en by a 
formula by taking the average rate of change ^ 
* over smaller and smaller intervals. * ■ 

The prpcedjtfe described in the previous section is 
useful when we want to -calculate an average speed. The 
result is^almost certainly more .prec\ise than any we could 
read from a graph. But what we, are really a^fter is not the 
average speed over an interval, but the sjpeed al some par- 
ticular instant of time (what we haye called the instantane- 
ous sp^ed) . In Appendix 1, we calculated instantaneous % 
speed *by finding average speed over shorter and shorter* * 
intervals. What we did there^ can be summarized in the folr * 
flowing way, * 



c 



%The process that defines instantaneous speed may be 
* ; ^^presseci * in. compact form by the following word equations: 

>alue approa^ed by [--agej ..'^o^ [i^tantaneous] 

or . 

value approached by £ v av J as At 0 - v. 

Numerical Approximation of Instantaneous Spewed 

t 

We will again approximate the speed of a falling body 
at t 1.5 seconds, but this time we will begin with the 
formula for distance of fall as a function of time, h * 16t 2 , 
instead of starting with the graph as before.. ^ 

, Throughout the calculation, h^ 36 ft and t^ = 1.5 
sec. Table 1 shows the results. The first column gives the 
Jtime at which the instantaneous speed is to be .found, 1.5 
sec. The second column* gives the interval of time over, 
which the average speed is*to be calculated. As we move 
down the table, wa let t\\e value of At get smaller and 
smaller. Column 3 gives t]\e value of t? 2 = + At, which is 
the time at»the end of the int e r val. Solving ^this simple 
equation for t 2 gives ' % * v 



*2 = 



t 1 + At,. 



Jn our, calculation t-Zis always, equal to 1.5 sec, so we 
can write this last expression , 

/ < **2~~~ + ^2 an ^ At* both in seconds). 

Calculation of Row 1 o^Table 1 ^ r% 

— \ I— • ~ ' 

Follow the calculation of the numbers in Table 1 by 

examining in detail the second row (which is typical of all 

rows in the table). Column^ of^this row is, of "cOulsse, 

t^ = 1.5 sec. In Column 2 of the second row we have t??e 

number 0.80 sec. This mean's that the Second row corres> 

ponds to the choice At ■ 0.80 sec. Column 3 gives the 

value of t 2 which is gotten from the last equation above: 
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* # TABLE** 1 

Calculation of Average Speed 
*\0ver Shorter and Shorter Time Intervals 



/ (1) • 


' (2) 


(J) 


(k) 


(5) 


(6) 


Time at 
Beg i nn i ng + 
of -Interval 
(sec) 


Length of 
. „ Time ■ 
1 nterva 1 
(sec) * 


Time at 
■ End of 
Interval 
(sec) 


Height at 
End of 
Intervaflflfe 
(ft) lP 


Change 

i n 
Height 
(ft) 


Average 
Spewed Over 
1 nterva 1 
(ft/sec) 


\\ 


At # = 


_ ^2 


h 2 1 ' 


~ Ah 


Ah/At 


1.5 — + 


\ .00 


■2.50 


. 100 


6/f * 




*T5 + 


0.80 


2.30 




48.64 , 


60.8 


^ 1.5 + 


0.60 


2.10 


70.56 * • 


34.56 


57.6 


1.5 + 


o.Ao 










' 1.5 + 


0.20 


1.70 


46.24 


10. 24 


51.2 


J. 5 + 


0.10 


1 .60 


kb.se 


Jk96 


^ 49.6 


t 1.5 


0.01 


1.5K 


36.4816* 


0,4816 


• 48.16 


1.5 . + 


0.001 


1.501 


36.04802 t 


0.04802 


48.02 


< 

putting At 


= 6.80 we 


find t 2 = 


s 1.5 + 0.80 = 2.30, 


the'num- 



ber giv,en*in the second bow of Column 3 of the table. 

Column 4 of Table 1 is h 2 ,. the distance the object has 

'fallen by - 2.3 set/ 'This va^ue is^obtained* from our- 

formula for h in terms of t: c h = 16t 2 . Thus when 

♦ 

t = t 2 = 2.3 sec, we get *^ ' * 

• .* 

-h = h 2 = 16(2. 3) 2 = 84.64 ft 

and you can see this number in the second row of Column 4. 
These* numbers , as well as those we will # calculate next, are 
shown in Figure 5.3, which is exactly like Figure 2.2 
except that it corresponds to At = 0.80 $ec rather than 1 
sec. t , - 

Column 5 gives h which, \by definition is * ' 

* • 
Ah = h; 



197 



198 



i 

A- 34 





hi 




h (fatt) 
























































*- 


























































































s 


N 
























HA 


























, 


























































































7 
















J 


























/ 


V 










































/ 


/ 












































'/ 












































// 


























cn 
OV 










I 






/j 






































\ 




/, 






























4 




















































• 
















0 3 










L 


_ 










> 














































































































































































: — r 


s\ 


































































































































Msec) 



1.0 



15 



2.3 2.5 



Figure 2.3 f * 

Since h- ■ 36 feet throughout our calculation this equation 
becomes^. » - 

' Ah = h^ - 36 (h^ and Ah both in feeH)* • 

From the previous calculation (Column 3) we know that 

h 2 * 84.64 ft. • It follows then that for row 2 of 'Table 1 

Ah ■ 84.64 - 36 ■ 48.64 ft u p • 

and this result' appears in Column 5. ^ 

The last column, of the* table (Column 6) 4 is the average 
speed of the falling object over* an interval of time . » 
At ■ 0*80 sec* beginning at t ■ t^ ■ 1.5 sec. Using 
At = 0.80 sec (Column 2) and Ah, - .48 . 64 *£t (Column 5) we 

0 ' ' 1 

; Ah 48.64 ft _ , n VfWc, N ^ 
\ S Jt O.ffO sec " 60 ' 8 ' ft/seC .' . 



'the number* listed in Column 6 o*f Table 1. 
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1.93 



Air* of Table 1 is made* this way, each value of At 
(Column 2) being chbsen sm^ler than the value in the pre- 
vious row. Spot check a few)of the numbers in^other rows of 
Table .1 to be sure you understand, how each is, calculated . 



Exercise 1. The fourth row of Table .1 (corresponding to At =* 0.40 sec) 
has been left incomplete. Follow the procedure discussed above and 
fill in the row. Consult Table 2 to see if the value you obtain for 
Ah/At Is correct. 

■» * ' 

Results of the Calculation- as At ■» 0 ' 

Now that we have Table 1, we can 'make n secondT much 
-simpler table, by omitting everything except the At and, 
Ah/At columns, as is dorfe in'TabJ.e 2. Since we are. 
interested in the average speed, Ah/At, for smaller, and^ 
smaller values of At, these^are the important columns of 
Table 1; the other columns were ^ut in only to assitt us in 
making the calculations. 



^ TABLE 2 9 • 

The 'Average' Speed 
Over Shorter and Shorter Time intervals 





Length of 


^Average Speed 




Time interval 


Over interval 




(sec) 


(ft/sec) . 




At 


Ah/At 




1 .00 






y 0.80 


6G\8 




*0.60 * 


57.6 


• 


0.40 


5k.k 




> 0.20 


51 .2 




07.10 






- 0.01' 


48.16 


* 


0.001 


48.02 


» 
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The values .of average speed from Table 2 are pJLotted 
in Figure 2.4. By looking at either the table., or the 
graph it is possible to see that, as At grows smaller and 
app^foahces 0, Ah/At gets closer and closer to a va~lue at 
or near 48 ft/sec. In a word equation: ^ 

value approached by j£ as At + 0 = 48 ft/sec. ' % 



Since we already have stated that 
value approached by ^ as At * 0 = 



[instantaneous] 
speed " J 



■ y» 



the value 



v-= 4fr ft/sec 

is our approximation f&r the instantaneous speed of fall * 
at l t> 1.5 sec. 



At , 



55 



54 



Avtragt 

SpM 

53 

tnttrval 

52 



51 



50 



49 



48 



47 



: cm jljtll 

: mill 





















































.1 j2 * A 

Length of Inttival At 

Figure -2. A*. Some values of Table 2 plotted to see what happens to 
Ah 



At 



as At * 0. « 
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Exercise 2 . Compare^ the results of the calculations shown in Table 2 
with your calculations of average speed using the tangent line done 
in Table 1 o^ppendix 1. ^ 

Exercise 3 . Use the 'formula h = l6t 2 and the numerical method f just' 
described to find (at* least Approximately) the speed, of the falling 
object at t ■ 2 sec. • * 

Exercise* k (Calculator Exercise) , IT you have a calculator a,t your 
disposal, redo the< calcul a't ion of Exercise 3 to find the speed of the 
falling-object at t = 0, 0.5, 1, 1.5, 2, 2.5, 3 sec, eto. Plot the 
speed as a function of- t. What kind of curve is\our graph? 



Another Example; Calculation of a Slope at a Point 

Let -us try to find the elope pi the graph of the * 
function y = x 3 at x = 2. You will see that although in 
the previous example we calculated an instantaneous velocity, 
and in this second exampLe we wiLl calculate a slope 3t and 
the' procedures involved are identical. * • 

In the last example, we calculated instantaneous ' 
velocity from average velocity.' In this example, we carets 
late' the instantaneous rate of change o,f a function (the 
slope of its graph at a point) -from the average ra.te of 
change. Definitions* and procedures are identical and may 
b.e expressed in, the following form:* N 

•/ 

value approached by 



ayerage 
rate^of 
'change 



'as Ax + 0 * [ instantaneous 1 
|_rat.e of change] 9 



, value approached by \&\ as = slope 'of curve 

|_AxJ at a point 

Figure 2.5 is a graph of the function y ~ x 3 . You 
should make a few calculations yourself, and verify that . % 
the curve in Figure 2. 5^ really is a graph of*y = x 3 . (Fig- 
ure 2.5 shows ,y = $ 3 only for positive values of x. Can 
you iiraw a graph of y' = x 3 for negative values of x, say 
from 0- to -5?) * 1 
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12,34-5 
Figure 2.5. Graph of y ^ x^* 

Numerical Calculation of the Slope at x 2 

To find the slope of the curve in .Figure 2,5 at x = 2,» 
we proceed as we did above in calculating instantaneous ■ 
velocity. That is, we calculate ■ . 

Ay _ y 2 " ?l 

Ax " x 2 - x x ^. 

\ ' - 

for smaller and smaller values of Ax (where Ax s x 7 - x« ) . 
In our example, x^^ = 2 and so = 2^ *= 8. Hence, 



1 .Xv " Ay 



>Ax Ax Ax 

1 

The results of such a calculation are shown in Table 3.* • 
This table is jmade exactly the way we made Table 1, so we 1 
will not discuss -it in detail*, but you should spot cfoeck a 
few of the -numbers. (We. have rounded off the numbers in 
Tjible 3 so* as t& 'keep no more than two figures after' the 
decimal point in^rWATx.) * * • 




' We *now have a method for finding instantaneous rates' 
of change. Speed, or more correctly instantaneous j£ee$ 9 is J 
an example of this. Practice this method on the next few 
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TABLE 3 

Approximating the Slope of y = x^ 
at x * 2 by Letting Ax Approach Zero 



3 



1 


+ 


Ax 


* 2. 




Av = v - R 

a 7 72 0 


Ay/ Ax 


*2 


+ 


1 .0 




0*7 
£■/ ■ 


1 

19.0 


19 


2 * 


+ 


0.50 ' 


2.50 


15.625 


7.625 


15.25 


2 


+ 


0.20 


« 2.20 


10.648 


2.648 


13.24 


2 


+ 


i).10 


= 2.10 


-9-261 -> 


1 .261 


12.61 


. 2 


+ 


0.05 


= 2.05 


8.6151 


\o.615 


12.30 


2 * 




0.02. 


= . 2.02 


8.2424 


0.2424 


12.12 


,2 


+ 


0.01 


St 








2 




0.005 
/ 
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Exercise 5 . Complete the last two rows of Table 3 to verify that 
Ay/ Ax approaches 12 as Ax approaches zero. 

Exercise 6 . Verify the answer of Exercise 5 by drawing a line tan- 
gent to the curve in Figure 2.5 at x ■ 2, agd measuring its slope. 

m 

Exercise 7 . Adapt the method used above to find the slope of 

2 1 
y = x at x * 1 . 

Exercise 8 (Calculator Exercise). Use the methods of this unit to 
find the slope of y = x 2 at x = 0, 0.5, 1, 1.5, 2, 2.5, and 3 sec. 
Graph the slope versus x. What kind of curve do you get? 



Some Further Comments 

It may have seemed that we found both instantaneous 
speed and the slope at a point in a rather roundabout way. 
In each case we calculated an' average rate of change (Ah/At *\ • ' 
or Ay/Ax) over an interval (At or Ax) which got closer and 
and closer to zero. But why be so sneaky? Why not just 
set At or Ax equal to zero instead of letting it approach 
zero? 
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' To anstf«£ this question, consider -the IgLnd of^ * • 
measurements we make in order to' get the numbers displayed 
in. Table 1. Put «in its simplest terms, we determine the 
distance, Ah, which an object falls "in a'fcime At. We then 
calculate Ah/At. Putting At - 0* directly would mean that 
\ we wou^d have to determine the, distance/ the object falls 
: in a time "interval" of zero leng,th^ But in.zero time . 
the object' moves zero distance so .that corresponding "to 
At = 0 we would have Ah 55 0 , and the ratio Ah/At = 0/0. 
The additional row for Table l'whicji woul<l correspond to 
this -calculation is- shown in Table 4.* < 

TABLE k 

* An Additional Row for Table 1 • 

• Corresponding to At » 0 *- * 



tj(sec) 


+ t(sec) 


",t 2 (sfec) , 


, h 2 (ft) 


Ah(ft) 


Ah/At (ft/sec) 


1:5 


+ ' 0 


= 1 .50 

* * 


36 \ 


0 


** » 

0 1 
/ 0 • 




There are tw^ 'things .wrong with the "result", 0/0 n 
First, ix^BWiifc same result we would get for the 'object 
whatever its true speed 4 might be. Whether it is moving so 
•fast that all we see is a blur as it passes, or so slowly 
th^at it is^ at the proverbial "snail 1 s pace," ou 4 f "answer" 
will be 0/0. Obviously a quantity whith does n'6t distin- 
guish between something moving rapidly' and something* . 
- moving slowly can hardly bemused as the definition of * 

instantaneous velocity. * 

> / 

The second difficulty with* 0/0 is a mathematical one. 
0/0 is, mathematically speaking \ £ meaningless symbol; if 
is impossible to ascribe a unique numerical val.ue to it. 

* V 

4— 

All thi£ helps to explain why we must "-sneak up" on 

the instantaneous velocity* by determining the value Ah/At 

approaches a£ At ^approaches zero.. We can summarize our 
x t . >* 

^ * ' A- 



discussion of instantaneous velocity by saying that 
Ah 



At/ 



v as* At 0. 



Here v is thej instantaneous velocity at some specific 

time, At is an interval which begins at that time , and Ah 

is the distance the object, m6ves during that interval. 

ANSWERS -TO "EXERCISES 









— 5 ■ 


tj + At t 2 


h 2 


Ah 


Ah/At 


1.5 + 0.A0 ± 1.90 


57.76 

r 


21.76 





2. The two results should 
slopes by the tangent 1 



3. At t = 2 sec, 
v » (>k ft/sec. 

k. The graph is a 
straight line. 



100 



90 



80 



70 



60 



50 



40 



30 



.20* 



agree except for small errors in finding 
ine method. 

44+44 



(ft/*c) 




tin* (Mcondt) 
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■ 1 ; 

Xj + Ax' ■ x^ 


y 2 


Ay ' 


' A / A 

• » Ay/ Ax 


» 

2. + 0.01 = 2.01 , 

% 2 + 0.005 = 2.005 
• 


8.120601 
8.0601501' 


*e — 

o.r2oeor 

0.0601501 . 


: 1.2,06.01 
12. 03002s 



6. Thertwo results should agree except for smal'l errors in finding 
slopes by the tangent 1 ine, method. v * . 



^ 7. 



X, + 


Ax . = 


s * 2 \ 


• y 2 


Ay = y 2 - 1 


'Ay /Ax 


1 + 
1 + 
1 . + 


Oil - 
0% = 
0.001 = 


1*1 

«/ 1.01 
- 1.001 


l f .21 

1.0201 

1.002001* 


0.Z1 

0.0201 

0.002001 


2.1 

2.01 

2.001 



8. The graph is a straight linei 




207 



.APPENDIX S 



Differentiation Formulas from Calculus 



1 

If u and v are dif ferentiable functions of x, then 
the sum and product functions, u + v and u • v, are 
dif ferentiable functions of x, and their derivatives axe. 
given by the" formulas 

^_( U ♦ V ) = g + .Jv (Sum Rule) . 



and 



-g^-(uv) = u^~ + u^-. % (Product ^ule) 



Further, if, v(x) ? 0, then the quotient function u/v is 
dif ferentiable at and its derivative .is*, given by the 
formula 1 ^£ 

£(H] . v(du/dx) - z ujdv/dx) . (Quoti§nt * RuleX 

Finally, if y is a dif ferentiable function of u, and *u 
is a dif ferentiable function of x, then y*is a differen- 
tiate function of x and 
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APPENDIX 4 



TJie formulas for the derivatives of the trigonometric 

functions 0 of x are listed below. ^ 
f ' 

Bin x r * cos x, ^ cot x = -esc 2 x, 

cos x » -sin x. ^ sec x = sec x tan x. 

^ tan x ■ sec 2 x. ^ esc x = -esc x cot x. 

The formulas for the derivatives with .respect to 
x for trigonometric functions of'u where u is a 
diff erentiable function of x are listed below. 

d du d ' 2 du 

• £ sin U = cos u gj. ^ cot u = -esc 2 u ^. 

d du d du 

cos u ■ -sin u sec u = sec u tan u . 

d 2 du d * *. du 

£ tan u = sec 2 u ^ ^ esc u = -esc u cot u ^. 



< ■ . >, 

V> ANSWERS TO* M ODEL EXAMS - " 
Answers to Model Exam from -Unit 158 

„./ \1. 4 Caflculus is not needed. * The answer is 
* t h = 1200 tan 0.3 * 371.2 meters. 

2. ' Calculus is needed .'^ The first and second 

. deriVativ 4 es -of voyage function 'will be f o uid 
and the second derivative test w-ill \>e us^ed. 
The maximum Voltage is v = 1, wfien t = ~. 

3. Calculus is not needed. The formula for the 

r area of a sector will give 'approximately 51 in. 2 . 

4. falculus ±1, needed*- The area of the triangle, 
A = 18 sin 6 cos 6, is obtained by right triangle 
trigonometry. Then, as in problem 2, the second, 
derivative test will be used. The answer is 

a - * 

e - 4-.. • • ; « , 

5. A problem similar to any of the .four problems 
given in the unit or in Problems 2 and 4 of this 

x exam would be an acceptable answer. ' 
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Answers to Model Exam from Unit 159 

a * 37 ^ sin " cos x 
m 3x ' cos x ) = " Sln x * 

2. radian. . 

9 

3. Yes. The result wilTfcot hold if degree measure 
is used'.- 



5. 



a. -0.7 
b- -0.9 






* 


t 


i 


Ax 


■ Ay = sin(0.4+ Ax) - sin0.4 


Ay/Ax 


.1 


.0900 


.900, 


.01 


.t)0919 


.919 


.001 " 


.000921 A 


.921 


.0001 


.000921 


.921 


-.1 x 


-.0939 


-.939* 


5 . 01 

• * 


-,00923 


.'923 


-.oca 


-.00092^ c 


* S .921- 


- .0001 

V 


- .0000921 


-.921-* * 



The value of the derivative y = sin x at x =0.4' 
is approximately 0.9,21. 

* « 
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Answers to Model Exam from Unit 160 

■ — . * 

1. * S *i ^ is sandwiched between two other expressions, 
< n • ' 

each of which has a limit of 1 as h + Of. 



2. 


l im cos h - -1 . 
h+0 


4 




3. 


cos (~) ' 






4. 


-sin (f§£), - 


* 




9* 


2x cos (x 2 - 3) . ' 0 


V 




6 . 


2 rn< y — *? c i n ^ y 






7. 


*»2 sin x cos x 






8. 


-sin 2 x + cos 2 x,. 






* 9. 


-(6x 1) sin (3x 2 - x) 






io; 


cos 3 x - 2sin 2 x cos x 








-cos x "+ c. 






12. 


\ 

j sin 3x + c . ' > 






13. 


5x + 2cos x + c . 






14. 
• 


dF _ -KW (K cos 8 - sin 


8) . 


* 


^' (K sin 6 + cos* 6) 2 








> 





-Answers to Model Exam frofa Unit 1^1 
d _ d /sin x 



X'. 
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-2-(tan x) = 4-( sin x > =» cros X • cos x - sin x ( »sin x) 

<, , QXr COS. X • » 9 ^ 

cos' x ^ fc 



d: 



cos 2 x + sin 2 * & 1 v 

* *cos 2 x -cos 2 x 



sec x. 



ax*" ' dx^sin 



sin x ^ 0^<- »1 ■ cos ,x - 



-cos x 



sin' x 



sin ' x 
cos x 



siji x sin x 



-esc x cot x . 



CSC' X, 



* 3. 

5. j}£ = 3 sec («3* + .5). 



6 ' 2x = 2 cos 2x + 2x s ^ c * x2 '- 



V 



x sec* x + tan x. 



*g h_ = -sin x (j + tan 2x) - 2 cos x tan 2x 

- " d< + tan 2x)V4 



/ 
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STUDENT FORM 1 
Request for Help 



Return to: 

EDC/UMAP 

55 Chapel St. 

Newton, MA 02160 



Student: 



— + — u7 r If#yo ¥ have trouble with- a specific part* of this unit, please fill 
out this form- and take it to your instructor for assistance. The information 
you give will help the author to revise the unit. ) 



Your Name 



Page 



T 



0 -Upper 
OMiddle 
O Lower 



Unit No. 



OR 



Section 



Paragraph 



i i ^ , |_ 

Description of Difficulty: ' (Please be specific) 



OR 



Mo^el Exam 7 
'Problem No. 



Text 
Problem No. 



Instructor: Please indicate your resolution of the difficulty in this box. 



o 



Corrected errors in materials. List corrections her6: 



o 



Gave student better explanation*; example, 'or procedure than in uniw^ 
Give brief outline of your addition here: 



Assisted student in acquiring general learning and protflem-solving 
skills (not. using examples from this unit.) 



Instructor's Signature 
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Please use reverse if necessary. 



* STUDENT FORM 2 . . , . ' E^AJMaT 

'* Unit Questionnaire 55 Chapel St. 

„ • Newton, MA 02f60 

Name •. • . » 
• . 1 ■ : _Unit No. Dat e 1 

' In8 titution_j . : C ourse No : ' 

Check the choice for each questiol tha*Ws closest to your personal opinio*. 
1 - How useful Vas the amou nt of detail in "the unit ? , 

, Not enough detail to understand, the unit ; A 

Unit would have been clearer with more detail » 

* —^Appropriate amount of detail . fi 

J**' Wa * occas .i°nally too detailed, but this was not distracting 

, - Too much detail; I was often/ distracted. 

2 ' How helpful were the problem answers ? «. 

?S?^ SO J U f ± r iS t0 ° br±ef; 1 C ° Uld not do intermediate steps 

Sufficient information was given to solve the problems 

Sample solutions were too detailed; I didn't need/ them 

' v 3V T^J° T fulflUln « ^e prerequisites, how much did vou use other sources (for 
,\ jgiffl- Ple, instructor , frien d s, or other books) i n orde/to understand the unit? " 
__1A Lot ^Somewhat ' . A L ' ittle . Not at 

How long was this unit in comparison to the amount of time you generally spend on ' 
a lesson (lecture and homework assignment) in a typical math or science course? 



4, 



Much Somewhat About ■ Somewhat Much 

.Longer Longer the Same , Shorter Shorter ^ 



P~ Were any of the following parts of the uftit confusing or det racting? (Check 
' as many as apply.) : \ — TP — T tt 

P rerequisites ■ * ' \ 

Statement of skills and concepts. (objectives) 



^Paragraph headings 
^Examples ^ 
^Special Assistance Supplement (if pres^nff 
JJther^ please' explain^ ' , \ 



Were any of the following parts of the unit barticulLly helpful? (Check as many 
as apply,) ^ ; c J 

Prerequisites * \ •* 

Statement of skills and concepts (objectives) * - 

Examples 

Problems . 

Paragraph headings 

Table of Contents 



Special Assistance Supplement (if present) 
J)ther, please explain^ " 



Please describe anything in the unit that you did not particularly like. 



•4 

please describe anything that you found particularly helpful. <PleaWlise^he back of 
^this sheetN if you need more space.) *■ 
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DETERMINING CONSTANTS OF INTEGRATION 



4 Ross L.. Finney . 
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Cambridge, .MA 02139 



INTRODUCTION 



As you know, there are times when we have information 
about the derivative of a function and wish to conclude 
front it ^information* .about tjie function itself. 

n * DERIVATIVE FUNCTION 

Velocity Distance 

Acceleration Velocity 

-Marginal cost Cost 

Ra<e of growth Size of the „ * 

of a population population 

The reason tHat -we can often succeed in determining func- 
tions from their derivatives is that whenever two functions 
have the same derivative on an interval, the functions dif- 
fef only by a constant on the interval. Thus, if we can 
find even one function that has the given derivative, we 
know that the function we seek cannot" differ from it by more 
than a constant. The basic fact is this: 

IF f» (x) = g ! (x) v for ALL VALUES OF x - 
.IN SOME INTERVAL, THE^J FOR SOME CONSTANT C 

" gOO = C or f(x) » g(x) i C 
FOR ALL VALUES OF x IN THE INTERVAL. 
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y slope 2x 



y = x 



slope 2x 



"For every value of x, 
the two functions 
f(x) « x 2 + 1 and 
g(x) « x 2 - 2 
have' the derivative 
f f (x) a g' (x) ^ 2x. 
Wotice that 
'f(x) = g(x) + 3 
for 'all x. The 1 
value of C in the 
rule stated above 
is C - 3. 

To get the graph 
of f , we may si ide the 
graph «of y ■ x 2 tip 1 
unit. To get the • - 

c ' V 

graph of g,, we sli<£e it down 2 units. The three graphs have the same 
slope at any x. 

'Functions whose derivatives are equal only at isolated 
points, however, do not have to differ by a constant. 

EXAMPLE 1. The difference of the functions f(x) = 2x 2 and 
^g(x) ■ x 2 is 2x* - x 2 - x 2 , and not a constant. However, the 
derivatives of these two functions have the same value at x = 0, as 
you can see in the following table. 




THE 
FUNCTIONS 



fix) 

g(x) 



2x 2 

v 2 



THEIR 
DERIVATIVES 

f(x) =,4x 
g'(x) = 2x 



THEIR DERIVATIVES 
AT x 0 

/\ 

0 
0 



Exerci ses 



Find two values of x at which the derivatives of f(x) * 2x 3 and 
g(x) » 3x 2 are equal. , 

/ 1 \ 



. 2 
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2. Suppose that f (x) and g(x) are two functions that have deriva- 
tives on some interval, and that 

f(x) - g(x) = C 

* * on the interval. ( J* 

a) Differentiate both sides of the preceding equation 
to show that dif ferentiable 'functions' that differ by 
a constant on an interval have the rtme derivatives 
on the, interval . 

b) -Show that 

2x 3 - 3x 2 
is not constant on any Interval. 

3. Find two more functions whose difference is not a constant 
but whose derivatives agree at one or more points. 

In Exercises k and 5, use the coordinate axes provided to graph 
the given functions.. 

• •» 

k. Graph the lines y * x, y * x - 2, 
y » x + 1, and y = x + 3, 



1 



71 1 1 1 f 



-j I 1 1 h 



I 



2Sl 



5. Graph the cubic curves 
y * x 3 , y « x 3 - 0.6, 
.and y » x 3 + 0.4. 



1-» 



/ 



H 1 1 1 h 



I I I I 



2. INDEFINITE INTEGRALS , v 

t 

Since the derivative of 5x 2 is lOx, ariy function f(x) 
that has- the derivative 

$'(x) = lOx. 1 

must have the forny 

f (x) = 5x 2 + C - 

for some constant C. Without more information we ca 
learn the Value of C, but at least* we have detevrtined f up 
to a constant, as we say* We call the family of functions 
5x 2 + C the indefinite integral of lOx, and we show this by 
writing * 

jlOx dx = 5x 2 + .C. 



cahoot 



The constant C in this formula is called the constant of 
tntegration. 1 • * 

Many— indef inite integrals may be found by reversing 
derivative formulas we already know. Here are some examples. 

"4 
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DERIVATIVE 'COMPANION 
FORMULA INTEGRAL FORMULA • # 

U d (kx)wk . r.<*kdx«kx + C 

dx * % % 

If we change all the r 1 ^ in formulas (2), and (2. f i to 
r + 1, and the^i divide both sides of the formulas so 
obtained by (r + 1), we get formulas (3} and (3 1 ) shown 
below. m Formula (3 1 ) tends to be more useful than form- 
ula (2»). 

3. ± (±11) „ / V. f* r dx « + C 

dx V + V * > r + 1 

-Formulas (3) and (3 1 ) don't work when r .= -1, but the ' 
next formula takes care of this case. • 

h. d_ 
dx 



(In |x|) = i V. j^dx - In |x| + C 1 

5 ' £ (e X ) -e x 5'. fe X dx f e x + C 

Notice how nice a furrction e x is! 



\ 

EXAMPLE 2. 



J-5 dx - -5x + C 




EXERCISES 

Complete the equations ?n Exercises 6 - 21. 
. ** 

6. J* dx - 7. J-2 5 dx 

•••J— 



dx » 17x + C 9. dx » -3x + C 



10. |x 2 dx - _^ 11. |x 5 dx - 



f 



dx - 



i*2.[f -ln|x| + ^ 13. (I 

I*- [ dx - £ + C . 15. [y.dx - _ 

v r 102 { 

16. M03x dx - * 17. dx 

{ dx = e x + C 19. | ^ x = x> + C ' 

20. | d x « jtC 21 . j_ 



x 17 

V + c 



_dx ■ -103. Sx + C 



In Exercises 22-25, the lette.rs a, b, k and m are constants. 
Complete each formula. 

22. jk dx » _J 23. Ja dx = 

2! »- I. dx - mx + C 25. L dx - -bx + C 

So far we have used x as the only variable of integration, but 
other letters are commonly used." Complete the equations in 
Exercises 26-33. 



|32 dt - ^ ^27. f 

28. Jp 2 dP - . 29. d® = r 3 



dt » at + C 



dB ■ B + C 



30. Jss* ds - ^31. |v dv - 

32. je y dy - i N 33. j , d R « In |r| + C 

In Exercis£^34-39, the letters with subscripts are constants. 



Complete each formula. 
3^ |t 0 dt - 35. j_ 

36. |a 0 dt - » -37. | Vq dt = _ 



dt » 32t 0 t + C 



38- I dy - y 0 y + C 39. j d t » s Q t + C 



1 ^ 
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In Exercises 40-43, use the coordinate axes provided to graph 
the* three curves selected from each family. 



40. Graph y -2x + C 
r for C » 0, 3, and -4. 



ay 



H h 



H 1 I- 



H r- 



X 



41, * Graph y » + C for 



* > 



C » -0, -1 , and 



9 „ 



H H — I 1 h 




H 1 1 H 
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k2. Graph y » e x + C for 
C » -2.718, - 1, and 0. 



H 1 1 1 h 



H 1 1 1 h 



43- 1 Graph y*ln|x|+C, x?*0 
for C - 0, -1 , and 1 . 



4y 



1* 



-+ 1- 



H h 



-h h 



H h 



3. INTEGRALS OF LINEAR COMBINATIONS TDF PACTIONS 

Tne rules about differentiating sums of functions 
and* multiples, functions lead to the ^following two 
integration formulas : A „ - 
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(6) SUM RULE |[f(x) + g(x)] dx + |g(x) dx 

(7) SCALAR MULT I PIE RULE * jk f (x) dx » k |f (x) dx 

an y constant). 

The fact that the integral of the negative of a function 
is the negative of its integral is an ^immediate conse- 
quence of (7). We just take k - -1 to get 

« 

, J-»f (x) dx * j-1 ■ f(x) dx - -1 • jf (x) dx « rj>(x) dx. 



— EXAMPLE 3- 



^ J-x* dx « -|x 2 dx » - y- + C 



The' sum and scalar multiple rules are often used 
together, as in the next example. 

e/aj^le k. 

J (y* - J- + 9)dy - fy'dy - J} dy + |9dy j^i RULE' ^ 



,5 



-^--1n|y|+9y+C SCALAR MULTIPLE RULE 

In general, the sum and scalar multiple rules allow 
us to break problems into parts we know how to solve (we 
hope) . „ 

When we integrate a sum or difference of functions; 
one constant of integration is enough to generate the whole 
family of possible solutions* * 

EXAMPLE 5. * ' 



|(3t 2 + 12e t )dt - |3t 2 dt + uje' dt 
\ ' = t'+Ue* + C. 



' ' t 

We dp not need to wri'te^the answer ast 3 +C^ + 12e + 12C 2 

The formulas t 3 > 12e t + C and t 3 + C 1 + 12c 1 + 12C 2 may 

generate different functions for different values of the C 1 

but the family of functions generated by either formula is 

< 

9 
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.• the same as *the family generated by the other. We are 
therefore free to use the simpler formula, which -is what 
we do. 



EXERCISES . . 

Complete the equations in Exercises kk - 55. 



50. 
52. 
5*. 



(x + l)dx - y + 



-e dx 



jdx - — - x + C kl 



• | dx » -In | x | + C 



(mx + b)d 



y 



ds 



(1 . * + i,e*)dz 



49. J-t 2 dt - - 

51- |(3y 2 - 5y)dy = 

A 

55 



s 3 +C 53. J(lx -i e x )dx 



| dx « x 3 - x 2 + 7x + C 



4., FRACTIONAL EXPONENTS AND NEGATIVE EXPONENTS 



In this section we return to formula (3») of 
Section 2, which we now call the power formula. 

(8) POWER FORMULA * | x r dx . r f -] 

We do tfiis to point out that the exponent r in the 
formula does nat. have to be positive. It also does 
not have to be an integer. 



EXAMPLE 6. 



[**dx - 4 



|x* + C 



EXAMPLE 7. 



fx-dx 



USE THE POWER FORMULA 
WITH r * J. 

SIMPLIFY 



WRITE THE INTEGRAND WITH A 
NEGATIVE EXPONENT. USE THE 
POWER FORMULA WITH r - -k. 



10 
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+ c 



INTEGRATE 



+ C 



POSITIVE DENOMINATOR 



3x3 



+ C 



EXAMPLE 8. 



- -15 jxT 2 dx 



•15- 



+ C 



POSITIVE EXPONENT 



USE THE SCALAR MULTIPLE RULE 
TO SIMPLIFY THE INTEGRANO. 

WRITE THE INTEGRANO IN EXPONENTIAL 
FORM. TAKE r - -2 IN THE POWER 
FORMULA. 

A 

INTEGRATE 



15x" J + C 



SIMPLIFY 



II + 



POSITIVE EXPONENT 



EXERCISES r 

Complete the equations in Exercises 56-80. 
56. 



58. 
60. 
62'. 
64. 

'68. 
70. 
72. 
Ik. 



At 



-12»T ds 



Jlx-Jf+C 



Jt - fr + c 



-dx - 

x - 



dx- x + In |x| + C 



(x 2 + ^)dx 



>57. 
59. 
61. 
63. 
65. 
67. 
69 .« 
71. 
73. 
75* 



dx ■ 2x + C 



5*^x dx ■ 
'j^/y dy 
dx ■ 



dx - - + C 
- x 



jds « -j- + C 



dv 



X 



dv - v + — + C 
- v 



(x +1) dx ■ ' 



11' 
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? 6 ' } dz-e 2 + H 2 s | +C 77. {(e X -^)dx-"l 

80. J| /y dy - 



S_. DIFFERENTIA^ TO CHECK 

To check an indefinite integration there 'are^o 
steps to follow: ' . / * 

1. Make sure a constant of integration 
is there. 

2. Differentiate, to see if you get back 
the function you integrated. 

EXAMPLE 9. 

The equation 

• . |(2x 2 - 5x")dx - | x 3 - x 5 + C 

is correct because 
' (1) C is there, and " ■ 

<2) ^(fx 3 -x^ + C). 3 -fx^5xS^O 
- 2x 2 - 5x\. - • 

i 

EXERCISES 

True, or false? Differentiate to find ogt. - 

31. Jx'dx - £ + C 82. |x 2 dx * x 3 + C 

, 83. jt J dt - | t* + C 84. foj dy - y* + C 

85. |(e s -'l) ds - e s - 1 + C 86. Pj. dx - - I + C 



)x' x 
Find the missing integrands in Exercises 87-98 



87. | dx - x 5 + C 88. | dt - ft 

* »• (- f»V " 7 ♦ C * 90. | dt „ 

, (fiere, v 0 and s Q are constants.) 



I6.t 2 + c 



,6t l + V + H 
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3K I -Jy - »n Ivl ♦ C 92, | 

, dt " 32 1 * v 0 

' (v 0 is a constant.) * 



93- | dx - mx + b $k. J dx - x* 



dx ■ x - x + C 



(m and b are constants.) 
95- f IdT - M + To g 6 . | - 23 e* ♦ C 



(T 0 is a constant.) 
57. |_J__dx» -!20x 2 -7px+C 98. I . d v - I v* 



+ C 



6. NUMERICAL CONDITIONS vTHAT -DETERMINE ^ CONSTANT 
e OF INTEGRATION fc ' 

Every function r f whose derivative is given by 
formula , 

f'(x) = lOx 
is a member of the family of functions 

jlOx dx* » 5x 2 + C 
, But to determine just which one . f is, we need more in- 
formation. The information can ^supplied in* various 
wavs *^F°r instance, we might know the value of. f at a 
particular value of x, as in the following example. 

EXAMPLE 10, Find f if f (x) - lOx and f(l) * 3. 
SOLUTION ' * 8 

1. f (x)* ■ 5x 2 + C for some C. BECAUSE 

f'(x) - 10x. 

2. 5(i) 2 + C - 3 -* 

BECAUSE 

■ . - 5 + C " 3 " '0>-3. 

c - -2 x ^ 

Conditions like f(l) » 3 that select a particular solu- 



tion from a family given by an indefinite integral are 
called initial conditions. 



13 
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EXERCISES 

Copy and complete the table below. The first row, already com- 
plete, shows the following Information: 



The derivative 
of a function 



3. An Initial 
condl t Ion 



5. The function deter- 
mined by the deri- 
vative together with 
the Initial condi- 
tion. 



\E1 f* T ~ T ^ tft-S) " 291 



2. The indefinite 
integral of 
this derivative 

i 



T 



if (x) «x z + g 



4. The value of C de- 
termined by the given 
Initial condition 





DERIVATIVE 


INDEFINITE 
INTEGRAL 


INITIAL 
CONDITION 


VALUE 
OF C 


PARTICULAR 
SOLUTION 




f'W 


jf ■ (x)dx 




C 


f(x) 




2x 


x 2 + C 


f(-5) - 29 




f(x) - x 2 + k 
* 


99. 




In I x | + C 


f(e). = -3 






100. 


X 


• 


f(l) - 2 






101. 


-x 




f(D « o 






102. 


x 2 + 6 




f(l) - 10 






103. 


32 




f(0) « 0 












f(0)'« 5 


* — 




.105. 


32x + v 0 




f(0) - 0 







106. Find the*potential energy U(x) of an object as a function of 
Its position x, when the magnitude F(x) of the force acting 
on the object is given by 

F(x) » -kx. 



Assume that U(0) - 0 and that U(x) » F(x) dx. 
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7. GRAPHICAL CONDITIONS THAT DETERMINE A CONSTANT 
OF INTEGRATION 



The graphs of the functions 



f(x) » J(4 - 3x 2 )dx 

make a family j>f npn- 
overlapping curves in 
the plane. There is 
one curve for each 
value of C. 

Choosing a 
function from the ■ 
family amounts to 
choosing one of 
these curves. 

One way to pick 
out a curve is to 
name a point on it. 
We might say, for 
example, "Take the 
curve that passes 
through the point 
(0 t 3K M This says 
that, of all the 
.curves y ■ f (x) ; we * 
want the one that 
satisfies the ini- 
tial condition 
f(0) = 3, 



Figure 2. The graphs 
of 4x - x 3 + c, 
for C « -3, 0, 2, k. 



4x 



+ C 




EXAMPLE 11. Find f If - ; 

1. f'(x) « ^- 3x 2 . 

2. The graph of f passes through (1,5), 



SOLUTION 



'1. The values of f are given 1 by. the~formu1a 

f(x) - j(J»\3x 2 )dx - 4x - x 3 + C 

2t, The graph of f has the equation 
9 y ■ Ax * x 3 + C. 

3. Because 1 (1,5) Hes on the graph, . 
4(1) - (D 3 + C - 5 

3 + C - 5 , ' 
- C « 2. 

A. Thus, f(x) « Ax - x 3 + 2. 



EXERC l,SES 

In Exercises 107-112, find the value of C that makes the curve 
y « kx - x 3 + C pass through the given point. * 

107. (0,0) x 108. (2,0) 109. (-2,0) 

n °- <M) HI. (2,7) 112. (3,1) 



Copy and complete the table on page 17. 
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* * 


• • 




✓ 


FORMUtA 
, FOR f ' (x) 


• Jf'(x)dx 


A POINT ON 
THE GRAPH 
OF f • 


fORMULA 
FOR f(x) 


v»' 

2X 


x 2 + C 


(5,20)- . 




113. 5* 




(-2,1) 




114. 8x 




(o^5) >. 




115?** -4x + 3 




(-1,1) 




116. 9.8x 


• 


(1,3) 




117* e x - 2 




'to; 7) 




us. , i , 

A 

/ 


\ 


(1.4) 


0 


119. 




(2,0) 


- 1 


120. i/x 


• 


J * 

(3,0) 




li) Exercises 121-124, 


graph the^ function 


wffbse derivative 


is given 



and that satisfies the -given initial condition. 

121. 32t + 10 

dt . . % 

s(0) » -4. 



122. f -9.1 



123. 



v(0j - 0 

dt 56 
y(0j - 7 



124. Jf - 
1 s(0) - 0 
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8. MODELING; INITIAL CONDITIONS FROM - 
* PLAUSIBLE ASSUMPTIONS 

When initial conditions are^not stated explicitly, 
they can sometimes be. inferred from other^nformation or 
^>ased oh plausible assumptions. 

EXAMPLE 12. "To sample the upper atmosphere a rocket is fired 

straight up from the ground. The rocket engine accalerat.es 
the rocket at Wsec , and has enough fuel to' burn for 2 
minutes. 

1) How high is the rocket 1 minute after launch? 

2) How fast ijs it climbing then? 

3) How high will the rocket be when the engine stops? 

4) How fast will it be climbing then? 

ANALYSIS - 

A 

The questions on the list are not as formidable as they might 
seem at first glance because we can answer them all 'by finding for- 
mulas that describe. the rocket's heightfand speed as fuhctions of 
time. . * % 

To begin, let s(t) denote the rocket's height in meters as a 

function of time measured in seconds. The choice of meters and seconds ' 

is a natural one to make, because the rocket's acceleration is given in 

those terms. The use of the letter s Is traditional. J 

« „ 

Then s'(t) gives the rocket's velocity and s M, (t) its accelera- 
tion, so that while the engine is on 

s"(t) - k m/sec 2 . 

If we measure* time with t - D at the~time of ignition, and assume that 
the engine gives full thrust fronj the very start, then; 

s M (t) - k D < t < 12D 



and 



s'(t) = fs"(t)dt; 
[j»dt 

At + C, D S t < 12D,, 



■J* 



IS 
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meters per second being understood. Since C ■ s'(0) is the initial 
velocity of the rocket, we usually wri,te v Q in place of C, as in the. 
next equation: 

s'(t) » kt + v 0> o < t < 120. ' 

If we assume that the rocket, is fired from rest, then 

v Q o 0 and ' s'(t) « 4t when, 0 £ t £ 120*. 

j To f i/id s(t) for the two-minute interval the engine is on we 

integrate again. This gives ' 

€ *s(t) - Js'(t) dt 



J*t dt 



■ 2t 2 + s fl meters, 0 - t - 120. 

Notice that s^sjs s(0), the so-called initial distance* To assign a 
value to it we assbmeAhat distance ,is measured up from the launching 
pad with s (0) * 0. Accordingly, 

ft s(t) - 2t 2 , 0 * t * 120, 

# 

and the rocket's motion is completely described for the first two 
minutes of f 1 ight . 

We will take up what. happens to the rocket 'Sfter 'burn- 
out when we getTTcTthe next exercises. For the moment, let 
us look again at the decisions and assumptions we have made/ 
and how they enabled us to calculate the rocket's height as 
*a futicli^jhjwf' time. ^ ** 

# We first decided on a notation s(t) for the rocket's 
height as a function of time. In terms of this notation we 
wrote s'(t) for the velocity and s ,f (t) for the acceleration. 
Then we made assumptions about how the rocket worked and 
how it was launched, and decisions about how time 'and dis- 
tance were to be measured.. These translated into numerical 

\ 
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data about s(t) and \ts derivatives, and lead to, a descr 
tion of the _motion_.during._thl "burn" period. 



ASSUMPTIONS 
DECISIONS 



GENERATED 
DATA - 



CONCLUSIONS 
ABOUT THE MOTION 
DURING THE BURN PERIOD 



Time is 
measured in 
seconds and 
s(t) in meters. 
The en g i ne is 
on for 
0 < t <. 120. 
The engine 
gives full' 
thrust while 
on. 

The rocket 
Is fired 
from rest. 



Distance is 
measured up 
f rom^the 
launching ' 
pad. 



s"(t) = k 
0 < t < 120 



s 0 - 0 



s!(t) - kt + 



0 £ t £ 120 



s* (t) = kt + 0 » kt 
s(t) *» ly 2 + s 0 
0 £ t < 1^0 



s(t) « 2t 2 / 0 » 2t 2 

j 

0 </t '< 120 



— Ar 



^ 'w^npw_use the* equations for s(t)' and s' (t).to answer, the 
questions with which we began this section. The" table that 
follows" .shows how the fl^t two cftesYlons can be' rephrased In 
terras of s(t) and s'(t)<and then answered with simple Calcula- 



tions. Do'the same for the remaining questions. 
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REPHRASED IN 

QUESTION » TERMS OF THE ANSWERED 

MODEL * 



T 



How high is the rocket s (60) » 2 C60) 2 it 

1 minute after launch? ■ s(60) » ? - » 720O m 

« 7.2 km 



How fast is it climb- 
ing one minute after / ^ 

1aunch? , s'(60) - k (60) m/sec 

a) in m/sec $'(60) » ? -2*0 m/sec 

b) in km/h » 864 km/h 



125. How high will tfce 
. rocket be when the 

engine stops? 
^a) in meters 
b) in kilometers 

126. How fast will it 
s9 ^ be climbing when 

the engine stops? 

z 

127. When will the 
rocket be 20 ki^F .For what t is 



s T - - 

above the launch s(t) * 20,000 m? 

site? 



128. How long does it 
take the rocket 

to reach a velocity * y 

of lOOm/sec? 



123. How loncj did It 
take the rocket 
to rise the 
first SO m? 
Can«a goojd 
runner run >S0 m 
that ,fast? 



130. Hqw long di # d It 
take the rocket 
' to travel the 
next SO m? 
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As you gain -experience, you will not write as 
much as we did when we discussed the rocket problem. 
You might not need to write more than t 

s n (t) = 4 m/sec 2 

s 1 (t) = 4t + v 0 » 4t m/sec 

s (t) = 2t 2 + s 0 = 2t 2 m 
Remember to write d^rn th^jun^t&^-though, if only fori" 
later reference. i " ' -^ v * ^ 

It was' conveniient^to have ha'd the init ial^Wlocity 
v 0 and the initial distance s 0 both equal. to zero. -This 
allowed us to describe the velocity and distance traveled 
by the rocket during the burn period simply, by the 
equations s'(t) * 4lj and s(t) = 2t 2 . It will not always s 
be possible to'make-Vo and s 0 both zero 'in describing a « 
motion, however, norj will making them zero always be 
desirable. Exercise! 131 below- is a case in point. * " 
}. : s 

EXERCISES , / 

r 

There is more tj be learned about the flight of the rocket.' 
The rocket coasts upvjards for a while after the engine shuts off. 
For how long? And hjjbw high? 

To answer the/se questions we need a mathematical model that 

is different fronythe one we have been using. The reason for l^his 

is that when theyengine shuts down the force acting pn the rocket 

changes. The acceleration of the rocket is no longer k m/sec 2 

upwards provided by the engine, rather 9.8 mVsec 2 'downwards , ' 1 - v 

/ t - - » 

provided by .the. earth's .gravitational attraction. * * s 

Jf we continue to measure distances* 1 as before, tfut reset our 
. . / . ' ^ * a « a t 

clock toystart with t * D again, -bhe equation^for the~ acceleration 

y • ' * - »*,*•*' 3 n?"^ 

becomes/ * ■* * » * ' 4 X 

/ 7 /- * * v 

/ s»»(t) - -9.8 m/sec 2 * . 0 <° t ■ ' W ft** - 

i4at i ■ * 



so 



s» (t) = 9.8t + v 0 m/sec 0 < t 

s(t) « -4.9t 2 + v 0 t + s 0 m 0 < t 
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131. a) What numerical values should v^ and s 0 have? 

ib) Rewrite the equations just given for s'(t) and s(t) using 
the Initial- values from Part a. • 

Now complete the following table. 

* * • 
< •'. REPHRASED IN* 

QUESTION TERMS OF THE ANSWERED 

/ MODEL 



132. 


How long does the 


> 




rocket coast upwards * , 






after burnout? 


*- 


• 

133. 


How high does the ' 






rocket go? , % 




134. 


When do the equations 






of motion predict the 






rocket wi 1 1 crash? 





135. What is the rocket* s 
predicted speed 
just before it 
crashes? 



136. Would you expect a real 
rocket to behave as 
predi€t(M in Exercises 
134 and 135? Explain. 



V 



* * 9. REPEATED INTEGRATION 

As you saw in the preceding section, when we 
integrate more than once to solve a problem we need a 
corresponding number of fnitiaU conditions to determine 
Vthe constants of integration.. Here ,are two more exam- 
ples. 

* * 
EXAMPLE 13 . Fin/I f(x) if f"(x) = 12x - V (0)* - 5 and f(0) » -3? 

SOLUTION The initial conditions can be used one at a time as 

. >. 
* ' 23 
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V 

f '(x) ->2x - 14 



- ' the Integration proceeds, or 'at the end of the inte- 
gration. 
METHOD #1 

1 . Integrate f"(x) 
to get f» (x) - 6x 2 - Hx + C r 

,2. USe the condition f'(0) - 5 to find Cj. 

6(0) - 1M0) + Cj - 5 

Ci - 5 

3. Integrate f 1 (x) - 6x 2 - l4x + 5 

to get f(x) « 2x 3 - 7x 2 + 5x + C 2 . 

k. Use the condition f(0) = -3 to find C 2 . Jfa 
2(0) 3 - 7(0) 2 + 5(0) + C 2 - -3 
C 2 --3. 

5. Write the completed formula for f(x). 

<?M - 2x 3 - 7x 2 + 5x - 3- 

METHOD #2 

1 1. Integrate twice. f"(x) = 12x - 14 

f ' (x) - 6x 2 - 14x + C 
f (x) - 2x 3 - 7x 2 + C x x + C 

2. Substitute x * 0 

in the last two Cj « f*(0) * 5 
equations' to find C 2 » f(0) * -3/ 
x Cj and C 2 . 

3. Write the formula 

forf(x), • f(x)-2x 3 -7x 2 + 5xr 3 

' EXAMPLE 14. Find f (x) if f"(x) » Ax - 2 and if the graph of f 
passes through the point (1,0) with slope 3. 

SOLUTION • s 1. Integrate f H (x) ^ 

fc . * f"(x) « kx - 2 

twice. 

* f ' (x) - 2x 2 - 2x + C t 

f (x) - j x 3 -- x 2 + + C 2 



24 
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2. Determine C l from the 
fact that the slope of f ! (l) ■ 3 
the graph is 3 when 2(l) 2 - 2(1) + C x - 3 
X - 1. 

3. Substitute this* value 
of Cj in the expres- f (x) - j x 3 - x 2 + 3x + C 2 »' 
si on for f(x). 

4. Determine' C 2 from the 

fact^at f(1) - J>. 2 M , 3 „\ 2 - t " ° 

\ j(l) 3 - (I) 2 + 3(0 + C 2 - 0 



ar 

3 

C 



v + C 2 - 0 



5. WriteVf(x). f( x ) - | x 3 - x 2 + 3x"- $ 

* • 3. 



140. 



k M (y) » ^y; k' (9) * 10, k(0) = 8. 



Ul. >"'(t) = 6; p»(o) » -8, p'(0) » 0, p(0) -5. - 

^2. f M (t) » 1 - 6t* and the graph of f passes through the 
point (2,0) with slope 0.. 




143, 



rJ'(x) Vs » y-, and the graph of r passes through the point 
(4,4) with slope 3, 



144. g"(x) » e x , and the graph of g passes through the origin 
w*th slope 2. 

145. . h M (x) » - x > t), and the graph of h passes through 



/ (1,2) with slope 0. 

i 



EXER^SES 

Find thej function determined by each set of conditions. 

137. fjW - 2-6x; f (0) - 4, f(0) « 1, 

138. g"(t) - 30t; g'(l) - 0, g(l) » 10. 

139. h"(x) » e X ; h'(0) » h(0) r 1. 



25 
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Section 2 

i 

6. kx + C 7. "25x + C 

8. 17 * 9. -3 



10. | + c * ■ n. f + c 



12'. C •* * - 1.3. In |x| + 

14. *x 7 1S. + C 

16. -xJ 03 + C 17. x 16 * 

18. e x , 19. Ax 3 

20. x 3 ' ' 21. -103.5 



22. kx + C 23. ax + C 

24. is. -b 

26. 32t + C * ' '27. a 



28. I + C . ' 29> ?3B 2 



3 

30. s 5 + C . 3K* + d 



7 



32. % e y + r 33. £ - 

34. t 0 t + C % 35. % 32t 0 

c 

36. a 0 t + C ^ * 37. v 0 t + C 



38. yo m ♦39. 



s 



0 



Section 3 







x + C 




l, 

MX - 




46. 


* - 1 


A?/ 


_ X 

X 




48. 


-e x + C 


49. 

• 


C 


• 


o 50. 


x 2 

m -j- + bx + C 


5'- 


y3 „ 



52 - 25 - 3s 2 53.-^«.'l e x + c 



,2 



54. 3 1n|z| + Ae 2 + C 55. 3x 2 -,2x + 7 

Section Jf 



56. | t* + C 



" 2 # * 
58. j x 2 + C 

60. -8s^+ C 



62. - + c 



64. ^ 
66. -I!- 

>t* 

S * 



68. 



4 



70-. in|xf+ c , 

72. 1 +1 ■ 

x 

74. ^-i + c ' 
3 x 

76. e 2 + i " ^ 
z 



^ 78. 200 ln|z| + C* 

1 i 



80. y* + C 



4 



4 



ERLC 



,7 •• -1 

57. x or jff" 

* 


» 












» 


65. - -V • 




67 ' 8 - 




69. 1 + C 
* v 




71. -2 ln|x|-+j c - 




73 x 9 \ -,i T . 

X 2 " * * 

75- y t + ln|x) + C 




77.. e* + 1+ c. 

'». fi.c :/ 


r 




i 

29 

i 
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Section 5 



J 81 . 


True 






Fal se 


83. 


True 




84. 


False 


L 


ra 1 se 




oo • 


True 


87. 








32t 


89- 


1 






32t + v 0 


91. 


1 
1 

y 




92. 


32 


93. 


m 




94. 


2^? " 1 


. 95. 


600 
- yy- 




96. 


23 e x 


97. 


-240x* -> 70 




98.. 


1 

3/v 



INDEFINITE INITIAL VALUE PARTICULAR 

DERIVATIVE INTEGRAL CONDITION OF C SOLUTION 



f»(x) Jf'(x)dx f(x 0 ) 



f(x) 



2x ' x 2 + C ♦ ; fT-£) - 29* 4 fM = x 2 + k 



§9 % - 1n|x| + C f(e) - -3 -* f(x) » ln|x|- A 

100, . ~ ln|x[ + C f(l)*- 2 2 F(x) = In |xj + k 

* '* * ^ 2 • 2 * 

;!pl\ r '--x **T~ +C , f(D - 0 i f(x) 

)$2-j x*. +'6 ^ |i+6x + C- f(1) - 10 y ' f(x) - j-+ 6x + j- 



103. 32 32x + C f(0) - 0 0 f(x) - 32x 

• **> * * 30 
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v 0 * + 


f(o) « 0 


0 f(x) » v 0 x + 5 


105* 


,32x + v o 


l6x 2 + v fl x + C 


f (0) * 0 


0 f(x) = 16x 2 +*v 0 x 


106. 


U(x) - - 


kx 2 
2 , 


. 


* > 












107. 


0 ■ * 


Section 7 

108. 


0 


109. 0 


, ( 110. 


k 


1 1 1 


7 


112. 16 


*> 


FORMULA 
FOR f *(x) 


Jf (x)dx 


-A- POINT ON- 
THE GRAPH 
OF F 


FORMULA 
FOR f(x) * 


: 


2x 


< x 2 + C " 


(5, JO) 


x 2 - 5 

* * 


113. 


5 


5x + C 


(-2, 1) 


5x + 11 ' 

* 


* 

n't. 


8x 


Ax 2 + C 


(0, /?) * 


kx 2 + 

X 


115.' 


-Ax + 3 


-2x 2 + 3x + C 


(-1. 1) 


-2x 2 + 3x + 6 



116. "9.8x " i,.9x 2 + C (1,3)-* 4.9x 2 * 1.9 



117. 


e x - 2 


e* - 2x + C 


■ (0, 7) 


e x - 2x + 6 


118*. 


6 


6 ?n|x| + C 


(1. 4) 


6 ln|x| + 4 




3 

X 7 " 


X 


(2, 0) 
1- 


— , * 

M ..' 


120. 




fx* + C 


(3, 0) 


fx*-**,' 


* 






* 

♦ 
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QUESTION 



How high Is the 
rocket 1 minute 
after launch? 



How fast is It 
cl imbing one min- 
ute after Jaunch? 

a) in meters <*f 

b) in ki lometers 



REPHRASED IN 
TERMS OF THE 
MODEL 



s(60) 



,s'(60) - ? 



ANSWERED 



s(60) - 2(60) 2 m 
» 7200 m 
- 7-2 'km 



s'(60) « 4(60) m/sec 

a) » 2k0 m/sec 

b) = 864 km/h 



125. How high wi 1 1 the 
- rocket be when 
the engine stops? 

a) in meters 

b) in ki lometers 



s(120) 



sXl20) 
a) 
b) 



2*(l£0) 2 m 
28,800 m 
28.8 km 



126. 



How fast *ill it be 
cl Imbing when the 
engine stops? (The 
speed of sound in 
air at sea level is 
about 335 m/sec. 
Speeds of about 500 
m/sec are typical 
for oxygen molecules 
at room temperature. 



1^7. When will the rocket 
* *be 20^km above the 
launch site? 



'(120) 



For what t is- 
s(t) » 2^000 m? 



(120) 
a) 
b) 



Ml 20) m/sec 
480 m/sec 
1728 km/h 



t 2 ». 10,QOO . 
100 sec 



128. How long does ft 

take the rocket to 
reach a velocity 
of 100 m/sec? 



For what t is 

s 1 (t) = 100 m/sec? 



4t « 100 
t ■ 25 sec 



-129. How long did it take 
the roqket to rise 
. the first 50 m? Can 
a good runner run 
50 m that fast? 



For what t is 
s(t) - 50 m? 



2t 2 
t 2 
t 

The 
brd 
seconds. 



« 50 " * 1 

- 25 
• 5 sec 

current world rec- 
for lOO.m is 9.9 
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QUESTION 



REPHRASED IN 
TERMS OF THE 
MODEL 



ANSWERED 



130. ; How long ( d id it 
. ^ake the rocket 
to travel the 
ne*t 50 m? 

* • • • 



Find ffie t for which a) 2t 2 ■ 100 



s(t) = 100 m. Then 
subtract 5 sec. 



^ t z - 50 

. t * 5/2 sec 

b) 5/2" - 5"" 5(J?-\) 

* 5(0. 414) 

• = 2.07 sec 



1)1. a) v 0 „* 480 m/sec, and s 0 * 28,800 m 

•) «'(t) » -9.8t + ,480 m/sec 
^ s(t) = r4.9t 2 +-"480t + 28,800 m 



xercises 126 a 



and 125). 



135- 



r 

QUEST 1 Ok 


REPHRASED IN 
TERMS OF TH£^ 
MODEL 


p 

ANSWERED „ 


How long does the 
rocket, coast upwards 
after burnout? 

* 


Fo t r what t is 
s'(t) = 0? 


0 = -9.8t + 480 
- - t = 480/9.8 sec 
= 49 sec * 


How high^loes the 
jacket gcH 
4 ' 


s(49) «•? 


, ' s(49) * 40,555 ra 


When do the equations 
of motion predict the 
rocket wij 1 crash? 


'For what t is 
*(t) * 0? 


J> ^ .1 40 sec 
^Wtfter burnout. ■ 
^ • » 


* * 

What is the rocket's 
predicted speecf just 
before^it cra$hes? 


s»(l40) * ? 


-892 m/sec, or 

892 m/sec downwards 

(=3,21 l w . 2 km/K) 


« H 

No. The predictions are sure to he underestimates, they neglect 
the air resistance that wilbslow the rocket ' s fal l y assume the 
rocket does not fragment /and so forth. ^ _ 
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Section 9 



137. f(x) = -x 3 + x 2 * AX + 

138. g(t) - 5t! : 15t '+ 20 

139. h(x) - e x 

li,0. k(y) -ig.yl - 8y ♦ 8 

1*1- P(t) = t 3 - At 2 + "5 

'*2. f(t) t ' ♦ lOt - 

1*3. r(s) - 

1**..- 9(x) "- e x + x - 1 ' 

1*5. h(x) - ln|x| - x + 3 ' 




